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Preface 


This report provides a summary of results for the past year f s 
effort under contract NAS9-14970. The contract called for work on a 
x?ide variety of separate and distinct but related tasks* Below is a 
list of the tasks as contained in the original Work Statement* 

As a result of this contract a large volume of results has been 
generated* Technical and research reports previously submitted or 
presently being published are listed in Table I* 

Because of the diversity present in the task list for the contract, 
each major subdivision of this report has been written to be relatively 
self-contained* We hope this will facilitate use of the report by 
readers with different interests* 

Many of the tasks below are continuing; as a result, the discussion 
presented relative to them constitutes a progress report* Task 2.1 and 

2.2 are not continuing, thus the discussion relative to them is more 
complete. 

The various tasks have been managed by various Purdue staff members 
during the year. It is appropriate that the contributions of these 
people be recognized. 

Task Purdue Task Manager 


2.1 Test of Boundary Finding/Per Field 
Classification (ECHO) 

2*2 Stratification of Scene Characteristics 

2.2a Stratification by Machine Clustering 
2.2b Digitization and Registration of 
Ancillary Data 

2.2c Crop Inventory Using Full-Frame 
Classification 

2.3 LACIE Field Measurements 

2.4 Scanner System Parameter Selection 

2.5 Technology Transfer 

2.6 Large Area Crop Inventory Design 

2.7 Forestry Applications Project 

2.8 Regional Applications Project 

2.9 Interpretation of Thermal Band Data 

2.10 Super Site Data Management 

2.11 Soil Classification and Survey 

2.12 Improved Analysis Techniques for 
Multitemporal Data 


P . H . Swain 
P.H. Swain 


M.E. 

Bauer 

P.A. 

Anuta 

J.C. 

Lindenlaub 

M.E. 

Bauer 

R.P. 

Mroczynski 

R. A. 

Weismiller 

L.F. 

Silva 

L.F. 

Silva 

M.F. 

Baumgardner 

P.H. 

Swain 


The efforts of Dr. J.D. Erickson, the contract Technical 
Monitor and Mr. M.C. Trichel are especially to be noted and 
greatf ully acknowledged . 
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Table I. Technical and Research Reports 


052977 

042777 

120776 

091576 

070676 

111076 

110976 

090776 

082776 

062176 

052576 

051576 

*032576 

*031276 

052075 


Davis. The Focus Series: A Collection of Single-Concept Remote 

Sensing Educational Materials. 

Russell and Lindenlaub. Disseminating Technological Information 
on Remote Sensing to Potential Users. 

Vanderbilt, Silva and Bauer. A Laser Technique for Characterizing 
the Geometry of Plant Canopies. 

DeWitt and Robinson. Description and Evaluation of a Bidirectional 
Reflectance Factor Reflectometer . 

Kristof, Russell, Cary, Lube and Weismiller. Determining Land Use 
Patterns Through Man-Machine Analysis of LANDSAT Data - A Tutorial 
Simulation. 

Russell. Systematically Disseminating Technological Information 
to Potential Users. 

Landgrebe, Simmons and Biehl. An Empirical Study of Scanner 
System Parameters. 

Svedlow, McGillem, and Anuta. Analytical and Experimental Design 
and Analysis of an Optimal Processor for Image Registration. 

Montgomery, Baumgardner and Weisnuller. An Investigation of the 
Relationship Between Spectral Reflectance and the Chemical, Physical, 
and Genetic Characteristics of Soils. 

Wiersma and Landgrebe. The Use of Spatial Characteristics for 
the Improvement of Multispectral Classification of Remotely 
Sensed Data. 

Lindenlaub and Lube. Matrix of Education of Training Materials 
in Remote Sensing. 

Bauer and Davis. Stratification of Landsat Data by Clustering. 

Todd, Mausel, and Baumgardner. Urban Land Use Monitoring from 
Computer- Implemented Processing of Airborne Multispectral Data. 

Mausel, Todd and Baumgardner. An Analysis of Metropolitan Land- 
Use by Machine Processing of Earth Resources Technology Satellite 
Data. 

DeWitt and Robinson. Description and Operation of a Field Rated 
ERTS-Band Transmissometer . 


* Papers reporting earlier SRT work which has just appeared in print. 
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OTHER REPORTS 
Biehl and Simmons. 

Field Measurements Data Library Catalog, Volume I, 1974-75 Crop. 
October 29, 1976. 

Biehl and Simmons. 

Field Measurements Data Library Catalog, Volume II, 1975-76 Crop. 
December 1, 1976. 
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2.1 Test of Boundry Finding/Per field Classification 

INTRODUCTION 

One might wish to adopt a new classification technique for one of 
three reasons: 

• It provides greater accuracy than the currently Implemented technique. . 

■ The new technique is less expensive than the current technique. 

• The classification results produced by the new technique are easier 
to use and interpret than those currently produced. 

The ECHO* classification algorithms developed at the Laboratory for 
Applications of Remote Sensing (LARS) have shown promise in all three areas 
listed above. * 

Contemporary classifiers for analysis of remotely sensed data compare 
the spectral measurements from each feature of each point to class statis- 
tics, computing a likelihood of discriminant function associated with each 
class, and categorizing the point according to the class with the largest 
discriminant function value. Each point is classified individually on the 
basis of its spectral measurements alone. One premise of this technique 
is that the objects of interest are large in comparison to the size of a 
point. If this were not so, a large proportion of poinfs^would be compo- 
sites of several classes, making statistical pattern classification 
unreliable since pre-specif ied categories would be Inadequate to describe 
actual states of nature. From this premise it follows that objects are 
represented by arrays of points, and that a statistical dependence exists 
between consecutive points. Contemporary classifiers fail to exploit 
the statistical dependence between adjacent points when assigning classes. 


*ECHO stands for Extraction and Classification of Homogeneous Objects 
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Giving no consideration to the spectral response of adjacent points 

contributes to a "salt and pepper" effect in classification results. While 

the results produced may be accurate, they may be too complex, and can 

even be confusing when attempting to produce a useful analysis product 

2 

such as a forest inventory map. Consequently time and effort have been 

invested in developing smoothing programs which "clean-up" the "salt and 

2 

pepper" effects in classification results. 

The ECHO processors benefit from spatial information by first aggre- 
gating into groups points whose spectral responses are not significantly 
different in a statistical sense, and then applying a maximum likelihood 
classification rule to these homogeneous groups. Homogeneous objects 
are identified in a three step process. First, cells are formed by 
systematically partitioning the data into N by N sized blocks of pixels. 

The statistics of each cell are then compared to a homogeneity criterion^ 

Points which do not comprise homogeneous groups are classified on a point- 
by-point basis, just as contemporary classifiers classify all points. 

Statistics of adjoining homogeneous cells are then compared to annexation 
thresholds. Adjoining cells which appear to belong to the same statistical 
population are combined into a single object. 

Two separate ECHO algorithms have been developed. The first, Supervised 
ECHO, makes use of pre-specif ied class statistics to identify homogeneous 
objects. The second, Nonsupervised ECHO, identifies homogeneous objects 
without the use of class statistics. Consequently, those objects identi- 
fied by the Nonsupervised algorithm may be used to aid in the training 
process. 

Past Work at LARS 

Much of the background research on the ECHO concept was performed at 

3 

LARS during 1975 and is documented in the Final Report for 1975, R. L. Retting' s 
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doctoral thesis, a LARS Information Note, and in symposium proceedings. 

The partitioning of a data set into homogeneous objects can occur 
in one of two general approaches. One can utilize a disjunctive approach, 
starting with large areas and repeatedly subdividing them until each unit 
appears to be homogeneous. Alternatively, one can choose a conjunctive 
approach where adjacent units which appear to belong to a single object 
are combined. One advantage to the conjunctive approach is, if second 
order statistics are to be used, the second order statistics of the combined 
area can be calculated from the second order statistics of each constituent 
unit. The statistics for each subdivision of the initial partition made 
under the disjunctive approach must be recalculated starting with the 
individual pixel values, however. Hence the conjunctive approach is poten- 
tially faster than the disjunctive approach. 

Prior to the 1975 work, Rodd's conjunctive partitioning algorithm** was 

combined with a minimum distance sample classifier and an improvement in 

classification accuracy over conventional classification was observed, however 

4 7 

processing time was increased. Gupta and Wintz, added a second order 
statistics test to Rodd's first order test, but obtained essentially the 

g 

same results at an increase in processing time. Robertson implemented a 
disjunctive partitioning algorithm with a minimum distance classifier which 
acheived the same classification accuracy as the conventional point-by- 
point classifier with an order of magnitude increase in computer time. 

During 1975 new statistical criteria were applied to the partitioning 
algorithm producing a "Supervised" field extraction routine to add to the 
Nonsupervised routine. In addition, the sample classifier utilized by ECHO 
was altered from a minimum distance strategy to a maximum likelihood 
strategy. Experimental results were obtained for two aircraft and two 
LANDSAT-1 data sets which indicated that ECHO algorithms improved classifi- 
cation accuracies for the four tested areas and that they reduced the CPU 
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time required to perform the classifications. 

FORTRAN programs were produced implementing the Supervised and 
Nonsupervised algorithms in two phases. Channel selection, cell size, 
and cell homogeneity criteria were specified in the first phase of the 
Supervised processor, annexation in the second. For the Nonsupervised 
processor, only channel selection and cell size were specified in phase 
one, cell homogeneity and annexation thresholds were specified in phase 
two. The Nonsupervised processor did not have the capability to "split" 
cells; i.e. the constituent points in a singular cell had to be classified 
as a small sample rather than point-by-point. 

Documentation of these programs consisted of commented program listings 

3 

and a thesis. 

Objectives of the 1976-1977 Task 

The objectives of the FY77 work were to: 

1) Deliver to NASA documented FORTRAN programs implementing the ECHO 
algorithms, and to make improvements to those algorithms, where appro- 
priate. 

2) Systematically test the algorithms on MSS data for agricultural 
regions as observed by the LANDSAT satellites, aircraft scanners, and on 
the simulated Thematic Mapper data. 

3) Provide products enabling the determination of the utility of the 
object maps to a LACIE Analyst Interpreter in the selection and labeling 
of training fields. 
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DESCRIPTION OF WORK 
A. Documented Fortran Programs 
1 . Supervised Processor 

FORTRAN listings and program abstracts for the Supervised ECHO 
processor are presented in Appendix A. This processor requires class 
statistics for field extraction as well as for classification. The 
general flow of this algorithm is presented in Figure 1. The software is 
designed to function in either a single phase or a two phase mode. The 
two phase approach accomplishes channel selection, cell size, and cell 
homogeneity criteria in phase one; annexation tests take place during 
phase two. The intermediate results between phases are stored on magnetic 
tapes. The final .results are written either on disk or on tape. Func- 
tioning in the two phase mode allows the analyst to vary the annexation 
threshold and, thus, to produce several results outputs without repeating 
the expensive cell homogeneity test. 

The single step approach allows the analyst to specify all the neces- 
sary input parameters on a single control card deck and produce classifi- 
cation results without the necessity of an intermediate tape. When the 
analyst wishes to produce classifications with varying annexation thresholds 
for a particular area, it is less expensive in terms of computer time to 
utilize the two phase approach. When a single annexation threshold is 
desired, it is more efficient to use the single phase approach. The 
Supervised ECHO processor determines whether a single phase classification, 
the first phase of a two phase classification, or the second phase of a 
two phase classification is to be generated based on the control cards 
specified. 



FIGURE 1 

SUPERVISED ECHO FLOW 
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Input and output requirements for each processing mode of the 
Supervised processor are summarized in Table 2 . 1-1 and Table 2 . 1-2 . 
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TABLE 2.1-1 
Supervised ECHO 
Single Phase Mode 


Input 

Cell Width Parameter 
Cell Homogeneity Threshold 
Annexation Threshold 
Channel Selection 
Specification of the Area 
to Classify 

Multispectral Image Storage 
Tape (MIST) 

Specification of Results 
Destination 


Output 

Results File 

Classification Map (optional) 
Singular Cell Map (optional) 


Class Statistics 
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TABLE 2.1-2 
Two Phase Mode 
Supervised ECHO 

Input Output 

Phase One 

Cell Width Parameter Intermediate Results Tape 

Cell Homogeneity Threshold Singular Cell Map (optional) 

Channel Selection 
Specification of the Area 
to Classify 

Multispectral Image Storage 
Tape 

Intermediate Tape and File 
Specifications 
Class Statistics 

Phase Two 

Intermediate Results Tape 
File 

Annexation Threshold 
Specification of Results 


Results File 

Classification Map (optional) 


Destination 
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Figure 2 

Effects of Cell Homogeneity Threshold 



Perpoint 


Low Homogeneity Threshold 


87.5% 

7.3% 

1200 Seconds 
34.8% 


Field Center Pixel Performance 
RMS Proportion Error 
CPU Time 

Classification Variability 


90.0% 

2.5% 

754 Seconds 
29.9% 



91.4% 

2,. 5% 

494 Seconds 
22.5% 


Field Center Pixel Performance 
RMS Proportion Error 
CPU Time 

Classification Variability 


90.7% 

2 . 8 % 

370 Seconds 
18.6% 
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Photographs of classification results of aircraft data collected over 
Williams County, Kansas. Notice the "salt and pepper’* in the perpoint 
results (above). The ECHO results (below) have both higher classifi- 
cation performance and lower classification variability. The ECHO 
classification was produced using 40% of the CPU time required to 
produce the perpoint classification. 






2. Nonsupervised Processor 

Output products for this task included those allowing the determi- 
nation of the utility of the resulting object map to the LACIE Analyst 
Interpreter for identifying and labeling training fields. Since the 
Supervised processor utilizes class statistics for the identification of 
homogeneous objects, it has only limited value as a training aid. The 
Nonsupervised processor, on the other hand, does not require class sta- 
tistics for object extraction, and thus, may prove quite beneficial as 
a training aid, serving, in a sense, as a spatial clustering algorithm. 

The Nonsupervised processor utilizes a homogeneity test which com- 
pares o^/lJLj to a threshold t^; where i refers to feature or channel 
and j refers to the cell. If the standard deviation divided by the mean 
for feature i exceeds the user-specified threshold t^, the cell is "singular” 
and elements of the cell will be classified on a point-by-point basis. 

Cells are annexed to fields on the basis of a two-step test, 1) that the 
channel variances of the f ield are equivalent to the channel variances of 
the cell and 2) that the channel means of the field are equivalent to the 
channel means of the cell. Should either of these criteria not be met 
for any channel, annexation will not take place. 

The research software was designed to run in two phases. The first 
phase merely calculated the cell mean and covariance matrices and wrote them 
on tape; the second phase proceeded to perform field extraction followed 
by classification. This process has two disadvantages. First, although 
an object map could be produced, it had to be produced in phase two, the 
same phase that required a class statistics deck for input in order to 
perform the classification. Second, since only cell mean and covariance 
matrices were written on the intermediate tape, when singular cells were 
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identified in phase two, they had to be classified as small samples, a 
cell at a time, rather than on a point— by-point basis, because data values 

i 

for individual points were not available to the phase two classifier. 

The software was restructured to move the field extraction algorithm 
into phase one of the processor and produce an intermediate tape which 
contains : 

1. class means and covariance matrices for each homogeneous object 
identified. 

2. an object map containing the mean for channel i of object j in 
every pixel of object j and the original data values for those 
pixels belonging to singular cells. 

3. a tag array identifying whether a pixel falls in a singular cell 
(and should be classified individually), or in a homogeneous object 
(and hence should receive a class assignment based on the sample 
classification of the object to which it belongs). 

The second phase of the restructured Nonsupervised processor reads 
the intermediate tape and the class statistics deck, performs a maximum 
likelihood sample classification (utilizing second order statistics where 
possible) on the objects identified, and a maximum likelihood point-by- 
point classification on points falling in singular cells. Figure 4A 
presents the general processing flow for phase one (field extraction) of 
the Nonsupervised processor. The general processing flow for phase two 
(classification) is presented in Figure 4B. Inputs to outputs from the 
Nonsupervised ECHO processor’s two phases are described in Table 2.1-3. 



FIGURE 4A 
GENERAL FLOW OF 
NONSUPERVISED ECHO 
PHASE 1 
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TABLE 2.1-3 
NONSUPERVISED ECHO 
PROCESSOR 

INPUT OUTPUT 

Phase One : 

Channel Selection Intermediate Tape containing 

Cell Width object map, object statistics 

Cell Homogeneity Thresholds and pixel tag array. 

(a 1 /y i for channel i) 

Annexation Thresholds 
(mean and variance) 

Multispectral Image Storage 
Tape 

Intermediate Results Tape 
and File Specification 

Phase Two : 

Class Statistics Results File 

Intermediate Tape from 
Phase One 

Specification of Results 


Tape and File 
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B. Tests of the ECHO algorithms 
1 . Data Sets 

The second obejctive of Task 2.1 is to test the ECHO algorithms 
on MSS data for agricultural regions. Data sets are to include LANDSAT, 
aircraft, and simulated Thematic Mapper data. Ten LANDSAT, three aircraft 
and eight simulated Thematic Mapper (2 sites and 4 resolutions) data sets 
were selected for the analysis. One bitemporal data set was examined 
separately. These data sets are summarized in Table 2.1-4 
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TABLE 2.1-4 




Data Sets 




(LANDSAT) 



CHANNEL 

NUMBER 



WAVE BAND 
yM 


1,5,9 



.5-. 6 


2,6,10 



.6-. 7 


3,7,11 



.7-. 8 


4,8,12 



.8-1.1 


Area 

Channels Used 

Data 

Collected 

Data Set 

Graham County, Kansas 

9, 

10, 11, 12 

5/26/74 

LACIE/SRS 

Grant County, Kansas 

5, 

6, 7, 8 

5/19/74 

LACIE/SRS 

Haskell County, Kansas 

9, 

10, 11, 12 

5/27/74 

LACIE/SRS 

Kearny County, Kansas 

9, 

10, 11, 12 

5/27/74 

LACIE/SRS 

Huntington County, Indiana 

1, 

2, 3, 4, 

7/7/73 

CITARS 

Shelby County, Indiana 

1, 

2, 3, 4 

9/7/73 

CITARS 

White County, Indiana 

1, 

2, 3, 4 

8/21/73 

CITARS 

Livingston County, Illinois 

1, 

2, 3, 4 

7/16/73 

CITARS 

Fayette County, Illinois 

1, 

2, 3, 4 

8/21/73 

CITARS 

Lee County, Illinois 

1, 

2, 3, 4 

8/5/73 

CITARS 
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TABLE 2.1-4 (Continued) 
Simulated Thematic Mapper 


TM CHANNEL 
NUMBER 






WAVE BAND 
yM 


1 






.45-. 52 


2 






.52-. 60 


3 






.63-. 69 


4 






.74-. 80 


5 






.80-. 91 


6 ' 






1.55-1.75 


7 






10.4-12.5 


8 






.74-. 91 


Area 

Channels 

Used 

Date 

Resolution 

Williams County, ND 

1, 

2, 

3, 

6, 

8 

8/15/75 

30m 

Williams County, ND 

1, 

2, 

3, 

6, 

8 

8/15/75 

40m 

Williams County, ND 

1, 

2, 

3, 

6> 

8 

8/15/75 

50m 

Williams County, ND 

1, 

2, 

3, 

6, 

8 

8/15/75 

60m 

Finney County, Kansas 

2, 

3, 

5, 

6, 

7 

7/6/75 

30m 

Finney County, Kansas 

2, 

3, 

5, 

6, 

7 

7/6/75 

40m 

Finney County, Kansas 

2, 

3, 

5> 

6, 

7 

7/6/75 

50m 

Finney County, Kansas 

2, 

3, 

5, 

6, 

7 

7/6/75 

60m 
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TABLE 2.1-4 (Continued) 

Aircraft 

Wavebands Used 


Area 

(VM) 

Data Collected 

Williams County, ND 

.4-. 49, .59-. 64, .65-. 69, 

8/15/75 

Finney County, Kansas 

.82-. 88, 1.53-1.62, 10.1-11.0 

7/6/75 

Tippecanoe County, IN 

.52-. 57, .61-. 70, 

8/13/71 


.72-. 92, 1.0-1. 4, 1.5-1. 8 



Bitemporal 


Area 

Wavebands Used 

for Each Date Dates 

Data Set 

Grant County, Kansas 

.5-. 6, .6-. 7, 5/9/74 

LACIE/SRS 


.7-. 8, .8-1.1 6/14/74 
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2 . Training Procedures 

a. LANDSAT Data Sets 

The LACIE training sets were created using ground truth information 
provided by JSC for both test and training fields. Odd-numbered fields 
appearing in each subclass were used for training; even-numbered fields 
were used for test. For Graham and Grant counties, statistics were generated 
by using the STATISTICS processor. For Haskell and Kearny counties sta- 
tistics were secured by clustering the training fields of each class into 
four subclasses which were then pooled on the basis of output from the 
SEPARABILITY processor. It should be noted that both training and test 
fields for the LACIE/SRS data are large enough to inset field boundaries 
approximately two pixels inside estimated field boundaries. This inset 
allows for any image misregistration which may occur between any two dates 
on the multitemporal input runs. Proportion estimates for the 1974 
LACIE/SRS segments were provided in ground truth packets provided by JSC. 

The CITARS training sets were originally created by a supervised 

procedure using ground truth information provided by the Agricultural 

Stabilization and Conservation Services (ASCS) to select both the training 

q 

and test fields. Those training sets were used without change in the 
ECHO tests. 

Six counties in the CITARS experiment were used as test sites. A 
data set free of clouds which occured late in the growing season was re- 
quired for each of the six counties. Dates from July, August or early 
September were selected. The classification results for all of the CITARS 
experiments are catalogued on a series of LARS tapes. The desired data 
sets were located on the catalogued CITARS tapes and the statistics which 
had been used for the CITARS experiment were obtained by using the LARSYS 
PUNCHSTATISTICS processor. The pooling of classes was determined by running 
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the PRINTRESULTS processor and requesting training field results. This 
combination of requests produces a table of available classes in the 
statistics deck and also the informational names under which they were 
classified. By using this list it was possible to reconstruct the com- 
binations of spectral classes and pooled spectral classes which were 
needed to reproduce the original classification. 

Training and test fields for each CITARS classification were secured 
and appropriate control cards added to evaluate the ECHO classifications 
of CITARS data sets. The proportion estimates used for the CITARS 
evaluation are estimates of the proportions of the various classes for the 
entire county made by the SRS. 

b. Simulated Thematic Mapper Data Sets 

Training sets from the simulated Thematic Mapper tests performed at 
LARS in 1976 were used for ECHO analysis. These training sets were generated 
by selecting fields of known cover types ana clustering each informational 
class separately to define subclasses. 

Color infrared photographic mosaic prints were made from photographic 
data collected concurrently with the scanner data. Informational class 
information provided by ground observations was transfered to clear plastic 
overlays on the mosaic print. The analyst could then easily locate the 
corresponding fields in his cluster maps and assign the field coordinates 
to the informational classes. 

Statistics were calculated for each training area and compared using 
the SEPARABILITY processor. Similar classes were combined, where indicated, 
and the data set was used to classify the flightline. Training areas were 
not excluded from the test fields since the test fields had been pre- 
selected for the entire flightline. 
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Two to four subclasses were found in each informational class. The 
Kansas flight was an exception. Because of severe line-to-line changes 
in signal level in the original 6 meter scanner data, the analyst was 
forced to create more spectral classes to account for the within-class 
variations due to excessive noise. This was most apparent in the 30 and 
40 meter resolution data. The effect was reduced but not eliminated in 
the 50 and 60 meter resolution data. Alternate fields were used for 
training and test decks, respectively. The procedure was repeated for each 
of the four resolution sizes. As resolution size increased, the number 
of spectral classes decreased. 

The entire training set selection procedure was repeated for each 
resolution size so that any effects on training set selection which might 
be caused by data resolution would be included in the analysis results. 

An example is the increasing difficulty and eventual impossibility of 
selecting samples from small, or narrow, fields as the resolution size 
increases. 

c. Aircraft 

Six Meter NASA Aircraft Scanner Data 

The two 6-meter aircraft scanner data sets used to generate the 
simulated Thematic Mapper data set were used. The same training and test 
fields used for the Thematic Mapper were available in the six meter data 
set. Because of the vary large number of data points in the full data 
set, only the first two miles (one third of the total flight line) were 
used for the ECHO evaluations. 

These data sets were not corrected for sun and scanner angle effects. 
To compensate for these angle effects, training fields were distributed 
across the width of the flight line. At least six fields in each infor- 
mational class were used in the training set. Fields from informational 
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classes were clustered together into five spectral classes. All spectral 
classes from all informational classes were combined into a statistics 
deck and appropriate pooling was done based on SEPARABILITY results. The 
classification results were strongly influenced by angle effects. 

Flightline 210, Corn Blight Watch Experiment 

Training and test fields from the Michigan Scanner data set were 
located on the LARS display system using ground information collected 
during the Corn Blight Watch Experiment. At least six fields were selected 
to represent each informational class. Each informational class was then 
clustered to produce five spectral subclasses. All spectral classes were 
combined and processed with SEPARABILITY to select the subset of the 
original 12 channels which was to be used for classification. A subset 
of five channels was chosen and the training set was again processed with 
SEPARABILITY to determine which classes should be pooled. 

d. Bitemporal Data Set 

A bitemporal data set from Grant County, Kansas was created by selecting 
data from a spring date and a mid-summer date. Training and test fields 
were available from the LACIE data and were edited to reduce the number 
of single line training fields by about 30 percent. Data from the two dates 
were clustered together within each informational class with four spectral 
classes requested. Classes were grouped as indicated by the CLUSTER 
algorithm and classes with fewer than 50 points were deleted. 
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3 . Dependent Variables 

There are six variables which were monitored to evaluate the ECHO 
algorithm: 

• CPU time 

• Field center pixel classification performance 

• Training field classification performance 

• Full field classification performance 

• RMS proportion estimate error 

• Classification variability 

These variables are related to the stated reasons for adopting a 
new classification technique listed in the introduction. The CPU time 
required to perform a classification is one way to measure the cost of 
classification. Field center pixel, full field, and training field per- 
formances and RMS proportion estimate error are all ways to evaluate the 
accuracy of the classifier. Classification variability is a measure of 
"salt and pepper effect" in classification results. 

The CPU time required to execute each of the ECHO classifications has 
been recorded so that the effects of varying the cell homogeneity and 
annexation thresholds may be monitored. The CPU times required to perform 
the perpoint classifications have been adjusted to reflect the increased 

efficiency of the LARSYS perpoint classifier which is coded in 
assembly language. Thus, the CPU time recorded for a perpoint classifi- 
cation is what a FORTRAN classifier would have required to perform the 
classification. 

The indices of classification performance were applied in several 
ways. Classification accuracy (identification) was evaluated utilizing 
field center pixel, "full field" and test field sample performance for 
all data sets. Proportion estimation was carried out for the LANDSAT 
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and Simulated Thematic Mapper data sets. 

The training performance is the overall classification accuracy 
(number of training pixels correctly classified divided by the total 
number of training pixels) of the pixels used to calculate the class 
statistics. Field center pixel performance is the overall classification 
accuracy of field center pixels. Field center pixels are pixels inset 
at least one pixel from the field boundaries. For the registered LACIE/SRS 
data the field center pixels are inset at least two pixels from the field 
boundary. Although this procedure insures that the pixels examined are 
not mixture pixels, it has the unfortunate effect of eliminating smaller 
fields from consideration. The third measure of classification accuracy, 
"full field" performance, includes those pixels on the boundaries of the 
fields in the classification performance. The "full field" pixels were 
generated by expanding the field center pixel boundaries . one pixel in all 
directions. 

The RMS error of informational class proportion estimates for each 
flightline was found by calculating the percent of the flightline classi- 
fied as a particular class and comparing it with the ground-collected 
estimate using equation (1). 



where, N = number of' informational classes 

C_^ = percent classified as informational class i 

C'^ = percent of class i estimated from ground-collected data 
RMS error is calculated for the LANDSAT and Thematic Mapper data 
runs. The Agricultural Stabilization and Conservation Service (ASCS) 
provided the ground truth proportion estimates for the simulated Thematic 
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Mapper data set. Proportion estimates for the 1974 LACIE/SRS segments 
were provided in ground truth packets received from JSC. The SRS county 
proportion estimates were used to calculate RMS proportion error for 
the CITARS data set. 

Average variability is a measure of the rate of change from one 
information class to another. It should reflect the degree to which ECHO 
reduces the "salt and pepper" effect which is sometimes present in per- 
point classifications. Variability is calculated by systematically 
selecting 50 lines of the classified area, counting the number of informa- 
tion class changes, and dividing by the number of opportunities for class 
changes. 

Variability = NCC/ (50*(NS-1)) (2) 

Where 

NCC = the number of class changes over the 50 selected lines 

NS = the number of classified pixels/line 
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4. Results 

This section outlines the results of the tests performed on the 
Supervised processor. Results are discussed separately for each scanner 
type, LANDSAT, simulated Thematic Mapper, and aircraft, followed by a 
summary of results for all three data types. 

a. LANDSAT Results 

Training and test information for the 10 data sets comprising the 
LANDSAT test data were drawn from the 1974 Kansas LACIE/SRS data sets 
and the 1973 CITARS data sets. One bitemporal and four unitemporal 
analyses were performed on the LACIE/SRS data sets. Six unitemporal 
analyses were performed on the CITARS data sets. 

i. LACIE/SRS Results 

Figures 5 through 10 present the average results for the four LACIE/SRS 
sites examined at a cell width of two for six cell homogeneity and four 
annexation parameters. Figure 5 plots the average CPU time in seconds 
required by the perpoint classifier (represented by the line of P's) versus 
the CPU time required by the average of the four LACIE/SRS data sets at each 
of 24 ECHO parameter settings for 2 by 2 pixel cells. The horizontal axis 
is the cell homogeneity threshold. As this threshold increases, it becomes 
more likely a cell will be classified as a unit, less likely that a cell 
will be split and its constituent pixels classified individually. The 
dependent variable, CPU time, appears on the vertical axis. The cell to 
field annexation parameter for each cell homogeneity threshold appears 
on the plot. When two or more annexation thresholds acheive the same 
performance, a star appears in that position on the plot. It can be seen 
from Figure 5 that as the cell homogeneity parameter increases, the CPU time 
required to classify the area decreases. Additionally, as the annexation 
parameter increases from zero to four, the CPU time required for the 
classification is reduced. These results are expected because fewer 
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FIGURE 9 
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TABLE 2.1-5 

EFFECT OF ECHO PARAMETERS ON INDIVIDUAL 1974 LACIE DATA SETS 
ABBREVIATIONS: 


CPU 

CENTRAL PROCESSING UNIT TIME 

FCP 

FIELD CENTER PIXEL PERFORMANCE 

FULL 

FULL FIELD PERFORMANCE 

TRAIN 

TRAINING FIELD PERFORMANCE 

RMS 

RMS PROPORTION ERROR 

VAR 

CLASSFICATION VARIABILITY 

CELW 

CELL WIDTH 

HOM 

CELL HOMOGENEITY THRESHOLD 

ANN 

ANNEXATION THRESHOLD 


ANALYSIS OF VARIANCE RESULTS (ENTRIES ARE SIGNIFICANCE LEVELS) : 
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HOM 

ANN 
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HOM 

ANN 
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Graham County, 

Kansas 






CPU 

.1 

.1 

7.5 

.1 

24.7 

99.9 

FCP 

.1 

.1 

18.7 

.1 

.1 

18.9 

FULL 

.1 

.1 

99.9 

.1 

1.6 

33.6 

TRAIN 

.1 

.1 

.1 

.1 

.1 

.2 

RMS 

.1 

.1 

.1 

.1 

.1 

99.9 

VAR 

.1 

.1 

.1 

.1 

.1 

.1 

Grant County, ! 

Kansas 






CPU 

.1 

.1 

.8 

.1 

2.2 

99.9 

FCP 

.1 

.1 

.1 

.1 

.1 

99.9 

FULL 

.1 

.1 

.1 

.1 

.1 

99.9 

TRAIN 

.1 

.1 

.1 

.1 

.1 

34.9 

RMS 

.1 

.1 

99.9 

.1 

99.9 

41.8 

VAR 

.1 

.1 

.1 

.1 

.1 

40.3 

Haskell County 

, Kansas 






CPU 

.1 

.1 

.2 

.1 

.1 

33.0 

FCP 

.1 

.1 

.1 

.1 

.1 

25.5 

FULL 

.1 

.1 

.3 

.1 

.1 

33.5 

TRAIN 

.1 

.1 

.1 

.1 

.1 

29.5 

RMS 

.1 

.1 

.1 

.1 

.1 

9.1 

VAR 

.1 

.1 

.1 

.1 

.1 

27.1 



2.1-33 


TABLE 2.1-5 (Continued) 

EFFECT OF ECHO PARAMETERS ON INDIVIDUAL 1974 LACIE DATA SETS 
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ANN 
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Kearny County, 

Kansas 






CPU 
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.1 

.1 

.1 

.4 

99.9 

FCP 

.1 

.1 

6.5 

.1 

.1 

99.9 

FULL 

.1 

.1 

1.6 

.1 

.1 

99.9 

TRAIN 

.1 

.1 

.3 

.1 

.1 

13.0 

RMS 

.1 

.1 

.1 

.1 

15.3 

2.0 

VAR 

.1 

.1 

.1 

.1 

.1 

.1 
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EFFECT OF ECHO PARAMETERS ON FOUR LACIE DATA SETS 



Cell 

Width 

Cell 
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Threshold 

Annexation 

Threshold 

Cell Width 
X 

Cell Homogeneity 
Threshold 

Cell Width 
X 

Annexation 

Threshold 
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2,6 

5,15 
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CPU Time 


1 

1 

.1 

99.9 
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Performance 

5 

5 


.1 

.4 

Full Field 

Performance 

5 

1 

25 

.1 

.6 

Training Performance 
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.1 

.1 

RMS Error 
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classifications are necessary when there are fewer cells which must be 
classified on a point-by-point basis (due to the increased cell homogeneity 
parameter) , and fewer classifications are required when a scene is composed 
of a few large objects rather than many small objects (the effect of 
increasing the annexation parameter) . Both of these effects are statistically 
significant at a 1% confidence level. 

The statistical significance levels of the various ECHO parameters on 
each of the six dependent variables (CPU time, field center pixel performance, 
"full field" performance, training field performance, RMS proportion error, 
and classification variability) for each of the LACIE data sets are 
presented in Table 2.1-5. Table 2.1-6 presents the statistical significance 
levels of the ECHO parameters on the dependent variables for the four LACIE 
data sets combined. 

Parameter -Considerations 

For the LACIE data runs, cell widths were sampled at 2, 

^average field size/2 + 1, andVaverage fields size; hence cell widths 
sampled were 2, 5 or 6, and 8 or 9. The field center pixel performances 
at the 2 by 2 cell size tended to be superior to those at 5 by 5 
or 6 by 6; RMS proportion estimation errors of 2 by 2 cells were about 
equivalent to those of the 5 by 5 or 6 by 6 cells; and less computer time 
was required for the 5 by 5 or 6 by 6 sized cells than the 2 by 2 cells. 

Graham County was also sampled at cell sizes 3 by 3 and 4 by 4 to 
observe what effects intermediate cell sizes might have. Results for the 
optimal values of these parameters are presented below: 
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Graham County, Kansas 



Optimal Results 



CPU 

FCP 

RMS 


Range 

Performance 

Error 

Cell Size 

(Sec) 

% 

% 

2x2 

109-73 

93.9 

0.1 

3x3 

101-59 

94.3 

0.1 

4x4 

95-50 

93.9 

0.0 

6x6 

98-51 

93.3 

0.7 

Per point 

121 

89.9 

1.2 


These results tend to indicate that cell width may be optimized for 
the LACIE/SRS runs when it is in the neighborhood of the fourth root of 
the average field size. This value is between three and four for Graham 
County. 

For the 2 by 2 tests (performed on all four LACIE/SRS sites) the 
optimal cell homogeneity parameter, based on the three fields performance 
criteria, appears to be 78. However, on the basis of RMS error, the optimal 
homogeneity parameter tested is 39. 

Due to the two pixel inset used to avoid mixture pixels and registration 
error for the LACIE/SRS field coordinates, small fields were not used for 
training or test purposes. Since only relatively large fields are included 
in the training decks, the effects that increasing the cell homogeneity 
parameter may have on the small fields is not reflected in the field center 
pixel performance results. Consequently, for the LACIE results, 39 is 
probably closer to the optimal homogeneity parameter than 78. Support for 
this conclusion is provided by the results in Figure 8. The "full field" 
performance measures the performance of the large test fields, including 
the points which contain the field boundaries. "Full field" test fields 


are composed of the pixels contained in the field center pixel test fields 
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plus those pixels on the borders of the field ("fringe pixels") . Since 
for the LACIE data sets, the performance of the field center pixels improves 
by .5% between cell homogeneity values 39 and 78 while the full field 
performance only improved by .2%, it would appear that the performance of 
the "fringe pixels" is degraded between homogeneity parameters 39 and 78. 
Those small fields not included in the test decks will be composed primarily 
of "fringe pixels". 

It can be noted by examining Figures 6 through 9 that the effect of 
annexation increases as the cell homogeneity parameter increases. This 
result is expected since a cell A only has the opportunity to be annexed 
to adjoining cell B when both cells A and B are homogeneous. The higher 
the cell homogeneity parameter, the more likely it becomes that both A and 
B will be homogeneous. 

Figure 10 is typical of all the variability results: as the homogeneity 

and annexation parameters increase, the classification variability decreases. 

LACIE/SRS Conclusions 

1 . CPU time and classification variablity decrease as cell homogeneity 
and cell to field annexation parameters increase (1% significance level). 

2. The optimal cell width parameter may fall between two and 
average field size/2 + 1. The fourth root of the average field size may be 
an appropriate value. 

3. For these data sets, the optimal annexation value appears to be 
in the neighborhood of two based on field center pixel and full field 
performances. 

4. For a cell width of two, the optimal cell homogeneity parameter 
is in the neighborhood of 39, based on RMS proportion estimate error. 
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ii. CITARS Results 

Figures 11 through 16 present plots for the 2 by 2 cell size Supervised 
ECHO results achieved from the CITARS data sets. Characteristic differences 
between the CITARS and the LACIE/SRS data sets include: 

* The CITARS data set has a much smaller average field size than the 
LACIE/SRS set. 

* The information classes are different. CITARS information classes 
are Corn, Soybeans and Other; LACIE/SRS classes are Wheat and Other. 

* The ground truth proportion estimates used for the LACIE/SRS sites 
were for the area of the LACIE segment. The ground truth proportion 
estimates used for the CITARS sites are for the whole county in 
which the data set lies, not for that area of the county which 

was actually sampled. 

Analysis of variance results for the six CITARS data sets are presented 
in Tables 2.1-7 and 2.1-8. Table 2.1-7 presents results for each data set 
individually. Table 2.1-8 presents the analysis of variance results for 
the combined CITARS data sets. 

As with the 1974 LACIE/SRS data sets, the CPU time decreased as cell 
homogeneity and annexation parameters increase. 

However, though Figure 5 (CPU time for the LACIE/SRS data set) and 
Figure 11 (CITARS CPU time) have the same general shape, there Is one 
important difference: the CPU time required to perform the ECHO classi- 

fication is a higher percentage of the CPU time required to perform the 
perpoint classification for each cell homogeneity value. This is a 
reflection of the smaller field size in the CITARS data set. Since fields 
are smaller, more boundaries and fewer homogeneous cells exist. As a result, 
more cells must be classified on a point-by-point basis and more CPU time 
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TABLE 2.1-7 

EFFECT OF ECHO PARAMETERS ON INDIVIDUAL CITARS DATA SETS 
ABBREVIATIONS : 

CPU CENTRAL PROCESSING UNIT (COMPUTER) TIME 

FCP FIELD CENTER PIXEL PERFORMANCE 

FULL FULL FIELD PERFORMANCE 

TRAIN TRAINING FIELD PERFORMANCE 

RMS RMS PROPORTION ESTIMATE ERROR 

VAR CLASSIFICATION VARIABILITY 

CELW CELL WIDTH 

HOM CELL HOMOGENEITY THRESHOLD 

ANN ANNEXATION THRESHOLD 

ANALYSIS OF VARIANCE RESULTS (ENTRIES ARE SIGNIFICANCE LEVELS): 
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.1 
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TABLE 2.1-7 (Continued) 

EFFECT OF ECHO PARAMETERS ON INDIVIDUAL CITARS DATA SETS 
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EFFECT OF ECHO PARAMETERS ON SIX CITARS DATA SETS 



Cell 

Width 

Cell 

Homogeneity 

Threshold 

Annexation 

Threshold 

Cell Width 
X 

Cell Homogeneity 
Threshold 

Cell Width 
X 

Annexation 

Threshold 

Cell Homogeneity 
Threshold 
X 

Annexation 

Threshold 

Degrees of Freedom 

2,10 

5,25 

3,15 




CPU Time 

1 

1 

1 

.1 

19.3 

99.9 

Field Center Pixel 
Performance 

10 

10 

25 

.1 

.1 

34.1 

Full Field 

Performance 

10 

5 

25 

.1 

.1 

99.9 

Training Performance 

1 

1 

25 

.1 

.1 

12.0 

RMS Error 

99 

10 

25 

.1 

.1 

28.2 

Classification 

Variability 

1 

1 

1 

.1 

.1 

.1 


Entries are Significance Levels 
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is required to perform those point-by-point classifications. For the 
smallest cell homogeneity parameter, ECHO takes more time on the average 
than the perpoint classifier. This is due to a great majority of the 
cells being split. For a cell that is split, the ECHO processor must 
first gather statistics and perform the homogeneity test, at some expense 
in CPU time. Then each point of the singular cell must be classified on 
a point-by-point basis which takes approximately the same CPU time as the 
perpoint classifier. Hence, when the great majority of cells are singular, 

ECHO is less efficient, in terms of CPU time, than the perpoint classifier. 

By examining the CITARS field center pixel performance (Figure 12) 
in light of the LACIE/SRS results, two conclusions can be drawn: 

1. The improvement in classification accuracy the ECHO procedure 
provides is reduced for the CITARS data sets, as compared to the LACIE set. 

2. The optimal annexation parameter is smaller for the CITARS data 
sets than for the LACIE data sets. 

It is evident by comparing Figures 12 and 14, that the CITARS training 
set is not representative of the test fields for the ECHO data sets. The 
performance of the training fields using ECHO is markedly superior to the 
performance of those same fields using the perpoint classifier. However, 
the field center pixel performance of the ECHO algorithm is only margi- 
nally superior in magnitude to the performance of the perpoint classifier. 

This indicates that on the data which the statistics truly represent, ECHO 
performance is superior; when the first and second order statistics fail 
to properly represent the classes, especially in areas with small field 
sizes, ECHO's performance may be no better than that of the perpoint classifier. 

The average RMS error achieved by ECHO for the six CITARS data sets 
is inferior to the average RMS error achieved by the perpoint classifier 
using the same class statistics. On a flightline by flightline basis, 



2.1-46 


ECHO's RMS proportion estimate error was higher than the perpoint RMS 
error for three of the data sets and lower than the perpoint RMS error 
for the other three data sets. This may be partially due to one or more 
of the following: 

* The RMS error values are generally high. The proportion estimates 
used for "truth" are the SRS estimates for the counties as a whole, 
not for the individual areas classified. Hence, if an area classi- 
fied is not a good approximation of the county as a whole, the 
validity of the RMS results must be questioned. 

* The training fields do not appear to be a good representation of the 
test fields. If the training statistics do not represent the scene 
as a whole, ECHO performance may suffer to a greater extent than 
the perpoint performance. 

* When the average field size is small, there will be a high proportion 
of "mixture" cells (cells containing points of more than one class) . 
Should a cell containing pixels of differing classes pass the homo- 
geneity criteria, it is likely to be categorized into the class 
with the highest variance. Thus as the cell homogeneity criteria 
becomes larger, the RMS error would tend to increase. (An observed 
result - see Figure 15.) 

The effect of cell width on RMS error was not significant when all 
the CITARS data sets are considered at once. However, cell width produced 


significantly different results in terms of RMS error for each of the 
six data sets taken individually. The effects of increasing cell width 
were opposite for different data sets. 
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Conclusions 

* As the field size becomes smaller, the improvement in classi- 
fication accuracy provided by the supervised ECHO algorithm is reduced. 

At the same time, due to the increased number of singular cells, the CPU 
time which Supervised ECHO requires increases. 

* When the training set is representative of the data, ECHO's 
classification performance continues to be superior to that of the per- 
point classifier for scenes with average fields at or below 16 to 20 pixels. 

* Due to the mixture-cell, large-variance-class problem, ECHO may 
or may not degrade RMS proportion Estimation for areas with relatively 
small fields. 
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ill. Bitemporal Results 

Figures 17 through 22 illustrate the Supervised ECHO results for the 
Bitemporal analysis. Figures 17a through 22a are the results of a uni- 
temporal analysis of the same area. The training for the bitemporal analysis 
was done independently of the training for the unitemporal analysis. Com- 
paring Figures 18 and 18a it can be seen that the ECHO field center pixel 
performance is a greater improvement over the perpoint classifer’s perfor- 
mance in the unitemporal case. This may partially be a result of some 
misregistration between dates making object extraction more difficult and 
partially a function of the different training statistics used for the 
separate classifications. 

With the larger number of channels in the multidate classification, 
the optimal homogeneity threshold is higher than it would be for the single 
date classification. This is expected as, for the Supervised algorithm, 
the homogeneity criteria approximates to a chi-squared statistic with 
(number of channels) times (cell width squared) degrees of freedom. For 
the same significance level of test, the cell homogeneity criteria should 
be higher for the eight channel data than for the four channel data. 

Based on the field center pixel performance and the RMS proportion 
error results, the optimal settings of the homogeneity and annexation 
parameters appear to be in the neighborhood of 180 and 1, respectively. 

Note that although the multitemporal analysis does not greatly improve 
the field center pixel classification accuracy, it does (especially for 
the perpoint classifier) substantially reduce the RMS proportion estimate 
error. In the Bitemporal analysis the perpoint RMS error is only .2%, 
leaving ECHO little room for improvement. Note also that for several ECHO 
parameter settings in the unitemporal analysis, ECHO performs with an RMS 
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proportion estimate error of under .5% while the perpoint classifier is 
functioning with a 6.5% RMS error. 

In the multitemporal case, ECHO requires less CPU time than the 
perpoint classifier for all parameter settings. This is not true in the 
unitemporal case. 


Bitemporal Conclusions 

* The Supervised ECHO processor produced classifications superior 
to those produced by the perpoint classifier in terms of CPU time and 
field center pixel accuracy for the frame analyzed. 

* The optimal homogeneity parameter for bitemporal analysis appears 
to he much larger than for the unitemporal case. 

* The CPU time required and classification variability produced by 
ECHO are smaller than those of the perpoint classifier. 
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iv. Overall LAND SAT Results 

Table 2.1-9 summarizes the effects of the ECHO parameter settings 
on the six dependent variables. Figures 23 through 28 plot the average 

t 

ECHO performance at the various ECHO parameter settings versus the average 
perpoint performance for the ten LANDSAT data sets. 

The cell width setting effects are significant for all dependent 
variables for all LANDSAT data sets considered individually, and are 
significant for all dependent variables except SMS proportion error when 
all LANDSAT data sets are considered together, using data sets as blocks. 

When the LACIE and CITARS LANDSAT runs are considered separately, 
cell width is significant for CPU time, training performance, and classi- 
fication variability for the CITARS data sets and for CPU time, field 
center pixel performance, "full field" performance, and classification 
variability in the LACIE data sets. 

The statistically significant (5% significance level) overall cell 
width effects for the LANDSAT data sets are: 

As cell width increases, CPU time decreases. 

As cell width increases, field center pixel, "full field" and 
training field performances decrease. 

As cell width increases, classification variability decreases. 

Though not statistically significant, the trend is for RMS error to 
increase as cell width increases. Two things should be noted about the 
cell width overall results: 

* It is not uniformly true that all cell width increases, field 

center pixel, "full field" and training field performances decrease. 
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EFFECT OF ECHO PARAMETERS ON SIX VARIABLES USING DATA SET 




AS A BLOCK 

(ALL TEN LANDS AT DATA SETS) 
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* Since not all cell widths were sampled between two and the 

\J average field size/2 + 1, it may be true that the optimal cell 

width for a data set lies between these two parameters . Indications 

from the Graham County analysis are that the optimal cell width 

may be in the neighborhood of the fourth root of the average field size. 

The cell homogeneity threshold is significant for all dependent variables 
both when all LAND SAT data sets are considered individually and when all 
LANDSAT data sets are considered together. 

The cell homogeneity threshold is significant for CPC time, "full 
field" performance, training performance t and classification variability 
for the combined CITARS data sets and is significant for CPU time, field 
center pixel performance, "full field" performance, training performance 
and classification variability in the combined LACIE data sets. 

The statistically significant overall effects of the cell homogeneity 
parameter are: 

As the homogeneity parameter increases, CPU time and variability 
decrease . 

As the homogeneity parameter increases, field center pixel, "full 
field" and training performances decrease. Figures 24 through 26 
indicate, however, that for the first two or three homogeneity 
threshold values, as the homogeneity threshold increases, the 
performances also increase. 

As the homogeneity parameter increases, RMS error shows a trend towards 
increasing (Figure 27) though this trend is not significant at a 5% confidence 
level . 

The annexation parameter is significant in 86% of the 60 cases in the 
LANDSAT data sets. It is significant only for CPU time and classification 
variability when the CITARS or the LACIE or all the LANDSAT data sets are 
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considered together. This indicates that the proper annexation 
parameter is highly data dependent. As the annexation threshold increases, 
CPU time and classification variability tend to decrease. 

The interaction of cell width and the cell homogeneity threshold is 
significant in 95% of the 60 cases of 6 dependent variables and 10 
LANDSAT data sets. This interaction is significant for all dependent 
variables when the CITARS or the LACIE or all LANDSAT data sets are 
considered together. 

The interaction of cell width and annexation parameter is significant 
in 93% of the 60 LANDSAT cases. When all the LANDSAT data sets are 
considered together, the cell width-annexation interaction is significant 
for all dependent variables except CPU time. This is also true for the 
CITARS data sets. However, for the LACIE data sets, the interaction is not 
significant for RMS proportion error. 

The cell homogeneity-annexation interaction was significant in 30% 
of the 60 individual LANDSAT data set ANOVAs, being significant in 90% 
of the data sets for variability and 50% of the data sets for RMS proportion 
error. When the 10 LANDSAT data sets are considered together, the cell 
homogeneity-annexation interaction is significant for training performance 
and variability. When considering only the CITARS or only the LACIE data 
sets, the interaction is significant for variability alone. 
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ECHO-Perpoint LANDS AT Comparison 

When the optimal parameter setting for each dependent variable was 
selected from the ECHO results and compared to the perpoint results for 
the 10 LANDSAT data sets the following results were obtained: 

ECHO was faster than the perpoint classifier (.1% confidence level). 
ECHO’s field center pixel, "full field”, and training field 
performances were superior to those of the perpoint classifier 
(1.3, 5.7, .1% confidence levels, respectively). 

ECHO had a lower RMS proportion estimate error (4.4% confidence 
level) . 

ECHO had less classification variability than the perpoint 
classifier (.1% confidence level). 
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b. Simulated Thematic Mapper Results 

The analyses of variance for the effects of the ECHO parameters on 
the eight simulated Thematic Mapper data sets (two sites, four resolutions) 
are presented in Table 2.1-10. Results for the eight data sets considered 
together are presented in Table 2.1-11. 

Cell width is significant for all dependent variables and all Thematic 
Mapper data sets individually. When the Thematic Mapper data sets are 
considered together, using site as a block and resolution as a fixed factor, 
cell width is significant only for training performance and classification 
variability. As cell width increases, classification variability decreases. 

As cell width increases, training performance decreases. Trends which 

4 

can be noted, but which are not statistically significant at the sample size 
used include: as cell width increases, CPU time, field center pixel perfor- 

mance, and "full field" performance decrease. It may be true that a cell 
width between two and the square root of the average field size divided by 
two would be superior to those cell sizes tested in terms of CPU time, RMS 
error, and field center pixel, "full field" and training field performances. 

The cell homogeneity parameter is significant for all dependent vari- 
ables in all Thematic Mapper data sets individually. When the data sets 
are considered together, the cell homogeneity parameter is significant for 
CPU time, "full field" performance, training performance and variability. 

As the cell homogeneity parameter increases, the CPU time and classification 
variability decreases (1% significance level) . "Full field" and training 
field performance decrease over the range of homogeneity threshold tested (1% 
confidence level) . The graph indicates that RMS proportion error increases 
as the cell homogeneity criteria increases. However, this is not a 
consistent trend in the eight data sets and is not statistically significant. 
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TABLE 2.1-10 

EFFECT OF ECHO PARAMETERS ON INDIVIDUAL SIMULATED THEMATIC MAPPER DATA SETS 
ABBREVIATIONS: 

CPU CENTRAL PROCESSING UNIT (COMPUTER) TIME 

FCP FIELD CENTER PIXEL PERFORMANCE 

FULL FULL FIELD PERFORMANCE 

TRAIN TRAINING FIELD PERFORMANCE 

RMS RMS PROPORTION ERROR 

VAR CLASSIFICATION VARIABILITY 

CELW CELL WIDTH 

HOM CELL HOMOGENEITY THRESHOLD 

ANN ANNEXATION THRESHOLD 






CELW 

CELW 

HOM 


CELW 

HOM 

ANN 

X 

X 

X 





HOM 

ANN 

ANN 

Williams County, 

North 

Dakota - 

30 meter 

resolution 



CPU 

.1 

.1 

99.9 

.1 

99.9 

31.8 

FCP 

.1 

.1 

.1 

.1 

.1 

20.5 

FULL 

.1 

.1 

.1 

.1 

.1 

28.7 

TRAIN 

.1 

.1 

4.9 

.1 

2.3 

99.9 

RMS 

.1 

1.0 

4.4 

.1 

4.1 

44.7 

VAR 

.1 

.1 

.1 

.1 

.1 

46.3 

Williams County, 

North 

Dakota - 

40 meter 

resolution 



CPU 

.1 

.1 

99.9 

1.6 

99.9 

99.9 

FCP 

.1 

.1 

.1 

.1 

.1 

99.9 

FULL 

1.9 

.1 

8.9 

.1 

3.0 

40.5 

TRAIN 

.1 

.1 

99.9 

.1 

99.9 

99.9 

RMS 

.1 

.1 

99.9 

.1 

99.9 

27.1 

VAR 

.1 

.1 

99.9 

.5 

99.9 

46.0 

Williams County, 

North 

Dakota - 

50 meter 

resolution 



CPU 

.1 

.1 

99.9 

.1 

99.9 

99.9 

FCP 

.1 

.1 

99.9 

.1 

9.1 

99.9 

FULL 

.1 

.1 

99.9 

.1 

36.3 

99.9 

TRAIN 

.1 

.1 

5.6 

.1 

42.7 

34.2 

RMS 

.1 

.1 

.1 

.1 

.8 

.1 

VAR 

.1 

.1 

.1 

.1 

.1 

15.5 
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TABLE 2.1-10 (Continued) 

EFFECT OF ECHO PARAMETERS ON INDIVIDUAL SIMULATED THEMATIC MAPPER DATA SETS 






CELW 

CELW 

HOM 


CELW 

HOM 

ANN 

X 

X 

X 





HOM 

ANN 

ANN 

Williams County, North 

Dakota - 60 meter resolution 



CPU 

.1 

.1 

99.9 

.1 

99.9 

99.9 

FCP 

.1 

.1 

19.8 

.1 

8.7 

99.9 

FULL 

.1 

.1 

2.1 

.1 

.7 

99.9 

TRAIN 

.1 

.1 

14.6 

.1 

99.9 

99.9 

RMS 

.3 

.1 

99.9 

9.2 

26.4 

99.9 

VAR 

.1 

.1 

.1 

.1 

.1 

11.6 

Finney County, 

Kansas - 

• 30 meter 

resolution 




CPU 

.1 

.1 

39.7 

.1 

15.0 

3.4 

FCP 

.1 

.1 

.1 

.1 

.1 

29.4 

FULL 

.1 

.1 

.1 

.1 

.1 

38.6 

TRAIN 

.1 

.1 

.1 

.1 

.1 

99.9 

RMS 

.1 

.1 

16.7 

.1 

44.4 

24.9 

VAR 

.1 

.1 

.1 

.1 

.1 

37.5 

Finney County, 

Kansas - 

- 40 meter 

resolution 




CPU 

.1 

.1 

10.4 

.1 

1.5 

9.4 

FCP 

.1 

.1 

.1 

.1 

.1 

46.2 

FULL 

.1 

.1 

.1 

.1 

.1 

47.4 

TRAIN 

.1 

.1 

.1 

.1 

.1 

99.9 

RMS 

.1 

.1 

7.8 

.1 

4.4 

99.9 

VAR 

.1 

.1 

.1 

.1 

.1 

43.4 

Finney County, 

Kansas - 

- 50 meter 

resolution 




CPU 

.1 

.1 

99.9 

.1 

99.9 

99.9 

FCP 

.1 

.1 

9.6 

.1 

2.7 

99.9 

FULL 

.1 

.1 

.3 

.1 

17.3 

29.7 

TRAIN 

.1 

.1 

7.5 

.1 

8.7 

99.9 

RMS 

.1 

.1 

.1 

.1 

.1 

2.3 

VAR 

.1 

.1 

.1 

.1 

.1 

44.0 

Finney County, 

Kansas • 

- 60 meter 

resolution 




CPU 

.1 

.1 

99.9 

.1 

99.9 

99.9 

FCP 

.1 

.1 

.1 

.1 

.1 

47.0 

FULL 

.1 

.1 

.7 

.1 

.1 

99.9 

TRAIN 

.1 

.1 

.7 

.1 

.2 

99.9 

RMS 

.1 

.1 

16.9 

.1 

13.4 

47.3 

VAR 

.1 

.1 

.1 

.1 

.1 

99.9 
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EFFECT OF ECHO ON SIX VARIABLES USING SITE AS A BLOCK AND RESOLUTION AS A FIXED FACTOR 


(Eight Thematic Mapper Data Sets) 


Degrees of Freedom 
CPU Time 

Field Center Pixel 
Performance 

Full Field 

Performance 

Training Performance 

RMS Error 

Classification 

Variability 


Cell 

Width 

Cell 

Homogeneity 

Threshold 

Annexation 

Threshold 

Cell Width 
X 

Cell Homogeneity 
Threshold 

Cell Width 
X 

Annexation 

Threshold 

Cell Homogeneity 
Threshold 
X 

Annexation 

Threshold 

2,2 

5,5 

3,3 




10 

1 


.1 

99 

99 

25 


25 

.1 

99 

99 

25 

1 

25 

.1 

99 

99 

5 

1 


.1 

99 

99 


25 

10 

.1 

99 

99 

1 

1 

1 

.1 

.1 

99 
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Figure 29 

Simulated Thematic Mapper CPU Time 
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Figure 30 

Simulated Thematic Mapper Field Center Pixel Performance 
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figure 31 

Simulated Thematic Mapper Full Field Performances 
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Simulated Thematic tapper Training Performances 
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Figure 33 

Simulated Thematic Mapper RMS Proportion Errors 
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The annexation parameter is significant in 63% of the 48 cases involving 
6 dependent variables and 2 sites and 4 resolutions of thematic mapper data. 
It is never significant for CPU time in the individual data sets. When the 
Thematic Mapper data sets are considered together, the annexation parameter 
is significant only for variability. As the annexation parameter increases, 
classification variability decreases. 

The interaction of cell width and cell homogeneity is significant for 
all dependent variables both when the Thematic Mapper data sets are 
considered together or individually. 

The interaction of cell width and annexation is also significant in 
63% of the 48 cases in the Thematic Mapper data sets. It is significant 
only for variability when all the Thematic Mapper data sets are considered 
together. 

The interaction of cell homogeneity and annexation is significant in 
8% of the 48 cases of the individual Thematic Mapper data sets and is not 
significant for any dependent variable when the Thematic Mapper data sets 

K 

are considered together. 

Comparison of the ECHO and the Perpoint Classifier 

F-values for the analysis of variance comparing the effects of 
classification method (ECHO versus perpoint) and resolution size (30, 40, 50 
and 60 meters) to the six variables for the optimum parameter setting of 
the Supervised ECHO classifier for each of the dependent variables is 
presented in Table 2.1-12. 
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Table 2.1-12 



Method 

Resolution 

CPU Time 

1204.39* 

9.86** 

Field Center Pixel Performance 

147.63*** 

1.25 

Full Field Performance 

74.21*** 

1.47 

Training Field Performance 

15.21+ 

.62 

RMS Proportion Error 

2.62 

.21 

Classification Variability 

13547.44** 

7.23*** 

Degrees of Freedom 

(1,1) 

(3,3) 


* significant at 1% confidence level 

** significant at 5% confidence level 

*** significant at a 10% confidence level 

+ significant at a 25% confidence level 

Classification variability is reduced as resolution is increased. This 
result is expected since the field which is 10 pixels wide in the 30 meter 
data is only 5 pixels wide in the 60 meter data. The CPU time required to 
perform the classification is increased as the resolution increases. 

Results based on classification method are hampered by having only 
(1,1) degrees of freedom. ECHO requires less CPU time and has less classi- 
fication variability than the perpoint classifier at a 1% confidence level. 

ECHO field center pixel and "full field" performance appear to be superior 
to the perpoint classifier at a 10% confidence level. 

More repetitions are necessary before it can be deomonstrated that ECHO 
RMS error is significantly less than that of the perpoint classifier for the 


Thematic Mapper data. 



2.1-73 


c. Aircraft Results 

Tables 2.1-13 and 2.1-14 summarize the analysis of variance results 
for the effect of the ECHO parameters on the three aircraft data sets 
taken individually and taken as a block. Figures 35 through 39 present the 
CPU time, field center pixel performance, "full field" performance, 
training field performance, and classification variability. RMS values 
were not calculated because ground truth proportion estimates were not 
available for all three aircraft sites. The aircraft data sets were sampled 
only at the 2 by 2 cell size. 

As the cell homogeneity parameter increases, field center pixel and 
training field performance increase (5% confidence levels) . As with the 
LANDSAT and simulated Thematic Mapper data sets, classification variability 
decreases as annexation and cell selection threshold increase. 

The effect of the annexation parameter increases as the cell homogeneity 
parameter increases for training field performance, classification variability 
and CPU time. 

ECHO - Perpoint Comparison 

Though only three aircraft data sets were tested, the Supervised ECHO 
processor produced results which had significantly less variability than 
the results produced by the perpoint classifier (5% confidence level) . The 
comparison between the perpoint and the ECHO processor is approaching 
statistical significance for the other four variables. The following 
statements can be made f° r the aircraft data sets at a 15% confidence level. 

* ECHO requires less computer time than the perpoint classifier, 

* ECHO has higher field center pixel, "full field" and training 
performances than the perpoint classifier. 
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For the aircraft data sets, the cell homogeneity value appears 
to maximize performance for a value in the neighborhood of 70 and an 
annexation threshold of 4 seems to be optimal. 
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Table 2.1-13 


INDIVIDUAL AIRCRAFT DATA SETS 
(Entries are Significance Levels) 


Finney County. North Dakota* 


Annexation Threshold 


Cell Homogeneity 
Threshold 


CPU Time 

Field Center Pixel 
Performance 
Full Field Performance 
Training Field Performance 
Classification Variabili ty 


Williams County, Kansas* 


Annexation Threshold 


Cell Homogeneity 
Threshold 


CPU Time 

Field Center Pixel 
Performance 

Full Field Performance 
Training Field Performance 
Classification Variability 


aecanoe County, Indiana* 


Annexation Threshold 


Cell Homogeneity 
Threshold 


CPU Time 

Field Center Pixel 
Performance 

Full Field Performance 
Training Field Performance 
Classification Variability 


* Degrees of freedom for the Annexation Threshold test are 3,15 and Degrees of 
freedom for the Cell Homogeneity Threshold are 5,15. 
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Table 2.1-14 

RESULTS FOR ECHO CLASSIFICATIONS OF 
AIRCRAFT DATA, USING ALL THREE DATA SETS 



Cell 

Homogeneity 

Annexation 

Threshold 

Interaction 

Degrees of Freedom 

5,10 

3,6 

15,30 

Field Center Pixel 
Performance 

3.92*** 

2.516° 

1.609° 

Full Field Performance 

0.92 

2.88 

1.75*** 

Training Field Performance 

6.84* 

4.64*** 

5.47* 

Classification Variability 

18.80* 

63.33* 

21.78* 

CPU Time 

1.81° 

1.78° 

3.20* 

* sig 1% 

** sig 5% 

Table entries 

are F values. 


*** sig 10% v 

— Sig. level indicated by 
symbols . 



° slg 25% 
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Figure 35 
Aircraft CPU time 
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Figure 36 

Aircraft Field Center Pixel Performance 
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Figure 37 

Aircraft Full Field Performance 
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Figure 38 

Aircraft Training Field Performance 
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Figure 39 

Aircraft Classification Variability 
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C. Products for LACIE/AI Evaluation 

Magnetic computer tapes containing examples of object maps produced 
by the Nonsupervised field extraction algorithm are provided in Universal 
format. Object map tapes are multispectral image storage tapes which 
have had the data values altered for pixels falling within objects 
identified by the algorithm. The mean response for the pixels comprising 
the object replaces the individual response for each constituent pixel 
of the object. 


SUMMARY 

The ECHO processor developed at LARS has been implemented in documented 
FORTRAN programs. 

The Supervised ECHO processor has been tested over LANDSAT, simulated 
Thematic Mapper and aircraft data sets. Evaluations of those tests 
indicate that: 

* ECHO is less expensive than the perpoint classification algorithm 
(statistically significant for LANDSAT and simulated Thematic 
Mapper data sets at a 1% confidence level and at a 15% confidence 
level for the aircraft data sets) , 

* ECHO provides better field center pixel and "full field" classi- 
fication performance than the perpoint classifier (5% confidence 
level on LANDSAT data, 10% confidence level on simulated Thematic 
Mapper data, 15% confidence level on aircraft data) . 

* ECHO RMS proportion error is significantly less than that of the 
perpoint classifier at optimal parameter settings for the LANDSAT 
data. Proportion estimation for the Thematic Mapper data set was not 
significantly different for the two classification^ methods. 
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* ECHO classification results are less variable than the classi- 
fication results of the perpoint classifier. This result is 
statistically significant at a 5% confidence level for all data 
sets. 

The effects of the three ECHO parameters on six variables were monitored. 
In general, cell homogeneity had the strongest effect on the six variables. 
Cell width had a strong effect on the variables when the data sets were 
considered individually. However, when the data sets were considered 
together, the effect of cell width was not significant for RMS error and 
was less significant than the cell homogeneity parameter for field center 
pixel and "full field" performances, indicating the effects of cell width 
are sometimes opposite in direction for differing data sets. The annexation 
threshold had a significant effect on classification variability and CPU 


time but a weak effect on the other variables. 
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CONCLUSIONS 

ECHO successfully exploits the redundancy of states characteristics 
of sampled imagery of ground scenes to achieve better classification 
accuracy , reduce the number of classif ications required, and reduce the 
variability of the classification results. The information required 
to produce ECHO classifications are cell size, cell homogeneity, and 
cell-to-field annexation parameters, input data, and a class-conditional 
marginal density statistics deck. 

Future research should be directed towards developing methods for 
utilizing the information produced by the ECHO processors to aid in the 
training process. For example, the singular cell map produced by the 
Supervised ECHO processor (see Figure 40) provides information on training 
adequacy. Singular cells are represented on this map by 0*s. For the 
Supervised processor, cells are identified as being singular when the 
likelihood of the cell belonging to the most likely of the available 
classes falls below a threshold. Therefore, when a cell is categorized 
as "singular”, it is either because it contains pixels from more than one 
class, or because the spectral class of the field is not represented in 
the available class statistics. Large groups of contiguous singular cells 
will occur when one or more spectral classes have been omitted. The 
singular cell map may indicate where additional training statistics should be 
collected. 

The intermediate tape produced by the Nonsupervised ECHO processor 
contains the statistics for each homogeneous object identified by the field 
extraction routine. It is possible that this spectral-spatial cluster will 
contain all the information an analyst will need to produce a class statistics 
deck superior to what would be produced by traditional training procedures 
for a smaller investment of analyst time than is now required. 
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Example of the Supervised ECHO Singular Cell Map 
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Lastly, background research should be performed to determine what 
training techniques are optimal for a field classifier. ECHO performance 
may be enhanced by utilizing training techniques developed specifically 
for it, rather than techniques which are optimal for the perpoint 
classifier. 
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1. Module Usage 
SECSUP 

Calling sequences 
CALL SECSUP 

SECSUP is the supervisor routine for the supervised ECHO 
classifier. It is called with no arguments. 

2 . Internal Description 

The program SECSUP basically controls the program se- 
quencing for the many subprograms involved in the supervised 
ECHO classifier. Two subroutines are called by SECSUP: SECRDR 

and SECINT. 

SECRDR is the control card reader for reading the user input 
data cards and is called first. SECINT is the initialization 
routine for the SECHO function and is called to handle utili- 
zation of all input/output operations as well as classification 
processing. For more detailed descriptions of these modules, 
see their program abstracts. 

Two common blocks are used in SECSUP. GLOCOM, the main 
LARSYS common block and SECCOM, the common block for the super- 
vised ECHO function are included. 

3 . Input Description 

Not applicable 

4 . Output Description 

Two messages are produced and written to unit TYPEWR, the 
console. 

SUPERVISED ECHO FUNCTION REQUESTED 
Signifies beginning of the function. 
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SUPERVISED ECHO FUNCTION COMPLETED 
Signifies end of the function. 

5. supplemental Information 

Two COMMON blocks are included in SECSUP: GLOCOM 

SECCOM 
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FORTRAN IV G LEVEL 
FILE SECSUP 


20.7 


SECSUP 


DATE s* 77 129 


15016012 


PAGE 0001 


0001 

0002 

0003 


0004 

0005 

0006 

0007 

0008 

0009 

0010 


0011 


0012 

0013 

0014 


0015 

0016 
0017 


0010 

0019 

0020 
0021 
0022 


C SECSUP LARS XXXX SEC00010 

C SEC00020 

q***************************************************************«*******SEC00030 

C SEC00040 

C SECSUP - SUPERVISOR ROUTINE FOR SUPERVISED ECHO FUNCTION SEC00050 


C 

c 


c 

c 


c 

c 

c 

c 


WRITTEN 01/20/77 BY P.D* ALENDUFF 


SEC00060 

SEC00070 


SUBROUTINE SECSUP 
IMPLICIT INTEGER * 4 (A-Z) 


CHKOUT. COPFIL. CLASSRt CLASSX. 


COMMON /GLOCOM/ BLANK, CARD (20 1, umkuui, tuKMii cl 
CLUSTX, CONPUT, CPYOUT, CRDRDR, CRDSEQ, OAT APE, 
DUPLTP ♦ DUPRUN, ERRMSG, FBPNT * 

FILESV, FLDBND, HDATA, HEADtBo), ID1200), 1 IMAGEX, 
IMARK, KEYBD » HAPTAP, MAXCHA, MAXCLS, 

PAGSIZ, PNCH. POINT, PRESUX. PRNTR, READIN, 
RESTRT, RUNFtL, RUNTA8( 


i to 3 I 

SD ATAi ; SE PARX^; SE PTPX * S P XrI f 10 > ♦ TEMP AS (30), 


iuflifl* jtrMMA* dcrirA. 

TPSTAT (6) , TTFLDX ,TYPEWft, 
TOP, ARRAY ('12500) 


£*********^***♦*♦♦**♦*♦*♦*♦♦**♦*♦***♦**♦**♦♦♦♦♦**♦♦♦*****♦♦*♦♦*♦*♦♦*♦♦*$£000080 

SEC00090 
S6C00100 
SECOOHO 
SEC 00120 
SEC 00130 
SEC00140 
SEC00150 
SEC00160 
SEC00170 
SEC00180 
SEC00190 
SEC00200 
SEC0021 0 
SEC00220 
SEC 00230 
SEC00240 
SEC00250 
SEC00260 
SEC00270 

BLANK2 SEC00280 

CHKOUT SEC 00290 

BLANK1 SEC00300 

(DATSAV* 1D{ 1 ) ) , ( CURRUN, I D( 3 ) ) , ( FRQCAL (1) . ID( 5 1) ) , SEC00310 

(HED1 f 1 ) , HEAD (8) ) * (DATE! U,HEAD(26) ), ( HED2( 1 ) ,HEAD( 39) ) , SEC 00320 
( TIME ( 1 ) .HEAD (58 ) I , ( COMENT ( 1 1 , HE A0< 72 I ) , SEC00330 

(MAPSAV, TPSTAT ( 1 ) > , SEC00340 

( SEP SCR , TPSTAT (2 ) ) , ( DUPI N.TPSTAT ( 3 > ) * ( CAST AT , TPSTAT (4)), SEC 00350 
( C OP SER, TPSTAT (5 ) ) , ( TRAQUT,TPSTAT(6 )) , SEC00360 

t BLANK , BLANK2 , BL ANK1 ) SEC00370 

SEC00380 
SEC00390 

/SECCOM/ BUFROZ,CSELEC 1 CSET(3 t 30) ,CSET3 ( 3, 30 ),THRES ,CELSIZ, SEC00400 


REAL 

REAL 


* 8 ARRAY 

* 4 

INTEGER * 
INTEGER * 2 
LOGICAL * 4 
LOGICAL * l 
EQUIVALENCE 


I 

FRQCAL (5.30) 
4 COMENT ( 


16), DATE( 5 ), HEDK16), HED2U6I, TIME(5) 


2 

3 

4 

5 
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COMMON 


aC LtL> f COt I Uf ttoCl J \ Df DVJ i \ _ . . 

INTTAP, INTFIL, JPTS, LINES, NOCLS, NOFET , SEC00410 
NVR, OUTPUT £ POLN AMI 2^60 ) PRtFIXi 2 ) , SEC0042C 
■“ “ k SEC00430 


1 CEL«tH,INF0ll7) , INPUT i 

2 N0FET3 , NO POOL , NWORD, NW0RO2 * m»n, uuirui, ruLH«nu*oui# rt\ 

3 PTS, RQTAPE, RQFILE, SYHCNT, VARSZ3. CSELt3Q), CSEL3(3Q), 

3 FETVEC < 30 ) , FETVC3(30>, POLPTR ( 2 ,60 ] , P0LSTK(60>, 

4 CDFLAG, CLSMAP, CSET1 (3,30) , OBJMAP, PHASE1, PHASE2, POLNMK 60 ) , 

5 PRSTAT, SYM (60 > , SYMMTX160) 


REAL * 4 CSELEC, CSET, CSET3, THRES 

INTEGER * 2 CSEL, CSEL3, FETVEC, FETVC3, POLPTR, POLSTK 


LOGICAL 
l PRSTAT 


♦ 1 CDFLAG, CLSMAP, CSET1, OBJMAP, PHASE1, PHASE2, 


100 


C 

C 

C 


CALL READER AFTER INITIAL MESSAGE PRODUCED 
WRITE (TYPEWR.IOO) 

FORMAT ( 1 10000 SUPERVISED ECHO FUNCTION REQUESTED 
CALL SECRDR 


(SECSUP)* ) 


CALL PROCESSING ROUTINE THEN PRINT TERMINATION MSG 


CALL SECINT 

WRITE ( TYPE HR, 9000) 

9000 FORMAT (■ 10000 SUPERVISED 
~ RETURN * * 

END 


ECHO FUNCTION COMPLETED ( SECSUP)?) 


SEC00440 
SEC00450 
SEC0046C 
SEC00470 
SEC00480 
SEC00490 
SEC0050C 
SEC00510 
P0LNM1 , SEC00520 
‘ SEC00530 
SEC00540 
SEC00550 
SEC00560 
SEC00570 
SEC00580 
SEC00590 
SEC00600 
SEC 00610 
SEC00620 
SEC00630 
SEC00640 
SEC00650 
SEC00660 
SEC.00670 
SEC00680 
SEC00690 
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SECCOM 

Variable Description 
BUFROZ 

CSELEC 
CSET (3,30) 

CSET3 (3,30) 

THRES 
CELS I Z 

CELWTII 
INFO (16) 


INPUT 

INTTAP 

INTFIL 

JPTS 

LINES 

NOCLS 


INTEGER* 4 Number of rows of cells in 
buffer. 

REAL* 4 Cell selection parameter. 

REAL* 4 Channel calibration code vector 
from statistics. 

REAL* 4 Channel calibration code vector 
from channels card.' 

REAL* 4 Annexation threshold . 

INTEGER* 4 Number of data points in a 
cell. I 

INTEGER* 4 Number of rows in cell. 
INTEGER* 4 Array for storing field de- 
scription information and areas to 
be classified. 

INTEGER* 4 Unit number for input tape. 

i 

INTEGER* 4 Tape number of intermediate 
tape. 

INTEGER* 4 File number of intermediate 
tape . 

INTEGER* 4 Number of cells horizontally 
in data area. 

INTEGER* 4 Number of rows of cells in 
area to be processed. 

INTEGER* 4 Number of classes in statis- 


tics deck 
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NOFET 

*s 

N0FET3 

NOPOOL 

NWORD 

NWORD2 

NVR 

OUTPUT 

POLNAM (2,60) 
PREFIX (2) 

RQTAPE 

RQFILE 

SYMCNT 

VARS 2 3 


CSEL (30) 


INTEGER* 4 Number of channels in sta- 
tistics deck. 

INTEGER* 4 Number of channels to be 
used in classification. 

INTEGER* 4 Number of pools to be used 
in classification. 

INTEGER* 4 Integer array word counter. 
INTEGER* 4 Integer word counter. 
INTEGER* 4 Number of channels on data 
tape. 

INTEGER* 4 Unit number, of output file. 
INTEGER* 4 Names of statistics classes. 
INTEGER* 4 Prefix information used in 
writing results. 

INTEGER* 4 Tape number of the results 
tape . 

INTEGER* 4 File number of the results 
file. 

INTEGER* 4 Number of user supplied 
symbols . 

INTEGER* 4 Number of words to store 
half of a triangular covariance matrix 
for classes (channels * (channels + 
D)/2. 

INTEGER* 2 Channel select vector from 


statistics 
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GSEL3 (30) 


FETVEC (30) 


FETVC3 (30) 


POLPTR (2,60) 


POLSTK (60) 


CDFLAG 


CLSMAP 


CSET1 (3, 30) 


OBJI1AP 


PHASE 1 


PH AS r 2 


POLNMl(CO) 




uOUTY Crjj- 

PAGfi © 


INTEGER* 2 Channel select vector from 
channels card. 

INTEGER* 2 Array of channel numbers from 
statistics deck. 

INTEGER*2 Array of channel numbers 
from channels card. 

INTEGER* 2 Array of classification pools 
and pooling information. 

INTEGER*2 Array of stacked class num- 
bers organized by pool request. 

L0GICAL*1 Logical variable signalling 
input of statistics on cards in con- 
trol card deck. 

LOGICAL*! Logical variable requesting 
production of classification map. 
LOGICAL* 1 I.ogical array indicating 
storage of calibration values by user. 
LOGICAL*! Logical variable requesting 
production of singular cell map in 
phase 2 (annexation) . 

LOGICAL*! Logical variable requesting 
initial cell processing is to be carried 
out. 

LOGICAL*! Logical variable indicating 
annexation of cells is to be carried 
out. 

L0GICAL*1 Logical array signalling that 
the name for POOL (I) has been stored. 
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PRSTAT 


S YM ( 60 ) 
SYMMTX (60) 


LOGICAL*! Logical variable indicating 
request to print statistics information. 
LOGICAL* 1 Array of user supplied symbols. 
LOGICAL*! Array of ’default symbols. 
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SECCOH 


OATE = 77129 


15017040 


PAGE 0001 


FILE SECCOM 


C SECCOH LARS XXXX IjjfgoolS 

C** ************************************************************* ****$***SUP 00030 
c SUP00040 
C SECCOM - COMMON BLOCK FOR SUPERVISED ECHO §2222222 
C WRITTEN 11/22/76 BY P.D. ALENDUFF §2£222$8 
C SUP00070 
C *********************************************************************** SUP 00080 


BLOCK DATA 
IMPLICIT INTEGER 


( A— 2 ) 


N0FET3, NOPOOL. NWCJRD, NWORDZ , NVR* CUTPUl • rULNAni <!«oU It PKfcf k A \ 1 1 1 
PTS, RQTAPE, RQFILEt SYMCNT, VARSZ3, CSEL(30), CSEL3<30), 
FETVECI30 > • FETVC3(30) » P0LPTR(2 ,60 ) , POLSTK(60>, 

CDFLAG, CL^HAP, C$ET1($,30>, OBJMAP, PHASEl, PHASE2, POLNMK 60 ) , 
PRSTAT, S YM (60 ) » SYMMTX(60) 


SUP00090 

VARIABLES USED IN SECCOM - COMMON BLOCK FOR SUPERVISED ECHO SUPOOUO 

BLOCK DATA SUP00120 

SUP00130 

IMPLICIT INTEGER * 4 (A-Z) §H£22}£2 

SUPOOI50 

COMMON /SECCOH/ BUFROZ, CSELEC, CSET(3 , 30 ) , CSET3C 3, 30 ),THRES,CEL SI Z» SUP00160 

1 CELWTH, INFO ( 17 ( « INPUT , INTTAP, INTFIL, JPTS, LINES, NOCLS, NOFET , SUP00170 

2 N0FET3,N0P00L»NW0RD, NW0RD2, NVR, OUTPUT, POLNAMI 2, 60 I, PREFIX! 2 ), SUP 001 80 

3 PTS, RQTAPE, RQFIL6, SYMCNT, VARSZ3, CSEL(30>, CSEL3I30), SUP00190 

3 FETVECI30 ) , FETVC3(30), P0LPTR(2 ,60 ) , POLSTKC60), SUP00200 

A CDFLAG, CLSHAP, CSET1(3,30), OBJMAP, PHASEl, PHASE2, POLNMK 60 ) , SUP002IO 
5 PRSTAT, S YM (60 ) , SYMMTX(60) SUP00220 

SUP00230 

REAL * 4 CSELEC, CSET, CSET3, THRES SUP00250 

INTEGER * 2 CSEL, CSEL3, FETVEC , FETVC3, POLPTR, POLSTK SUP00260 

SUP00270 

LOGICAL * 1 CDFLAG, CLSMAP, CSET1, OBJMAP, PHASEl, PHASE2, P0LNM1 ,SUP00280 
I PRSTAT SUP00290 

SUP00300 
SUP00310 

DATA SYMHTX/' I* , * 2 •» *3 • , ‘4 • , • 5*, • 6* , • 7* , • 8* , •?• , • A* , »B* , §22221? 2 

l»C»,'0», , E , ,'F*,'6*, , H , ,*I , ,*J , f , K , , , L», , M , , , N*, SUP00330 

2*0' 'P' '«• *R' 'S' 'T'.'U' SUP00340 

SUP00350 

4*1 , , , 2 , , , 3 , , , 4 , ,*5*, , 6 , , , 7 , , , 8 , , , 9*, , A , ,*B , , , C*, SUP0C360 

S'D't'E'.'F'.'G'/ SUP0Q370 

* * SUP00380 

END SUP00390 


CSELEC, 


CSET3, THRES 


INTEGER 


CSEL, CSEL3, FETVEC, FETVC3, POLPTR, POLSTK 


DATA SYMHTX/' 1',* 2', *3 ' , '4 • , • 5*, • 6* , • 7* , • 8' 
L'C','0' f 'E' t ’F','G','H','I' ( ' J' f 'KS'LS'M' 
I'0','P'»'Q',*R' f 'S'»'T','U',»V','W',*X' f »Y* 


,'N', 

,'Z', 


5'D','E','F','G*/ 


f ' “ t W f’ A f T»i.t 
<' t 

,*8S , 9 i . # A*, t fi , t , C*t 


CO 

M 

O 

O 

0 
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1 

CTl 
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1. Module Usage 
SECRDR 

Calling sequence: 

CALL SECRDR 

. SECRDR is called without parameters to read user option 

selections and detect any errors on specified combinations 
of options. Upon error detection, either corrected data is 
requested or the function is terminated following an informa- 
tion message detailing the error. Output and/or intermediate 
tapes and the multispectral image storage tape are mounted 
and a list of selected options is produced. 

2. Internal Description 

SECRDR first performs initialization on all variables 
stored m SECCOM common block to set up any default options and 

i 

clear any previously selected options. Control card decoding 
is heavily dependent upon LARSYS control card processing rou- 
tines CTLWRD and BCDVAL. CTLWRD is called with a list of ac- 
ceptable keywords to read input control cards. The number of 
the matching keyword is returned to SECRDR and a branch is made 
to the code that will appropriately decode the remainder of the 
control card. Other special purpose LARSYS routines used for 
control card decoding include CHANEL for decoding the CHANNELS 
card, and POLSCN for decoding the CLASSES card. 

After encountering a DATA or END card, a check is made for 
a request of at least one channel on the CHANNELS card. The 
cell width parameter is then checked to be greater than or equal 
to 2. The cell selection option must also be non-negative, and 
annexation threshold TIIRES must be zero or positive. If both 



SECRDR-3 


an object map and classification map are requested, only an 
object map is produced. At this point, tape and disk parame- 
ters are checked. If OPTIONS INTERMEDIATE (start processing 
with the intermediate tape) is selected, a tape number must 
be supplied. If an intermediate tape number is supplied, a 
file number must also be specified. 

Upon specification of OPTIONS INTERMEDIATE, processing 
begins with the intermediate tape and a results file is pro- 
duced. If, however, no OPTIONS INTERMEDIATE request is made 
and both an intermediate tape number and a results location is 
specified, a conflict is recognized and EXECUTION is terminated. 
If a results location is needed and none is specified, more data 
is requested. A file number must be specified when a tape loca- 
tion is selected. Specifications of both tape and disk options 
for classification results information causes termination. 

After the error checking sequence is finished, tapes are 
mounted on the appropriate devices by calls to MTAPE. Initial- 
ization, if requested, is performed on the results or interme- 
diate tape. A list of the selected options is then produced. 
Both GLOCOM and SECCOM common blocks are used by SECRDR. 

3. Input Description 

Input data is read from unit READIN in GLOCOM which is 
either CRDRDR, the card reader, or KEYBD , the console keyboard, 
when cards are typed m by the user. This selection is accom- 
plished by the LARSYS system. 

4. Output Description 

Output is written to both PRNTR, the printer, and TYPEWR, 
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an object map and classification map are requested, only an 
object map is produced. At this point, tape and disk parame- 
ters are checked. If OPTIONS INTERMEDIATE (start processing 
with the intermediate tape) is selected, a tape number must 
be supplied. If an intermediate tape number is supplied, a 
file number must also be specified. 

Upon specification of OPTIONS INTERMEDIATE, processing 
begins with the intermediate tape and a results file is pro- 
duced. If, however, no OPTIONS INTERMEDIATE request is made 
and both an intermediate tape number and a results location is 
specified, a conflict is recognized and EXECUTION is terminated. 
If a results location is needed and none is specified, more data 
is requested. A file number must be specified when a tape loca- 
tion is selected. Specifications of both tape and disk options 
for classification results information causes termination. 

After the error checking sequence is finished, tapes are 
mounted on the appropriate devices by calls to MTAPE. Initial- 
ization, if requested, is performed on the results or interme- 
diate tape. A list of the selected options is then produced. 
Both GLOCOM and SECCOM common blocks are used by SECRDR. 

3. Input Description 

Input data is read from unit READIN in GLOCOM which is 
either CRDRDR, the card reader, or KEYBD, the console keyboard, 
when cards are typed in by the user. This selection is accom- 
plished by the LARSYS system. 

4, Output Description 


Output is written to both PRNTR, the printer, and TYPEWR, 
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the console device. A variety o£ error messages is produced 
by SECRDR, and brief list follows: 

ERROR IK TAPE OK FILE SPECIFICATION-TYPE CORRECT CARD 

A non-numeric character was entered as either a tape or 
file number. 

EBROR IN PARAMETER VALUE SPECIFICATION-TYPE CORRECT CARD 

A non-numeric character was specified as the cell size, 
cell selection, or annexation value. 

YOU HAVE ENTERED X SYMBOLS. THE MAXIMUM ALLOWED IS 60. EXCESS 
SYMBOLS WILL NOT 3 E USED. 

The maximum number of symbols that can be stored is 60. 
Only the first 60 can be used. 

CELL SIZE MUST BE GREATER THAN OR EQUAL TO TWO- DEFAULT Or 2 
ASSUllED-TYPE CORRECT CARD 

The CELL SIZE (X) parameter cannot be less than 2. A 
corrected card must be supplied. 

CELL SELECTION PARAMETER MUST BE GREATER THAN OR EQUAL TO ZERO- 
TYPE CORRECT CARD 

The CELL SELECT (X) entry cannot be negative, a new card 
is requested. 

ANNEXATION THRESHOLD MUST BE NON-NEGATIVE. TYPE A CORRECTED 
ANNEXATION THRESHOLD CARD. 

The THRES specified is negative. The user is requested to 
correct the error. 

BOTH SINGULAR CELL AND CLASSIFICATION MAPS REQUESTED. ONLY 
OBJECT MAP WILL BE PRODUCED 

Only one map can be produced during annexation. If both 
are requested, only the singular cell map will be printed. 
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>F THS 
POOR 

NO INTERMEDIATE TAPE SUPPLIED FOR OPTIONS INTERMEDIATE. TYPE 
IN INTERMEDIATE CARD 

When requesting that processing start from an intermediate 
tape, an intermediate tape number must be specified. 
INTERMEDIATE FILE OR INITIALIZE MUST BE SPECIFIED. TYPE 
ADDITIONAL INTERMEDIATE CARD 

When producing an intermediate tape, some indication of 
file number must be given. 

BOTH INTERMEDIATE TAPE AND RESULTS LOCATION SPECIFIED WITHOUT 
OPTIONS INTERMEDIATE-JOB TERMINATED 

It is , impossible to determine which part of the classifi- 
cation should be performed when this set of options is 
supplied. Execution is terminated. 

EITHER INTERMEDIATE TAPE OR RESULTS LOCATION MUST BE SPECIFIED 
WITHOUT OPTIONS INTERMEDIATE. TYPE IN ADDITIONAL CARD 

Some output file is needed for processing. Additional 
information is requested. 

EITHER RESULTS FILE OR INITIALIZE MUST BE REQUESTED. TYPE IN 
ADDITIONAL RESULTS CARD 

When producing a results tape, some indication of file 
number must be given. 

NO RESULTS DESTINATION SPECIFIED-TYPE IN RESULTS CARD 

Either tape and file or disk must be specified when results 
will be produced. Additional information is requested. 

BOTH RESULTS TAPE PARAMETERS AND DISK SPECIFIED-FUNCTION 
TERMINATED 

Either tape or disk can be selected for results but not 


EOTODUe©&lTY i 
ORIGINS PAGE IS 


both. 
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BOTH FILE AND INITIALIZE REQUESTED FOR INTERMEDIATE TAPE- 
FILE REQUEST IGNORED 

The intermediate output tape is initialized, only file 
1 can be initialized. 

BOTH FILE AND INITIALIZE REQUESTED FOR RESULTS TAPE-FILE 
REQUEST IGNORED 

The results tape is initialized. Only file 1 can be 
initialized. 

In addition to these messages, the' list of options selected 
is printed on the unit PRNTR. 

5. Supplemental Information 

Common blocks GLOCOM and SECCOM are used by SECRDR. 



6 . FLOWCHART 


secrdr-7 



Qf 

s ?OOB 















gEPRQDUQB&lTY OP THB 
0B1O1NA.L PAGE IS POOR 


SECRDR-9 
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SECRDR - CONTROL CARD READER FOR SUPERVISED ECHO 


SUBROUTINE SECRDR 
IMPLICIT INTEGER * 1 (A'Z) 


“ K lfcvmnsggssr!^!TIIsSs:^ g - CLASSX - 

3 FILESY, FlDBND. HOATA, HEA0I881, 10(2001, JNACEX, 

A IXARK, KEYBD, HAPTAP, KAXCHA, HAXCLS, 

i am : s3&Pk:°iiSf.ra, WNra> Se *° ih - 

l 88MOfer* fcK!M,m ' 


REAL * A FR0CALI5,30> _ 

INTEGER • A CCWENi 1 16 1 f OATEt 5), HEOII 161, HE02U6J, TlNEISl 
INTEGER * 2 BLANK? 

LOGICAL * A CHKOUT 
LOGICAL * 1 BLANK I 

jJRRMiHHi;* immUB{liiinaUKl!tt2l»tii. 

2 CTIHE(U,HFADI5B>>, t CDHENTt I ) , HE ADt 72 )) , 

3 IMaPSAV.TPSTATII) It 

A (SEP SC R,TPSTAT(2) I , ( CUPIK.TPSTATO I J . ( C ASTAT , TPSTAT ( A 1 1 , 

$ ICOPSER.TPSTATISl J , I IRAOUT, i PSTAT(6 1 1 , 

6 (BLANK, BLANK?. BLANKl) 

CONMGN /SECCOM/ BUFROZ , CS EL EC . CS ET(3. 30 J » CSET 3 < 3, 10 1 , THRES .CEL S I Z , 

1 CELWTN.INFOt ITT. INPUT, INTTAP, INTFIL, JPTS, LINES* NOCLS. NDFET, 

2 NOFE T 3.NQPG0L »NNCR0,NWflR0? • NVR. OUTPUT, POLNAM( ♦ , 60 I, PREFIX! 2 1, 
J PTS, RQIAPE * ROFILE, SVHCNt, VAft$Z3, CSEU301. CSE13(30I, 

J FEIVECI 3 D), FETVC 3 C 30 I , P 0 LPTR( 2 , 60 l » P 0 l$TK( 6 Dl, 

A CDFLAG, CLSMAP. CSEU(i.30>, O0JPAP, PHASE1, PHASE?, P01NNII601, 

5 PRSTAT, SYM60), $YHKTXl60» 

REAL * A CSELEC, CSET, CSET3, THRES 

INTEGER • 2 CSEL, CSEL3, FEtVEC, FETVC3, POiPTR, POLSTK 


local variable definition 
real • A RV0RXC6O1 

INTEGER • A SUPLST (121 » TAPLST(A), CElLST ( 2 ) . ANNLST, OPlLST, 
I PRNLSTI3), CROLST. (WQRXI60I 

LOGICAL • l LWORKin.LSYMm 

EOUIVAtENCE IINORKID.LNORKII I.RWCRM ID, tL5Yft(l»,SYM< ll I 
OATA SLPtST/'RESU'.MNTE'.'CELL*, *ANNE a , 'OPTP ,»PRIN*. 'CLAS*, 

DATA CELLST/^SIZE'.'SELE 1 / 

OATA ANNLST/*THREV t OPTLSTF • INTE V, C*DLST?«READ*/ 

DATA PRNLSI/’SING 1 , 'CLAS*, 'STATV 
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STOPFG • 0 
OUTPUT - 0 
SYMCNT » 0 
NOPOOL « 0 
CSELEC - O. 

CELHTK • 0 
OPTJNT « O 
RUT A PE - -I 
ROFILE « -l 
RQ 1 NI 1 a 0 
INTTAP * -| 

INTER - -l 
INMNT » 0 
THRES - 0.0 
N 0 FET 3 a 0 
DO 20 1 * l. 30 
CSELUII a 6 
FET V£ 3 ( 1 ) « O 
DO 10 J * I, 3 

csEhu.n i - 50000 . 

CSETl ( J, I > a .FALSE 

10 continle 

20 CONTINUE 

DO 30 I - l, 60 
SYR 1 1 ) * SYRMTXU) 

POLSTK ID « 0 

HUNII!,! : T 8“ E - 
10 ®2i? u i .« lSE . 

PHASE 2 • .FALSE. 

COFLAG • .FALSE. 

CLSHAP - .FALSE. 

D 8 JRAP a .FALSE. 

PRSTAT • .FALSE 

C BEGIN CONTROL CARO PROCESSING 


IF ICONPUT -EQ. 21 WRITE ( TYPEWR.5Q1 
SO FOaHAie !*♦♦* TYPE IN SUPERYlSEC ECHO CONTROL CAROS I SECRDA ) 


ASSIGN 100 TO GOVEC 
10 GO TO GQV£ C. 1 100,11001 
LOO CALL CUVROlCARO.CCL.su 


100 IF L (cScE W !EO i *llI C SPiREDl T i l I ,CO,:ETRe * OIK,6,,RCO^ ' , 

IF JER 8 CDR .EQ. 1 > CALL RTMAIN 


C ERROR OK CONTROL CARD 

C no WRITE I TYPEWR . 1 1101 CARD 
9110 FORWAYf 5 X, 20 A* ) 

11 $ CALL ERPRNI ( ERRNUH, * GO * 1 
120 IF lERKCOR .NE. 31 fRRCOR - I 
CO 10 100 


c INTERMEDIATE TAPE CARD 

& 30 D ENTRY • 1 

310 15 TSZ « 2 ♦ ENTRY 

31 $ IF (COL .EQ. 721 GO TO 90 

CALL CILPRW (CARO, COL. TAPLST, 1 , 
CO 10 ( 320 , 330 , 310, 3501 , CODE 

C TAPE SPECIFICATION 


SECRDR -10 
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ooa? 

006* 

ooas 

ooa& 


0087 

ooa« 

ooa-9 


0091 

0092 

0093 
009* 


00<Jf 

0090 

0099 


m 


010*5 
01 C* 
?IC7 


0109 
01 10 
01 11 


Oil? 
01 I* 
011 $ 
0116 


0117 
01 IB 
0117 
0120 
0121 
0122 
0123 


012* 

012$ 

0126 


C ERROR OH TlP£ SPECI P (CATION 

c am iKSfxj'iitsi ass 

” 4 %We"oSI1c| S £m( 

GO TO 120 
C Flit SPECIFICATION 
C 

330 SLL Z lML|C«O.COL,lWORR,V€CS2,C32n 

!f mi :IS: ? -iBBsn} 

IF CEKTRY .EG. 2> ROFILE * 1M0RM1) 


INITIALIZE REQUEST 

345 |F UKI** *19 1 


m ® 

C RESULTS disk 

C 350 OUTPUT - C USSR 
00 TO 315 

\ CELL CARO 
C *00 


2 I ROIN 


Hl^US^^c^ellst, 

GO TO (*IQ.*201, CODE 


2, CODE* £1201 


C CELL SIZE 

M ° CALL^lVALl CARD,C0L»lWQRKrVECSZ*C*111 
CELWTK - IWGRKII) 

OU TO *00 

C ERROR 
C 

*11 WRITE (TYPEWR,9*I1J CARO 

9*11 FORMAT t/5x!2iA*//7X.*6*** ERROR IN PARAMETER VALUE SPECIF1CATII 
1 - TYPE CORRECT CARO ISECRORl'l 
CO TO 120 

£ CELL SELECTION PARM 
C 

* 2 ° cALL Z FVAL[CARO,COL t RMOR>£.VECSZ,£*lll 
CSELEC « RhORKil) 

CO TO *00 

£ ANNEXATION CARO 

c 

500 


SUPO 1530 

suP0i5*r 

$UP01S5< 

;uPor 

jUPO ' 
SUPO. 
SUPO 
SUPO 
SUPO 
SUPO 
SUPO 
SUPO 
5UP0 
SUPOI 
SUPO 
SUPO 
SUPO 
SUPO: 
SUPO 
SUPO 
supo 

,UPO 
UPO 
UPO 
^UPO 
SUPO 

supo: 

SUPO 
SUPO 

supo: 
supo: 
supo: 
supo 

SUPO 
SUPJ" 
SUPL- 
super 
supoi 

SUPOI 
SUPO- 


m 

-610 

620 

630 

6*0 


Hi 

760 

790 

600 

aio 


900 

910 

920 

M 

-960 

970 

960 

888 


IF (COL EQ. 72) CO TO 90 

CALL CTLPRM( CARD* COL, ANNlSf, 1 , CODE, Cl 20) 

CALlS VAL (CARO i COL, R WORK , VECSZ , £* 1 1 1 
IF (VECSZ *EQ. 01 GO TO *11 
THRES « RVCRKIU 


GO TO 500 
C OPTIONS CARD 

C 600 IF (COL EO 721 GO TO 90 

CALL CTLFRRICARD.COL.OPTLST,!, 

optint • 1 


Supo 
supo, 
supo, 

SUPOI 
SUPO! 
SUPOI 
SUPO^ 
SUP02010 
SUP02020 
SUP0203C 

iOTS§ 

SUP020bO 
SUP02070 
SUP02080 
SUP02090 
SUP02100 
SUPD211C 
SUP02120 
SUPO2130 
SUP02 I* 0 
SUP021$0 
SUP02I60 
SUP02170 
SUPO210O 
SUP02I90 
SUP02200 
SUP02210 
SUP02220 
SUPD22J0 
SUP 022*0 
SUPD22S0 
SUPD2260 
SUP022 /O 
SUP02280 
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DATE ■ 77129 


FILE SECROR 


01 28 
0129 
0120 

0131 

0132 
OL33 


013* 

0135 


8 H? 

01*0 

01*1 


01*2 

01*3 

01** 

01*5 

01*6 

01*7 

01*0 

01*9 

sii? 


0152 

HR 


0167 

0163 

0169 

0170 
01 71 


c GO TO 600 

C PRINT CARO 

700 TF (COL EQ. 721 GO TO 90 _ 

CALL CTlPJtK (CARD.COL.PRNLST, 3.C00E.C 1201 
IF ICOOE * TO 11 CSJMflP * .TRUE. 

IF (CODE .EO. 2) CLSMAP * .TRUE. 

IF (CODE .EO 31 PRSTJ-- 

c GO TO 700 

C CLASSES CARO 


STAT * .TRUE. 


000 CALL PGLSCN IPOLNAH. POLPTRrPOLSTK.POLNMl, NOPOOL rSTXPTR, COL, CODE, 
I CB10I 
GO TO 90 

ERROR ON CLASSES CARO 


CARDS READSTATS 
050 IF 


IF (COL .EO. 721 GO TO 90 


,1 »CC0E,C1201 


GO TO B50 
' SYMBOLS CARO 


ICOVAUCAR 

yi 6 oJ 

SYKCNWI 


SYMCNT 
GO TO 90 
9*0 ERRNUM * **T 
GO TO 110 

C TOO MANY SYMBOLS SPECIFIED 
C 

9S0 WRITE 1TYPEVR .99501 K 


WRITE (TYPEVR, 99501* 
WR| TEfPRNTR .99501 K 


ID*. 13,' SYMBOLS 
U MOT BE USED 


CSETMI.J) « TRUE. 


ID IS 60-'/9X.» 

I * I - (K-601 
GO TO 502 

C CHANNELS CARD 

C 1000 CALL CHANELfCARD,C0L,NCR.CSEL3.C$ET3,FETVCJ,tl07O> 
00 106C X * 1,3 
GO 1060 J < 1,30 
1060 IF ICSET3(lrJJ NE. -50000 

artiM " ci 

c 107 '’ !§“TiIo 111 

C DATA OR END CARO 

l E0. II CALL RTMA1N 
C CHECK IF NUMBER OF CHAN GT 0 

C IF INQFET3 GT. 01 GO TO 1120 

ERRNUM » 3*3 
ASSIGN IloO TC GOVEC 
ERR COR - 3 
GO TO 11S 

C CELL SIZE GT 0 IF SPECfFIEO 


. THE maximum ALL0UESUPO27SO 
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P^iljtY OP THE 

is poor 


1170 IF IC6LWTH »GE. 2) CO TO IUS 
CELHIK - 2 
Mill. |TYPLMa.9i 1701 


Mill. | TYPLUa.Tl 170) StlPOJOBO 

Mill (PRN1R.Q1 |20) SUP0309C 

’•"".sersf’i iiiLF'-uiiihJ'iY 86 “ tA, “ 08 EOU ‘ l ,0 * - otF ‘^o o ;ifs 

1124 SiMCH l LlO0 TO CDVEC fuposijo 

ERR CO A • 1 SUPO3i*0 

GO ID 100 5UP03150 

C 5UP03160 

C BOTH MAPS SUP0317Q 

C _ SUPO310O 

U*5 IF l-NCT. CLSMAP OR. .MOT DBJHAP) GO TO 1150 SUP03190 

MITE mP£WR,911*5l SUP0320C 

MITE t PAMIR, 911*$ I SUPO32I0 

91U5 FORMAT I » E*** BOTH SI NGULAR CELL AMO CLASSIFICATION MAPS REQUESTf SUPO3220 
10 - ONLY SINGULAR CELL MAP»/9X**W ILL BE PACOUCEO (SECRDRl') 5UP0323C 

CLSMAP > .FALSE. SUP032*C 

C SUP 03250 


C BEGIN TAPE PARAMETER CHECKING FOR RESULTS AND INTERMEDIATE TAPE 


AND CUTPUT.EO.O I PHASE I 


,lM Iflftllff :& S , .«S!&;a , ! 1!§: a,™. 

l |MOP?INT.NE O.QR (RQTAPE.EQ -l.AND OUTPUT-EO Oil 
1 CO TO 1151 
PHASE 1 - .tRUE. 

PHASE? - .TRUE 
1151 CONTINUE 


ir Ik.lCLCk. .Lit. V 

MITE ITYPEMR » 91 1 
1*R. I IE (PRNTR ,9113' 
1130 FORMAT ( ■ £♦** C£ 


IL SELECTION PARAMETER MUST BE GREATER THAN OR */SUP03*6C 


TYPE CORRECT CARD 


C CHECK FOR ANN GT 0 

1 l AO IF 1THRES GT. 0.0 .OR .NQT.PHASE2) GO TO U52 
kRITEI IYPEWR.911A0) 

M|TEtPRNTR,9H*0) 

9 H AO FORMAT!' £♦•* ANNEXATION THRESHOLD MUST BE NOM-NEGAT 1VE. » 

1 /8X. 'TYPE A CORRECTED ANNEXATION THRESHOLD CARO (SECRDRl* I 


MlTEItYPEM. 911501 
MRITEI PRNTR, 911501 

91150 FORMAT I* (-«** NC INTERMEDIATE TAPE SUPPLIED FO', 

1'* OPTIONS INTERMEDIATE.' /8X, 'TYPE IN INTERMEDIATE CARD 

2 1 * I 

CO TC 1125 


CHECK FOR INTERMEDIATE FILE OR INIT 


VRI IE1 PRNTR.91 1551 
Ml TEI TYPE M.9 1155) 


»Kt It; I ITKCHK,YIK73J 

91155 fOKMAll* £*•* INTERMEDIATE FILE CR INITIALIZE MU$t BE SPECIFIED 
1 /fix* • TYPE ADDITIONAL INTERMEDIATE CARD ISECRDR 1* I 


SUP03WC 
SLP03A80 
$0203*90 
SUP0350D 
SUP03S10 
5UP03520 
SUP0353C 
SUP035AC 
SUP0355C 
* I SUPO 356 C 
5UP0357C 
SUP035HO 
SUP03590 
SUP03600 
SUP03610 
SUP0362C 
SUP0363C 
SUP03650 
SUP03650 
ISECRDRSUP036&C 
SUP0367C 
SUP036SJ0 
SUP03S90 
SUPO3700 
SUP037I0 
SUP0 3720 
SUP03T30 
$UP037*0 
IFIEO 'SUP0375C 
SUP037SC 
SUP 03770 
SUP0379 0 
SUP03790 
SUP0380Q 


1560 IF I. NOT. PHASE 1 .OR. INTFIL €0- -1 .OR. INTINT £0. 0) SUPO301C 

I CO TO 1600 „ SUP039ZG 

MITE! TYPE MR.91S60I SUP03630 

kRlTE(PRNTR t 915&0) SUPO38A0 

91560 FORMAT I* IC&OO BOTH FILE AND INITIALIZE REQUESTED FOR INTERMEDI ATSUP038SC 
IE TAPE -V8X.«FILE REQUEST IGNOAEC ISECRORIM SUP03B60 

INTFIL * -1 SUPO307O 

C SUPO3B90 

C CHECK FOR BOTH XNT AND RESULTS 5UP03B90 

C 1600 IFIIINTTAP.EQ.-1.CR I RQTAPE.EC.-I .AND OUTPUT EQ 01 1 .OR .CPTINT NE* SUPoll ?0 


INTFIL €0- -1 .OR. INTINT £0. 0) 


Ml (IfClKtYlnU 
kRl IE I PRNTR.91560) 

1 FORMAT f* IC&OO 80 


C CHECK FOR BOTH XNT AND RESULTS 5UP03A90 

C 160 0 IFIIINTTAP.EQ.-1.CR I RQTAPE.EQ.-I .AND OUTPUT EQ 01 1 .OR .CPTINT NE* SUPoll ?0 

l 01 GO 10 2100 5UP03920 

kRltE I TYPEtaft, 920001 SUP03930 

MITE TPftNTR. 92000 1 SLIP 039*0 

92000 FORMAT I* £**• BOTH INTERMEDIATE TAPE AND RESULTS LOCATION SPECIF ISUP03950 
1ED WITFDUI CPT10N5V9X,' INTERMEDIATE - JOB TERMINATED I SECROR > 'SUP03960 
21 SUP0397C 

CALL RTMAIN SUPO390O 

C SUPO3990 

C SEE IF EITHER TAPE WAS REQUESTED SUP 0*000 


KKl It IPKNIH.9ZIUUI iUP0*050 

92100 FORMAT! ■ £•*• EITHER INTERMEDIATE TAPE OR RESULTS LOCATION MUST BSUPO*OfiO 
IE SPECIFIED W|TH0UT*/6X* 'OPTIONS INTERMEDIATE - TYPE IN AEOI TI0NASUP04Q70 


t CHECK FOR RESULTS If NEEDED SUP04110 

C „ . SUP04120 

2200 IF l.NOT. PHASES) GO TO *000 $UPO*130 

IF (RQIAPE.NE.-1. OR. OUTPUT. NE 01 GO TO 3050 SUP0*1*0 

WRITE! TYPEkR, 930001 SUP0415C 

MITE! IRNTR. 93000) SUP04160 

93000 FORMAT I * E*** NO RESULTS DESTINATION SPECIFIED - TYPE IN RESULTS SUPQ4170 

I CARO .. ( SECRDRl '1 SUPO* 180 

CD TD 1125 SUPO* 190 


CD TD 1125 SUP0A190 

C SUP04200 

C CHECK FOR FILE AND INIT SUPO*?10 

C SUP0*220 

3050 IF CROFILE .EQ- -* .AND. RQlNlt .EQ 01 CO 10 3075 SUPO*23C 

SUP 0*2*0 

WRITE I PRNTR, 930501 SUP042S0 

93050 FORMAT !• £*•• EITHER RESULTS FILE OR INITIALIZE MUSI BE* , SUP04260 

1 * REQUESTED •/«X*»TYPE IN AD01TICNAL RESULTS CARD ISe£rOA1*1 SUP042TC 
GO TO 1125 SUP0428Q 

C SUP0*290 

C CHECK FOR BOTH FILE AND INIT SPECIFIED ON RESULTS REQUEST SUP04300 

C SUP04310 

3075 IF IRQFILE -EQ. -1 .OR. RQINIT .EC. 0) CO TO 3100 $UPO*320 


; CHECK FOR TAPE PARAM AND DISK 


■ ic irnniR.vjtwi 

93100 FORMAT ( • E«** 80TH RESULTS TAPE PARAMETERS ANO DISK SPECIFIED - 
UNCTION TERMINATED ISECRDR)*! 

CALL RTKAlN 

*000 IF (OUTPUT .EQ D) OUTPUT - KAPT*P 
IF ICHK0UT1 RETURN 


C MOUNT TAPES - result 


call MTAPE (ROT APE.ROFILE* 1 ) 
IF IRQFILE .CE. 0) GO TO *20i 
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DATE » 77129 


file secrdr 

0256 

0237 

0738 

0739 


0260 


oil 


0262 

0263 

026* 

0265 

0266 
0267 


*100 CHRtDR . i 

ASSIGN UOO in GDVEC 

cu la too 

£ INTLRKEOIAIE TAPE 

*200 IF (PHASE J -AND PHASES) GO TO *300 
IF IPHASE2 ) NODE * -2 

IF (INtlWI .NE. O .AND. PHASEl) INTFU « 0 
CALL HTAP£tlHttAP,INTF!L,*ODE) 
IFIINTFIL.CE.OI GO TO *300 


C OPTIONS 

C *300 IF UNTFIL .£0- OJ INTFJL 
IF (RQFILS -£Q. 01 RQFILE 

94300 ^ 


SELECTED THE FOLLOWING SUPERS ISIO ECHO OPTIONS 


94301 
9*302 ^ORRAT 


CONTROL CARO DECK* I 
ICAH^H RESULTS NAP* 1 



?*JU2 MJK**I I ’PHI NI Lt ASS l r 11AI 1 1JH KtS 

94303 F0RRATCT5. 'PRINI^SINGUlIr^CELL NAP* 

94104 FORMAT*?!, UeRF0RN P In[tIal 3 ?ELL PROLts: 

(F IPHASE2) WRITEIPRNTR f 9410SJ 
94305 FORMAT IT5r 'PERFORM CELL ANNEXATION AND . . 

94306 l ^ S b?lK!pSW? £ «!!fis R fIi 3 TgSJ STATISTICS USED* 


PROCESSING OF AREA*) 
PRODUCE CLASS1 


im 

5UP045IO 
SUP04B20 
SUP0AB3G 

fuSot'iO 

SUP0*«6C 
SUP Q* B7C 
SUPO40BO 
F (CATION RESCSUP0499C 
'UP 0*900 
UP0*9i0 
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MODULE ABSTRACT 

SECINT carries out required initialization of the rest of 
variables used by SECSUP and finishes reading any data cards 
as well as the statistics to be used with the processor. In 
addition to the array allocation performed, some option in- 
formation is printed and a loop is entered that carries out 
the desired stages of classification for each area to be 
classified. 
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1. Module Usage 
SECINT 

CALL SECINT 

The program SECINT is called with no arguments . Any 
variables to be used or changed are contained in common blocks 
GLOCOM and SECCOM, both required by SECINT. All information 
required to perform the classification is gathered by SECINT 
from the appropriate source (ie., statistics and areas to be 
classified) and array allocation is carried out using va- 
riables supplied to common block SECCOM by the control card 
reader SECRDR. Statistics are reduced by calling REDSAV to 
save statistics for only classes and channels requested and 
statistics information is printed by calling SECPRT. Initial 
records are then written to the output device. A loop is 
entered and, for each area to be classified, a sequence of 
tests is performed. If multispectral data is required, the 
correct run is requested and existance of requested channels 
is confirmed. Further allocation is performed and SECH01 is 
called to perform the needed processing for the area. 

2. Internal Description 

The first information handled by SECINT is the statistics 
file. If processing is starting from an intermediate tape, 
records 1 and 2 are read from the intermediate tape to re- 
establish information used in initial processing such as chan- 
nels selected and calibration used, and classes selected and 
pooling requested. The statistics file is then read from the 
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tape and written onto the unit SDATA, the statistics file on 
disk. The next card is checked for an "end" card and if any 
area description cards are required, the function is ter- 
minated. 

Otherwise, CLASSX, the storage disk file for all 
classification areas for this file, is rewound and processing 
loops until all areas are stored on disk. Some initial allo- 
cation is then performed as array storage is set up for the 
covariance and mean matrices for the entire statistics file on 
disk. REDSTA transfers the file into memory. CLSCHK and 
FETCHK then check for fatal errors in selection of pooling and 
channels. The addresses of reduced arrays for storing only 
needed information are then computed and REDSAV reduces the 
information to that required. 

Initial records are then produced and written *to OUTPUT. 
The entire statistics -file from SDATA is then transfered to 
the output file, and SECPRT is called to print the statistics 
information, if requested, and to produce the record type 4 
with covariance and means matrices of the classes. A loop is 
then started and executed once for each area to be classified. 
If processing begins with the data type, the correct input run 
is requested and calibration information is established. If 
processing begins with an intermediate type, all areas on the 
tape in this file are used. Field size and the allocation of 
all other arrays used by SECH01 are calculated. Then a record 
type 5, area identification record, is produced and SECH01 is 
called once for each area to be processed. When finished, the 
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final record is written to the output device. 

J. Input Description 

Input cell information for SECINT is read from CPYOUT if 
an intermediate tape is requested as the starting point of 
processing. The format is identical to that of a classification 
results file. SDATA, the statistics file on disk is also read 
to be copied into the output file. CLASSX, the areas to be 
classified as stored on disk, is also used to store the list 
of areas to be processed. 

4 . Output Description 

A copy of the statistics file is written to SDATA from 
the intermediate tape is used as input. CLASSX is written 
as the list of areas to be classified. The PRNTR is used for 
messages and output information. The file OUTPUT is the re- 
sult of processing with the same format as a classification 
results file. 

5. Supplemental Information 

This program uses common blocks GLOCOM and SECCOM. See 
LARSYS Systems Manual for classification results file format. 
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FILE SECINT 


0001 

0002 

0003 


OQ04 

0005 

0006 
00C7 
00C8 

0009 

0010 


oou 


0012 

0013 

0014 


0015 

0016 
0017 

0016 

0019 

0020 


0021 

0022 

0023 


C SECINT LARS XXXX SECOCOIC 
c SEC00020 

C*44* ••*♦*******•*****•*♦******♦*♦♦*♦♦♦♦*♦***** ♦•*•♦*♦♦♦*•♦♦*•♦*•******* SEC CCC3 0 
C SEC CCO<* 0 
C PERFORM ECHO SUPERVISED VERSION ON AREA SEC0005C 
C SEC0006Q 
G***4 +**♦«* »«*♦♦♦**»**»* »***#**+***♦*♦* ****#*♦♦*#*♦♦#♦♦***♦♦♦#♦*#♦♦****♦ SEC OC 070 
C SEC0C080 
C WRITTEN 03/4/77 BY P.D. ALENDUFF SECOC C9C 
C SEC00100 


SUBROUTINE SECINT 
IMPLICIT INTEGER * 4 <A-ZI 


COMMON /GLCCOM/ BL ANK, CARD120). CHKDUT « COPFIL. CL 
CLLST5W CONPUT , CPYOUT, CRDRDR, CRDSEQ, OATAPEt 
DUPUPf DUPRUN, ERRMSGt FBPNT * 

FILESV, FLDBND, HDAT A, HE AD( 83 ) , 10(200), IMAGEXt 
I MARK * KEYBD* MAPTAP, MAXCHA, MAXCLS, 

PAGSlit PNCH. POINT, PRESUX. PRNTK, READIN? 
RESIRT, RUNFIL, RUNT AB l 10 ,3) , 

SOATA, SEPARX, SEPTPX, SPARE(10l t TEMPASOO), 

TP STAT { 6 ) , TTFLDX , T YPEWR , 

TOP, ARRAY ( 12500 ) 


CLASSRf CLASSX , 


SEC00120 
SEC 0013 0 
SEC 001 4 0 
SEC00150 
SEC00160 
SECCC170 
SEC 00 1 80 
SEC 00190 
SEC00200 
SEC 0 02 1 0 
SEC 0022 0 
SEC 0C2 30 
SEC 0024 0 
SEC0C25C 
SEC00260 
SEC002 7 0 
SEC0028 0 
SEC0C29C 
SEC00 30 0 
SEC 0031 0 
SEC 00 0 

SECCC3 50 


c 

c 


REAL ♦ 8 ARRAY 

INrEGER i ** F 4 0 caMENf?lM • DATE<5) # HEDi (161, HED2(16)* T1MEI5J 
INTEGER * 2 BLANK2 
LOGICAL ♦ 4 CHKCUT 
LOGICAL * 1 BLANK1 

EOUI VALENCE ( CATSAV . I D ( 1 )) , < CURRUN, I CI3I ) ♦ ( FROC AL ( 1 ) • ID ( 5 1 ) ) , 

1 IHLD1 l 1 ) ,HEAO (8H . <0ATEU>fHEAD(26>), (HED2( lit HEAD! 39M.SEC00340 

2 ( T IM£ ( 1 ) , HEAD <58)1, ICOHLMl i ) ,HEAO( 72 ) ) , SEC 00 35 0 

3 IMAPSAVtTPSTATI lMf SEC0C36C 

4 ISEPSCR ,TPSTAT12> ) , t DUPI N ,T PST AT 13)). I DASTAT, TPSTAT ( 4 ) ) , SEC 0 C3 70 

5 (CCPSER,TPSTAT(5) ), ( TR AOUT, T PST AT (6 ) J t SEC00380 

6 (BLANK, 0LANK2,BLANKl) SEC00390 

SEC 00400 
SEC 004 1 0 

COMMON /SECCCM/ BUFROZ,CSELEC t CSETI3,30), CSET3t3,30), THRES ,CELSI 2 , SEC C042 0 

1 CELWTH,fNFO( 17) t INPUT, INTTAP, I NTF I L , JPTS, LINES, NOCLSt NOFET » SEC 0C4 3 0 

2 NOFc 13,N0PG0L,NRCRD, NW0RD2 , NVR, OUTPUT, POLN AW I 2 ,60 ) , PREF I X ( 2 ) t S EC 0C44 0 

3 PTS, RQTAPE, RQFILE, SYMCNT, VARSZ3 # CSELOO), CSE13(30>» SEC0C45C 

3 FE TVEC (30 ) * FETVC3130), PCLPTR I 2 ,60 f , P0L5TK I 60 ) t SEC C 0460 

4 CDFLAG, CL§MAP, CSETHJ*30>, OBJWAPt PHASE1, PHA5E2, POLNMK 60 ) • S6CCC470 

5 PR ST A!, SYM160). SYWNTX<60> SEC CC480 

SEC0049 0 

REAL ♦ 4 CSELECt CSET , CSET3, THRES SEC00500 

SEC0051 C 

INTEGER ♦ 2 CSELt CSEL3 , FE7VEC, FETVC3* PCLPTRt POLSTK SECCC520 

SEC C0530 

LOGICAL « 1 CDFLAG, CLSMAP, CSET1* 08 JMAP, PHASE1, PHASE2, P0LNM1 ? SEC 00S40 
1 PR S TAT SEC00550 

SEC 00560 
SEC 0057 0 
SEC 0058 0 
SEC005?0 
SEC 00600 
SEC 006 i 0 
SEC0062 0 
SEC 00630 
SECCC640 
SEC0C66O 
SEC0C66C 
SEC 0 06 70 
SEC0C68C 
SECCC69C 
SEC00700 
SEC00710 
SEC00720 
SEC0C73 0 
S6C 00 74 0 
SEC00750 
SEC00760 


LOCAL VARIABLE DEFINITION 

REAL * B BUFFER 1 6250 ) 

RFAL*4 WORK ( 60 ) 

INTEGER * 4 CLS0AT<5)» CARDK 19) , KEPPTSI 60 ) ,NCCNT ( 60 1, 

i i work m 

INTEGER ♦ 2 BINFLG 
LOGICAL ♦ 1 CHKFLG 

EQUIVALENCE ( WORK (1 1 , NOCNT 11 1 >,<CARD( U, CARDK 111, 

1 ( I WORK ( I) , WORK (II) 

DATA I EOS/ *E05 # /t€NO/ f END •/ 


S8 


INPUT ■ 
areafg 


CPYOUT 

■ 0 
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FILE SECINT 

002 * 

0025 

0026 


0027 

0028 
0029 


0030 

0031 


0032 

0033 


003* 
0035 
00 36 

0037 

0038 

0039 

0040 
00*1 
0042 
00*3 
00 ** 


00*5 

00*6 

00*7 

00*8 

00*9 

0050 

0051 

0052 

0053 


005* 


0055 

0056 
005 7 

0058 

0059 


0060 

0061 

0062 

0063 

006* 

0065 

0066 


0067 

0068 
0069 


CALL GTSERL ( SERIAL ) 

IF ( PHA SE 1 ) CO TO 200 
IF ( CHKOUT ) GO TO 300 

REAP STATISTICS TO DISK FROM I NTE RMEO I ATE -T AP £ 


READ (INPUT* END -9000)1.1*1*1 * I*!*SERIAL*MOOFLG 
IF IWOOFLG . NE . 0) GO TO 105 

READ < INPUT, END^COO) I ,1 , NOCt S* , NCFE T 3 * NGF LO , NGPCOL 

t (FETVC3- ‘ ' * 

2<FRQCAL 


NPUT*END a 9C00) I * I •N0CLSA*NCFET3*NGFL0 jNCPCDL, 

3(1 )* I *2. 1 NCFE T3 1 ,{CSEL3(£e TVC3 (111*1= I* NCFE T 31* ^ 

2lfRULALU,F£TVC3U ) ) . * NCFET3 J # < FRQCAL < 2 , FET VC3 i I ) ) * 1*1* N0FET3I • 

J UP CLNAM i, J l ♦ l=*l*2l»4*l*NQPCGLi * D PCLPTRI l • J ) * 1*1 * 2 > * 4*1 .MOPOQL ) 
**( POLS TK < I J ,1*1 t NCCLS*)t (ilSOAT( I ) * I * 3 • 5 ) , CLSD AT ( D •CLSDAT 12) 



GO TO no 

105 READ ( UPUT f ENO = 900OH*I,N0CLS*,NCFET3*NOFLD,NOPOOLf 

1 CFETVC3i‘ ‘ ‘ 

2IFR0CAUI 
3 « I 4 PLLNAM 

*?(PCL$lK( t » »i*i tnuvL 
5 CLSDAT(l) ,CL5DAT(2) 

CELWIH * tWORKU) 

IHRES * WORK ( 2 ) 

C 

C READY TO CCPY FROM INPUT TAPE TO DISK 

L 110 REWIND SDATA 

P6ADI INPUT ,END=9000) I , I »CAR01 tCRDSEQ 
BINFLG * 0 

IF IUR01U1) . 60. IONE) QINFLG « 1 
WRIIEI SDATA, *110) CARD1 tCRDSEQ 
9110 FORMAT! 18A*. IS) 

115 REACH INPUT ,END* 9 QO 0 ) I * I • CARD 1 tCRDSEQ 
FRI TE( <DATA, 9 ll 5 )CARDl*CRDSEQ 
9115 FORMAT (L3A*,ia> 

IF (CARDUD.NE. IEOS) GO TO 115 
GO TO 300 

£ CHECK STATS ON OISK 

^ 200 CHKFLG * .FALSE. 

UNIT = CRDRDR 

IF t.NCT. COFLAG) UNIT * SDA7A 

IF (UNIT . EG. CRDRDR. OR.. NOT. CHK0UT1 GO TO 210 
REWIND SDATA 
RE ADI SCA TA , 92 00 ) 1 
9200 FORMAT (20A*) 

IF (I.NE. IEOS) CO TO 210 
CHKFLG*. TRUE. 

£ STATS ON CARDS OR OISK 
C 

210 CALL SI AT I UNI T * NOCLS *NOFL0*NDFET * NUN* QINFLG) 

C STATS ARE ROW CN DISK » TRANSFER TO CORE 
C 

300 REWIND CLASSX 

IF IPHASE2-AND..N0T. PHASED AREAFG m l 
IF (CAftDIU .NE. END) GO TO 310 
IF (PHASED CALL ERPRNT ( 25* « • STOP* ) 

GO TO 350 

£ READ AREA CARDS 
C 

310 IF (AREAFG .60. I) GO TO 350 

CALL l ahead I INFO* CODE* RE AO IN *0*01 
IF (CODE .EQ. *) GO TO 350 

IF (CODE .NE. 0) CALL ERPRNT 1 373* ’(GOTO* 1 300 1 
IF (CHKOUT) GO TO 310 
WR I TE (CLASSX I INFO 
GO TO 310 
C 

£ CONTINUE PROCESSING 

C 350 VARSII « N0FET6 4 NOFET +1 ) /2 
CO VAR 1 • 1 

AVAR1 - C0VAA1 ♦ < VARSU*NOCLS*D /2 


SEC 0077 0 
SEC00780 
SECOC79C 
SEC 00800 
SECCC01C 
$ EC 00 82 0 
SEC00830 
SEC 008* 0 
SECC085C 
SEC 0086 0 
SECC087C 
SEC00880 
SECCC890 
SEC00900 
SEC 009 1 0 
SEC00920 
SEC 0093 C 
SEC 009*0 
SEC 0095 0 
SEC00960 
SEC00970 
SEC0098 0 
SEC 0099 0 
SECOiCOO 
SEC0101 0 
SECO 1020 
SECO 1030 
SEC 0 1 0* 0 
SEC01C50 
SEC 0 1 060 
SEC 0 1 0 70 
SECOIC80 
SEC01C90 
SEC 01100 
SEC01U0 
SECO 1120 
SECO 1130 
SECGl 1*0 
S6C01150 
SECO 1 1 60 
SECO 1 170 
SEC 01180 
SEC 0 1 1 9 0 
SEC012C0 
SECO 1 210 
SECO 1220 
SECO 1 230 
SEC 0 1 2* C 
SEC 01250 
SEC0126Q 
SECO 1270 
SEC 01280 
SEC0129C 
SECO I 300 
SEC 0 1 3 1 0 
SEC 01 32 0 
SECO 1 33 0 
SECO 13*0 
SEC013S0 
SEC 0 i 360 
SEC 01370 
SEC 0 I 36 0 
SECO l 390 
SEC 01*00 
SEC 01*10 
SECO 1 *? 0 
SEC 0 1 * 3 0 
SECO l A* C 
SEC0145C 
SECO 2 *60 
SECO 2 *70 
SECO 1*80 
5EC01*90 
SEC01500 
SEC01510 
SEC0152Q 
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FILE SECINT 


CL S ID 1 * AVAR1 ♦ (NOFET * NOCLS ♦ l)/2 
CAT0AS - CLSID1 ♦ NOCLS 
CORE * TOP - ( DATBAS *8 ) 

IF ICO&E *LE. 0) CALL ERPRNT t 371. * STO PM 

CALL REOSTAiARRAYlCOVARl ) * ARRAY ( AVAR! ) , ARR AY ( CIS 1 0 1) * KEPPTS* 

IB INFLG * NOC NT *NOCLS»CSET,CSEL» FETV EC, V ARS T 2 t NOFET > 

C STATS ARE IN CCAE - CHECK CLASSES AND CHANNELS IF NEEDED 
C 

C IF HE JUST GOT THEM OFF TAPE THEN THEY ARE OK 

C IF (PHAS£2*AND«*NCT*PHASE1) GO TO 380 

NO = 0 

CALL°CL$CHKt ARRAY ! CL S1DI > ,NOC LS, NCPCQL, POLPTR, POLSTK, POLNANt POLMN 

AL L^FE fcHK ( MAXCHA,FETVEC»NOFET *NCFET3» FETVC3, CSEL3»CSEL ) 

IF (N0FET3 *LE . 0) CALL ERPRNT ( 3 72 ♦ • STOP* j 
N * POLPIR !2,N0P0CL) ♦ POLPTR ( 1 , NOPOOL I - 1 
DG 360 1 1 - l >N 

360 N0FLD3 - NCFLD3 ♦ NOCNT ( POLSTK ( II )) 

C ALLOCATE STORAGE AND PRINT DATA 
C 

300 VARSZ3 * NGFET3*(N0FET3M)/2 
C0VAR3 * C C VAR 1 

AVAR3 = C0VAR3 t VAR SZ3 *NC POCL + 1 >/2 
SCR AR 3 * AYAR3 ♦ ( NOF ET3*NCP0CL ♦ li/2 
CETAR3 - SCRAR3 ♦ (NOPOOLMI/2 
C0NAR3 * DETAR3 ♦ (NCP00LMJ/2 
DATBS3 * CCNAR3 * !NCPC0t>l>/2 
END * CATBS3 ♦ ( NOFE T 3 * < NOFE T3«-2 ) ♦ 1 )/ 2 
CORE * TOP - END *8 

IF I CORE . LT.O ) CALL ERPRNT! 371 STOP* ) 

IF (.NOT* PHASED CALL SECRO ( ARRAY i C0VAR3 > , ARR AYC AVAR3 )) 

C 

C PRINT SUPERVISOR INFORMATION 
C 

WRITE! PRNTR, 93501 

9350 FORMAT!//* SUPERVISED ECHO INFORMATION***-*/! 

WRI TE(PRN TR, 935 1) SERIAL, NQPCOL, NO fET3,N0FLC 
S3SI FORMAT !5X, 'CLASSIFICATION STUDY SI9, 

1 /5X » • NO* OF POOLS *,I4, 

2 /5X , *NC* OF CHANNELS **NMt 

3 /5X , * NO* OF TRAINING FIELDS **H) 

C REDUCE STATISTICS 

C CALL REDSA V( ARRAY (C0VAR1) r ARRAY! AVAR1 I , ARR AY ( C0VAR3 ) • ARRAY! AVARS) 

1 KEPP TS* NOPOOL fNOFETt VARS 1 2? NOFET 3fPOLPTR»POLSTKtFETVECfFETVC3» 

2 VARS 23) 

£ WRITE OUT INITIAL RECORDS 

C IF !CHKOUTI RETURN 

DO 390 1*1. NOPOOL 
390 WORK ( I ) x 0*0 
II * 3 
JJ a 0 
KK « I 

I * inttap 

J = INIFIL 

PRFFIXU) » l 

PREF I X ! 2 ) * 0 

IF (PHASE2) I * ROTAPE 

IF ( PHASE2 ) J « ROFILE 

WRITE (OUTPUT) PREF l X , I , J, I I , JJ,SERI AL, KK , ( JJ, 

PREFIX!!) * 2 . 

IF { *NCT • PHASE2) 1WQRK ! 1 ) * CELWTH 

IF (.NOT. PHASE2) W0RK(2) = THRES 

WRITE (OUTPUT) PREF I X , NOCLS .NCFET3 ,N0FLD. NCPCOL, 

I (FETVC3U ).I«1 .N0FET3) • (CSEL3(FETVC3U) )lIxl t N0F€T3)t, . 

2(FRCCAL(1.FETVC3(I ) ) , I- l , N0FET3 ) . ( FRQCAL ( 2 . FETVC3 (I ) )* I«1*N0FET3) 
3# | f POLNAH (I,J)«I*1*2 J»1 iNQPOQL.1 * ! ( PQUPTR! I iJ), 1*1*2 4* i* NOPOOL 


SEC0153C 
SEC01540 
SEC0155 0 
SEC C 1 560 
SEC 01 57 C 
SEC0158C 
SEC01590 
SECC16C0 
SEC01610 
SEC 0 1 62 0 
SEC 0 1 63 0 
SEC 0 1 64 C 
SEC01650 
SEC0166Q 
I S ECO 1670 
SEC 0 1 68 C 
SEC 0 i 690 
SEC01 700 
SEC 0 1 710 
SEC01720 
SECO1730 
S EC 0 1 74 0 
5 ECO 1 7 SO 
SEC 01760 
SEC0177Q 
SEC 01780 
SEC 01790 
secoiaco 
SECOiaio 
SEC01820 
SEC01830 
SEC 0 l 84 0 
SEC 0 1 85 0 
SEC0I860 
SEC 0 1 0 7 0 
SECC18B0 
SEC 01890 
SEC 0 190 0 
SEC019I0 
SEC 0 1 92 C 
SEC 01 93 C 
SEC01940 
SEC 01950 
SEC01960 
SEC 0 1 9 70 
SEC01980 
SEC 01990 
SEC 02 COO 
, SEC 02010 
SEC02020 
SEC 02030 
SEC02040 
SEC 02 05 0 
SEC 02C60 
5 EC 02 070 
SEC02CaO 
SEC02C90 
SEC02100 
SEC 02 110 
SEC02120 
SEC02130 
SEC 02140 
SECG? 150 
SEC 02 1 6 0 
SEC02L70 
SEC02180 
SEC0219C 
SEC 0220 C 
SEC 02210 
SEC 02220 
SEC 0223 C 
SEC0224C 
SEC0225C 
1SEC0226C 
SECO2270 
SEC 022 80 
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0117 

0118 
0X19 
0120 
0121 
0122 
0123 


0124 

0125 


0126 

0127 


0128 

0129 


0130 

0131 

0132 

0133 

0134 

0135 

0136 

0137 
013ft 
0139 


0X40 

0141 

0142 


0143 

0144 

0145 

0146 

0147 

0148 

0149 

0150 


0151 

0152 


0153 

0154 

0155 

0156 

0157 

0158 


REWIND SOATA 
PREFIXUI * 3 

400 RCADI SC AT A , 9390, END*9 100 ) CARO 
9390 FORMA T (20 A 4 > 

hKIIE (OUT PUT ) PREFIX. CARD 

IF < C ARD ( 1 ) .EQ. IEOS) GO TO 410 

GO TO 400 


PRINT STATS, CALC INVERSE COVAR HTX.SCAPRD, VECTOR PRO* OET 


C PHASE 1 ONLY 


REWIND CLASSX 
PREFI X (21 * 0 

ENTER LOOP TO READ AND CLASSIFY THE REGIONS. 

425 READ (CLASSX, ENO « 530) INFO 

426 PREFIX12) * PREFIXI2D1 


SEC02290 
SEC02300 
SEC 02 3 1 0 
SEC 02 32 C 
SEC 02 3 3 0 
SEC 02 34 0 
SEC02350 
SEC 0236 0 
SEC 02 3 7 0 
SEC 02 38 0 

410 CALL S EC PRT(ARRAY(C0VAR3)| ARRAY CAVAR3 ), ARRAY! SCR AR3), ARRAY (DETAR 3 ) SEC 02 390 
1 .ARRAY (C0NAR31 , ARRAY 1 0ATBS3D ' SEC02400 

420 IF (.NOT. PHASED GO TO 458 SEC024I0 

••♦♦♦♦♦•♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦•♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦SEC02420 

SEC0243 0 
SEC 0244 0 
SEC02450 
♦ SEC 02460 
SEC 0247 0 
SEC 024 80 
SEC02490 
SEC 02500 
SEC 025 1 0 
SEC 02 52 0 
SEC 02 53 0 
SEC 0254 0 
SEC 0255G 
SEC 02560 
SEC02570 
SEC 02 58 0 
SECO 2590 
SEC 0260 C 
SEC02610 
SEC 0262 0 
SEC 0263 0 
SEC02640 
SEC 02 65 0 
SEC 0266 0 
SEC 02670 
SEC 0 2o3 0 
SEC02690 
SEC02700 
SEC 02 7 1 0 
SEC0272 0 
SEC 02 73 0 
SEC02740 
SEC 02 750 
SEC02760 
SEC 02 77 0 
SEC 02 7 8 0 
SEC0279C 
SECO2Q0O 
SEC 02 S 1 0 
SEC 0262 0 
SEC 0283 C 
SEC 02 84 0 
SEC 02 850 
kSEC 02 860 


C GET RUN NUMBER 
C 

IF(CURRUN.EQ.INFOU) I GO TO 460 

CALL GETRUMINFO! D ,DATAPE , 10 .ERROR, RUNT A0, I MARK) 
IF(ERRCR.GT.O) CALL RUNE RR ( ERROR, INFO! 1)1 
430 DO 440 1 • 1 »NCFE7 3 

IFIFETVC3! I).LE.ID!5) )G0 TO 440 
VRI TE I PRNTR.9430 )FETVC3 ( I ) 

WRITE! TYPE WR, 9430 IFETVC3I I ) 

9430 FORMAT UOX, f CHANNEL « t i 3 > 

CALL ERPRNT U70,'STO P* 

440 CONTINUE 

C SET UP CALIBRATION SET VECTOR. 

C 

DO 450 J * 1,30 
DU 450 I - 1,3 

450 IF ( CSE T3( I , JJ -EQ. -50000.1 CSET3!!,j) * FRQCAL ( I *2, J I 
C CHECK FOR NEW CALIBRATION 
FLAG * 0 

CO 455 I - 1 , NOFE T3 


.GT. 1.02) FLAG » 1 


DO 455 J * 1,3 

455 IF< ABS(C5ET< J,FETVC3< I ) ) -CSET 3< J* FETVC3! It M 
IF< FL AG. EQ. 1 ) WR I TE ( T YPEWR , 94 55 > 

I F ( FLAG.EQ.l ) WR I TE I PRN TR , 945 5 ) 

9455 FORMAT!* 10214 CALIBRATION VALUES USED FOR CLASSIFICATION ARE NOT 
1 SAME AS THCSE USED IN THE STATISTICS') 

GO 10 460 


C SEC 02 8 7 0 

C PHASE 2 ONLY - CHECK ALL SUCESSIVE AREAS CN TAPE SECO208C 

C SEC02890 

C ♦♦♦♦♦*♦♦♦♦♦♦**♦♦*♦♦♦*♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦*♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ ♦•♦♦♦♦♦♦SEC0290Q 

C SEC02910 

458 READUNPUT) PR6F I X , PTS, L I NES, I NFO , 1 D, CSET3 SEC0292C 

IF !PTS .EC. 0) GO TO 426 SEC02930 

C SEC 02940 

C COMPUTE FIELD SIZE AND SET UP STORAGE ALLOCATION. SEC02950 

C SEC 0296 0 

460 LSAMNO “ 10(6) - 6 SEC02970 

IF (INF0(8) .LE. LSAMNO) GO TO 470 SEC0298Q 

WRITE tTYPEWR ,9460 ) INF018), LSAMNO SEC0299C 

9460 FORMAT I* 10212 LAST SAMPLE NUMBER <«,I4*M EXCEEDS LAST SAMPLE ONSEC03000 
1 TAPE. LAST SAMPLE NUMBER RESET TO*,!** 1 .’) SEC030IO 

WRITE IPRN TR ,9460 ) INFOC0), LSAMNO SEC03020 

_ SEC03030 

SEC 03040 


1C irKN IK 

INFO! 8) - LSAMNO 


PAGE 0004 
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SET UP JPTS AS NO OF CELLS HORIZONTAL LV TS AS NUMBER OF TRUNCATED * TS |f£g|3|g 


A70 JPTS * IfNFa<0J-INFOm*t«FOl9n/<INFCl*>*CELWTHl 
FT S =» JPTS * CELWTH 

RESET INFO O) TO CORRESPOND TO TRUNCATION 
I NF 0(81 « INFOITI ♦ TPTS-1 >*INF0I 9) 

llNFOt 5?°* 6 !nEO^aV W > H J*U I NFO« 5>~INF0< A> ♦INFO^ I ) /J I - INF0I61 

NStt * PTS+6 

CELSI2 - C6LWTN*CEL*TH 

NVR IS NUN OF SAMPLES IN ROATA TO SKIP TO GET NEXT CHAN OF DATA 
FOR THIS POINT 

NWQRO t^MAXOtNOPCOl, ICELS I Z+ll/2 J 
NN0R02 m N*0RQ*2 


ALLOCATE ARRAYS USED IN ACTUAL CELL SPLITTING BY ECK03 IfcoIl&O 

""«* «M Waf.T^jttKWfOK P1XCU F«0* HON“HGNOGEK ||S|||8 

CELIKE f DrStOR I NG LIVELIHOOD COOE FOR EACH CELL SECC33C0 

P?Jc§R F FCft R CORRELAUcfi tOATA VALUE PRODUCTS) OF Pl*||co|32C 
PiXVAL FOR ACTUAL PIXEL DATA VALUES SECC3330 

£lhi8S £g£ iVtm !K«1|S 

%isk$ m s^i s o? c s?£^ ELL mum 


SCC03-C7C 
SEC 03 0*f C 
SECC3CUC 
SEC03 100 
SEC 03 1 1 0 
SEC0112C 
SEC031 30 
StCOMAO 
SFC03 150 
SEC03160 

ift8M8 

SEC0M9C 

SECO32Q0 

utim 

5EC03Z3C 

SEC03240 

Hg8?!IS 

SEC 032 7 0 
SEC0>28a 

mum 

SEC03310 


- OAtBST 

« NCFET3P 
T. PHASE 1 
DATBS3 
OD A t A ♦ 

- ROATA + 

- PfXCOft 

- PIXVAL 

: ESttSIt 

* CELCLS 

: SfHSl 

- PRfSUF 


PHASE 2 ALLOCATION 


Celsuh 
VARS23 
IF I » NC 
BOAT* - 
ROATA - 

SUSS? 

CEL COR 


(N0F£T3*I 0!S) ♦ 71/a 

<NVR*CEt WTH*1 >/2 

♦ [ VARSZ3*CELSIZ+ I >/2 

♦ {NOFET3*CELSIZ+l >/2 


♦ I N*0RD*JPTS<*W/2 

♦ IPTSOl/* 


IF I. NCI* PHASES) GO TO *90 
STACK - DAT0S3 
CLOSED * SIACK * (JPTSM1/2 
OPEN - CLOSED ♦ IJPTS ♦ 31/A 

ftttU : KHxt l J TSK|m„ 

0ATBS3 - FLOllK ♦ f JP tI^* NCPD0L*l)/2 


*90 JF IDAT8S3 .GT, TOPI CALL ERPRNT t 139* *5T0P* I 
BUFRCZ » 50000/ t JPTS*2*CELSI 2 > . t 

IF iBuFROZ ♦Lf.r U CALL £RPRNTm9,'StaP*l 

i PRINT OUT PROCESSING PARAMETERS 

WlJTEfPRNTR *9500 1 CE LVTM *N0FE73 • NCPD0L*THRE5# CSELEC* BUFROZ 

9500 FORMAT I // / / ♦ OPRCCESS I NG PARAMETERS*// 

2 * NUHB6R ^^CHANnIlS^- «.|3// 
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SEC 1 NTS 


DATE » 77133 


21QQ9Q09 


FILE SECINT 


0196 

0197 

0198 


0199 


02 CO 
0201 
0202 


0203 

0204 


0205 


02 06 
0207 


0200 
02 C9 
0210 
0211 
0212 

0213 

0214 

0215 

0216 
0217 


0218 

0219 

0220 
0221 
0222 

0223 

0224 

0225 

0226 

0227 

0228 
02 29 
0230 


4 1 CELL SELECTION THRESHOLD * *,iP E12-4// 
5» NUMBER GF CELL LINES IN BUFFER * * f I71 


C PRINT RUN HEADER 
C 

CALL GETIME (TIME! 

WRIT€(PRNTR,HEAD) # * = 

CALL HEADER(N0FET3,FETVC3 f CSEL3,CSET3 f ID,FRQCAL> A SEC03680 

q *4 4*4444*44*4*444*4*4444*444444*444444 *4****** *4************* ********** SEC 0 3 890 


S6C038L0 
SEC03820 
SEC 0383 0 
SEC03840 
SEC03850 
SECOJBbO 
SECQ38T0 


C WRITE RECORD TYPE 
C 


SEC 03900 
SEC 039 1 0 
SECQ3920 


c *********** ****** ***************** 4*4* *********«*************4********4 SEC 0 393 0 


520 LINES =* UNFQ(5)-INFQ(4)>/INFCI<6> ♦ 1 


! iftES (BY • , 12, •)•/ 

(BY •»I2 r MVl 


CALL SECH01 (ARRAY (C0VAR3) #ARRAY( AVAR3 )« ARRAY ( SCRAR3 ) t 
L ARRAY (DE T AR 3 ) f ARRAY(C0NAR3)f 

C ARRAY (BDA TAJ* ARR A Y ( RQAT A I . AR RAY ( PIXCOR ) • ARRAY ( P IXVAL If 
L ARRAY 1CELC0R) 9 ARR AY I CELSUH) , ARR AY ( CEL IKE) • ARR AY ( CELCLS >, 
£ ARRAY ( PI XCLS \ , ARRAY ( 5 TACK) f ARRAYtOPEN)* ARRAY { CLOSEO ) 1 
£ ARRAY (FLDSIZ), ARRAY ( FLOCLS > , ARRAY ( FLDL IK) « ARRAY < PRTBUF ) » 
£ bufffru > ) 

IF ( .NOT.PHASEl. AND.PHASE2 I GO TO 458 
GO TO 425 

END OF FILE RECORD 


SEC 0394 0 
SEC 03950 
5EC0 39b 0 
SEC0397Q 
SEC039B0 
SEC 0399 0 
SEC 04 COC 
SEC04C10 
S6C0402Q 
SEC 0403 0 
SEC 04 04 0 
SEC 04050 
SEC 04 06 0 
SECQ4Q7Q 
SEC 04080 


C CORRECT AREA INFO PUT ON INTERMEDIATE TAPE TQ REFLECT CELL LINES 
C 

527 PRErixm - 5 

WRITE (OUTPUT) PREF I X ,PTS , L INES, 1 NFO , ID? CSET3 
PREFIX(l) * 6 

C PRINT HEADER. 

C 

WR I TE( PRNTR t 9520) ( I NFOU) 9 1*4,9) 

9520 FORMAT! 'OAREA PROCESSED LT 

C T2 1 1 ' COLUMNS * , 1 4 , •- * , I 5 1 f 

C**4* 4 4 * 4 44* 44*4 4 * 44 4 4 44****** ********* 4444 4**4 444**************4*****4*SEC0^090 
C SEC 04100 

C CLASSIFY THE REGION. SECOAUO 

C S6CQAU0 

C******4*4****** 4***4* *4**4 ******** **** 4444444444444 ******************** S £C 04 130 

SEC 04 14 0 
SEC 04 15 0 
SEC 04160 
SECQ4L70 
SEC 04 1 80 
SEC04190 
SEC 04200 
SEC04210 
SECQ4220 
SEC 04230 
SEC04240 
SEC0425 0 
SEC 04260 
SECQ4270 
SEC 042 00 
SEC04290 
SEC 04300 
SEC 04310 
SEC04 J20 
SEC0433Q 
SEC04340 
SEC 04 35 C 
SEC 04 360 


530 PTS « 0 

PREF1XU) * 8 
PPEFl XI 2) * 0 

WRITE ( OUT PUT , ERR* 536) PREFIX, PTS, LINES, INFO, ID, CSET3 
WRITE (PRNTR,9$35) 

9535 FORMAT COECF RECORD WRITTEN, M 
CO TO 537 

536 WR I TE( PRNTR, 95361 

9536 FORMAT COERRCR WHILE WRITING EOF RECORD. 1 I 

537 IFIC'JTPUT .EQ. CLASSR) RETURN 


C FINISH TAPE 
C 


( OUTPUT | 


CALL TOPEF 

540 CALL TCPRF (OUTPUT) 

10 * 0 
13 * 3 
K » - 1 

FILESV « F l LE S ¥♦ l 
PAEFlXil) * 1 
PREFIX(2> - 0 

WRITE (OUTPUT) PREFIX, RQTAPE, F1LESV, 13, K, (10,1*1,8) 
CALL TCPEF (OUTPUT) ^ 

CALL TOPRF (OUTPUT) 

CALL TOPEF (OUTPUT) 

CALL TCPBF (OUTPUT) 


SEC 04 380 
SEC 04 39 0 
SEC 04400 
SEC 044 1 0 
►5ECQ442Q 
SEC 0443 0 
SEC 0444 0 
SEC 0445 0 
SEC 04460 
SEC04470 
SEC 04430 
SEC 04490 
SEC04500 
SECQ451 0 
SEC 0452 0 
SEC 04530 
SEC04540 
SEC0455C 
SEC04560 
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FORTRAN IV G LEVEL 
FILE SECINT 


20. T 


SECINT L 


DATE - 77133 


21009009 


0231 

0232 

0233 
02 34 
0235 
02 36 
0237 
0228 
0239 
0260 
0261 
0262 
0263 
0266 


„ CALL TOPBF (OUTPUT) 

8900 RETURN 

9000 WR| TE ( TYPE WR *99000 ) 

WRITE ( PRNTR * 99COO) 

,90 °° GO R TO T 990O*** EN ° ° F R£C0R0 INTERMEDIATE TAPE -•) 

9100 WRI T£ l T YPEUR *99 100 ) 

WRITE ( PRNTR. 99100 ) 

’lloo WRire T i'Tvi?Sn o!^S? ECTED EN0 0F FILE RE * 0ING STATS FRO* DISK 


9900 WRITE ITVPEWft, 99900) 

NR I TE I PRN JR ,99900) 

99900 FO««*TI* FUNCTION TERHINATEO 

wALL K f rlAi H 

END 


(SUP I NT)* ) 


SEC06570 
SEC 0658 0 
SEC 0659C 
SEC0660C 
SEC 06 61 0 
SEC 0662 0 
SEC 06630 
SEC06660 
-•) SEC 06650 
SEC 06660 
SEC06670 
SEC0668 0 
SEC04690 
5EC06700 
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MTAPE mounts and positions the results or intermediate tape 
(or a tape to be used as output for copying results files) . 
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1 . Module Usage 
MTAPE 

CALL MTAPE (RQTAPE, RQFILE, MODE) 

Input Arguments : 

RQTAPE 1*4 Tape number of requested tape. 

A tape number of 0 is a request for 
a scratch tape. 

RQFILE 1*4 File number of requested file. 

If RQFILE is « 0 , then the tape will 
be initialized by writing a record 
type 1 on the results tape with file- 
type = 0. 

MODE 1*4 Flag indicating usage of MMTAPE . 

MODE = -1 indicates MMTAPE has been 
called to mount and position a tape 
to be used for copying results files 
onto. Mode = 0 indicates that a re- 
sults tape is being mounted for reading 
a results file. In this case, the 
tape is mounted ring out. Also, if 
MODE = 0, RQFILE = 0 is invalid and 
will cause an error when an attempt 
is made to write on the tape. MODE = 1 
indicates a tape is being mounted for 
writing a new results file (or continuing 
a suspended classification) . MODE = 2 
indicates a tape is to be used as an 
intermediate tape for reading an inter- 


mediate results file. The difference 
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between MODE = -1 and MODE = +1 is the 
DSRN used for the tape. For MODE = 

-1 and -2 DSRN is CPYOUT and for MODE 
* +1, DSRN is MAPTAP. {DSRN is MAPTAP 
for MODE = 0) . 

Output Arguments : 

RQTAPE 1*4 When MODE = 0, set to -1 if re- 

quested tape file was full and user 
decided to use disk for results . 
Otherwise, remains unchanged. 

RQFILE 1*4 When MODE = 1, set to -1 if re- 

quested tape file was full and user 
decided to use disk for results. 
Otherwise, sends back current file 
position of tape. 

MTAPE checks the validity of the tape by reading the 
record type 1 from the tape and verifying the tape and file 
number as well as checking fro the correct type of file. Any 
attempt to overwirte an existing file causes MTAPE to ask the 
user (via the typewriter) if he wishes to overwrite the file, 
respecify a new results card, or terminate the function. Note, 
however, that if a request has been made to initialize a tape, 
no checking is performed on previous contents. 

2, Internal Description 

See output description. 

3. Input Description 


The record type 1 of the results tape is read for each 
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file up to and including the file needed. That is, if file 
4 is requested the record type 1 is read from files 1-4. 

4. Output Description 

The following information messages are issued under 
the circumstances listed. The term filetype means the file- 
type code from record type 1 of a results file (the program 
uses variable CHECK for this number) . 

10041 is typed when a tape has been mounted and before 

MTAPE positions it. This message is not typed when 
the tape is being initialized or when the correct type 
number was already mounted. 

10043 is typed when MODE = +1 and file type of the requested 
file = 0. 

10044 is typed when MODE = +1 and filetype of the requested 
file = 1 and the restart flag from GLOCOM (RESTRT) is 
not = 1. 

10045 is typed when the tape is correctly positioned. This 
is not typed when initializing a tape. 

After 10043 and 10044, the user is asked whether he wishes to 
overwrite the file, respecify a new results card with a new 
tape and/or file or disk option, or terminate the function. 

10 100 is typed to allow entry of the new results card. This 
occurs when the user requests to respecify the results 
card. 

10101 is typed to confirm usage of disk for results and occurs 
whenever disk is specified on the results card. 

The following error messages are typed under the conditions 


listed. 
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E361 is written when the tape is being filed forward and a 
file is encountered with filetype other than zero be- 
fore the requested file is reached and MODE = 0. 

E362 is written when the circumstances for E361 occurs 

and MODE = 1. It is also written when MODE = 1 and 
the filetype of the file requested is ~ -1. 

E36 3 is written if the RESTRT flag is = 1 and the filetype 
of the requested file is not = 1. 

E364 is written when MODE = 1 and the filetype of the file 
requested = 1. 

E365 is written when an EOF is read on the results file. 

* 

This should never occur with valid results files. 

For message texts refer to the User's Manual. 

5 . Supplemental Information 

This section deals with the handling of tapes by MMTAPE. 
Input: 

If a tape is mounted on the device and it is the incorrect 
tape number (as noted from the appropriate status words in 
GLOCOM) , TOPRU is called to unload the tape before the correct 
tape is mounted. If the correct tape is mounted, MTAPE will 
check for the ring in if MODE = +1. If the ring is not in, 
the tape is unloaded and MOUNT is called to mount the tape with 
the ring in. If the correct tape is mounted, MTAPE assumes 
that the file number (as recorded in GLOCOM) is correct and 
moves the tape backwards or forwards to find the requested file. 
Output: 

The tape is mounted with ring in for MODE = +1 and with ring out 
for MODE = 0 and -2. 
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The tape is left positioned at the beginning of the requested 
file. When the tape is initialized a TOPRW is used to do 
this . 
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MOVE TAPE 
TO CORRECT 
. FILE 
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FORTRAN IV C LEVEL 

FILE NTAPE 


0001 

Gob? 


DATE > T7U7 


OdOL 

0005 

OOCi 


0014 

0015 



I NON /GLCCQK/ BUNK, CARDI201, CHKOUT, CQPFIL. CLASSR. CLASSX, 
CLUSIX, CCNPUT, CPYOur, CROROP. CRCSEO, DATAPE, 

OUPLTP, DUPRUN, ERRMSG, F0PNT, 

F 1 LES V , FLD9N0. HOATA, HEAO(86>, 10(2001 , I MAGEX , 


I NARK « XEY8D, MAPTAP, MAXCHA, MAXCLS, 
PAGSI2, PACK, POINT. PRESUX. PPhffl, READ IN , 


l 10), TENPASUO). 


RUt * 8 ARRAY 
REAL *4 FRQ£AU5,30l 

INTEGER • 4. CONENI ( 16 ) * OATEI 51 , HED1U61, HED2U61, TIM£(5J 
INTEGER • 2 BLANK2 
LOGICAL * 4 CHKOUT 
LOGICAL * 1 BLANK 1 

iovl VALENCE (DAT5AV. TOIL n, I CURRUN, I CI3T > , I FRCCAL (II, 10(51)), 

- 


i 


(SEPSCR,TP5TAT 12)1, I OUPI A .TPSTAT 131 I » 
(CCPSER,tP5TAT(MI . I TRAOL I , T PS T A 1 1 6 1 1 , 
( BLANK ,BLANK2 , CLANK 1 1 


INTEGER « 4 PREFIXI2I, RESC00I3I, CL1CCD ( 4 1 

DATA PREFIX/0,1/. RESULT/ ’RE$U , RESCOD/ ’CVER ■ , ■ RESP* , ‘TERN* /, 
l ClACOO/'f Apf 1 t’FlLE** 1 IN IT * » ’DISK’/ 


MMT0024G 

wm 

NNTOaZBO 
HHT0029C 
) ,MKIC0300 
MMTOOilC 
Mp»rco32q 
(0ASTAT,TPSTAr(4) ) .MMT0CI3J0 

58:88158 

NrtrOQ36C 

MMT00J70 

53:82158 

HHTG040C 
KHT004I0 


. ERROR 

. 11 
15 


» REPLY FROM USER AS TO WHETHER TO CONTINUE WITH RUN 
» FLAG NOTING TYPE OF FILE ON RESULTS TAPE 

0 - RESULTS FILE, -1 « CHECK HECCRO (WRITTEN AFTER LAST 
RESULTS OF RESTART FILE CN TAPEl, 1 - RESTART FILE 
* ERROR CODE FROM CPFUNC (NOT USED I 

CCNSTANT USED WRITING INITIAL RECORD ONTO TAPE 


MNl 

HHTI 


■MT0042C 


T00430 

T00440 


NODE 


NMTOQ4SO 
PNTC04SC 
MMTCC47C 
FHT0O480 
MMT00490 
PHTC050C 
PNTGOXIC 
NNT0052C 
MKT005JC 
MMTCQ54C 
MKT0055C 
MNT0056C 
MKT0057C 
PKIOOYaC 
MKTCQ59C 
HHTO0600 

_ 'MKT0061C 

TAPE HAD RESULTS AND DISK MKT0Q62C 

WILL SE USED FOR RESULTS, MMT0063C 

* FILE NUK8ER AS READ FRO** RESULTS TAPE MMT0064G 

- TAPE NUMBER AS READ FROH RESULTS TAPE MKTCC65C 

■ DSRN CF TAPE BEING MOUNTED KM-T0066C 

ThIS WILL BE MAPTAP UNLESS MCUNT ING COPYRESUITS CUTPUT TAP MMT0067C 

WHT0063D 

1 CONTINUE KMTOO&7U 

KMT007I0 

INITFG « O HMI0072D 

XMtOOTJO 

IF USER REJECTS T APE/ FILE PREVIOUSLY REOUESTEO. OISK WILL BE USED PHT0D740 

PM 100750 

)0Q |F (RQTAPg ,EQ -I .AND, RQFILE ,EQ. -I) RETLKN PMI00760 


» USED AS OUFPY READING VERSION NUMBER FROM TAPE 
■ INDICATES FUNCTION OF nhtape, supplied from caller 

WISH TO GET TAPE READY FOR WRITING CLASSIFICATION 
RESULIS, - 0 MEANS WISH TO GET T3PE READY FOR PRINTING 
RESULTS, » -1 MEANS MOUNT TAPE FCR COPYING CLASSIFICATION 
RESULTS ONTO. 

ROFILE > REQUESTED FILE NUMBER IN TAPE TO BE REACHED. IF RQFILE 
IS INPUT AS A 0, THEN MMTAPE IS TQ WRITE AN INITIAL 
RECUR C ONTO THE TAPE. 

CQNUNLE 


TFO 
(NO 
UNI I 


FORTRAN IV G LEVEL 
FILE MTAPE 


DATE » 77129 


0017 

0018 
0019 


0020 

0021 

0022 

0023 

0024 

0025 
oo?e> 
0027 

0026 
0029 


0031 

0032 

0033 
0054 

sill 

0036 

0039 


0043 

0044 

004-5 


0047 
004 B 
004 9 
0U50 

0051 

0052 

0053 

0054 
0053 


DOSS 

0059 

0060 
0061 


UNIT « MAPTAP _ _ _ 

IF (MODE .IE- — L 1 GO TO 108 
FILNO * FILESV 

IF CORRECT TAPE IS MOUNTED GO TO POSITIONING 

IF (MAPSAV 
105 IF (MAPSAV 
_ GO TO 110 
108 UNIT ■ CPYCUT 
FILNO - COPHL 
. IF ICOFSER -NE- 
no if (Mtinr — ■ 


JE .EO. 1 CAU MOUNT ( RQT APE, MAPI AP,*R| * ] 

,, fPOOf ,|y, 01 CALL MOUNT < ROT APE, MAPTAP, *R0* » 

S22I “ l) CALL NOUNT (RQT APE, CPYOUi , •Rt’ J 

IF (MODE .EQ, -2 I CALL MOUNT I RQT JP-EttPYOUT » *R0* 1 

IF CALLED FROM COPYRESUITS, THE UNIT IS COPVOUT 
FILNO * L 

RQFILE OF 0 MEANS THE CALL IS TG INITIALIZE THE TAPE 

IF (RQFILE ,NE, 01 GO TO 120 
11 « l 
U * 3 


gIET$jgi!Mlf F,x ’ R0, * PE - '*• ,3 '” 0F "- E - 1 


,91 


TAPE*, 15,’ TO FILE' 


call toprw 

,118 iMtlhi&r-pasif 

1 14,* (MIAPEI • } 

c CALL TfipRWIUNITl 
c READ FIRST RECORD AND SEE IF CORRECT TAPE 1$ MOUNTED 

CALL TlPRllwin 3C * Ef<0 * 5001 PftfFIX * ™°* * f0 * K * CHECK 
IF I TN0 EC. ROIAPE) GO TO 130 

Eitt " * mw MP£, * E ' ,Rt "" 

CD 10 110 

c POSITICN RESULTS tape to correct file 

c no IF (RQFILE .GT. 1) GO TO 140 

l 5 ? ll! 5 .aS c {SlN?l?i u m £ ,HcJl s ».i I «?i:c"S 5 iT‘ ,E ' ,CUM,ED 

?«iHi«i 5 0 ;uSjiJji i " OOE • E0 - 2 ’ 60 10 135 


if. ii. 


CALL TCPRUIUNITJ 
GO TO 110 
135 FILNO « 1 

CALL TCPRWIUNm 
CO TO 200 

140 J • RQFILE - FILNC 


J IS THE NUMBER OF FILES THE IAPE MUST BE MOVED 
200, 170 


THE DO ISO LCCP BACKS THE TAPE UP TO THE CORRECT FILE 


00 ICO I - 1,J 
FllNO « FILNO -I 


Jn e o?°cM2 ^29? FIVES THE TAPE AND CHECKS 

TD 0t SURE THAT FILES ARE REALLY ThERE TO SPACE OVER 


HHI00960 
MMT 00970 
MHT00980 
MMT 00990 
HRTOIOOO 


HKT0101C 

MMT01020 

MMTOIOJO 


>1320 
11330 
11340 
11350 
>1360 
11370 
>1380 
>1390 
I140C 
IU1C 
H42C 
>1430 
>144 C 


170 DO 180 I « 1, J 
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FORTRAN IV C LEVEL 20.7 


HEAD ILNIT,ERR-S30,END-500I PREFIX, TND, TFO, K, CHECK 
IF (CHECK -HE . 0] GO TO 210 
F UNO - FH.NO ♦ 1 
190 CALL TCPFFIUNU) 

C IF FUND AND HGFILE ARE UMECUAL , THE PROGRAM HAS COQDFED. 

£ 200 IF IF UNO NE- ROFILEI CALL ERPRNT 1 100, 'STOP* I 
C MOL READ THE FIRST RECORD OF IHf FILE AND CHECK VALlOtlV 

RIAL) llNlt,E»R'51l),ENn«SOO> PRU EX, I NO, TIO, K, CHECK 
C CHECK FOR ERRCK OF RESTART OPTION BUT NOT RESTART FILE 

C IF IR.ESTRT .EC I .AND. CHECK NC. I) CALL ERPRW7 1363, 'STOP' } 

IF CHOCE NE- O-AND KGDE NE. -2 1 (.0 TO 225 
IF (CHECK Mb -It GO 10 220 
I * FUND - I 

C THE FOLLOWING CCOE PROCESSES ERROR FOR INVALID FILE TYPE 

c 

210 CALL ERPRNT ( 36 1 M ABS I NODE ) ■ * GO* ) 

VRI IE (PRNTR .92 10> I 


IF A restart file. 


AND RESTART NCT SPECIFIED, ERRCR 


220 IF (CHECK EC It CALL ERPRNT (366, 'STOP') 

225 IF [CHECK EG -1 .ASD. RESTKT ,EC 01 GO TC 300 
IF (R.ESIRT .NE It CO TO 230 
IF (CHECK .EC It GO TO 3Q0 

C ERROR CF RESTARTING AND NO RESTART FILE 

C call ERPRNT (260, 'STOP*) 

230 (F (CHECK .NE 01 GO TO 233 
VRITfc ITYPEWR, 9230 J 
WRITE (PRNIft, 92301 

9230 FORMAT I* 10063 FILE HAS RESULTS IN IT - 00 YOU WISH TC *, 

i F i«f ’ -«Sl s f SS rI ¥1 . i«S1R2i e l '* 

23J J IF (CHE^K .NE^irii’io 235 
WRITE ITYPEUR, 92353 
WRITE (PRNTR, 9235) 

9235 FORMAT (• 10066 FILE IS A RESTART FILE - CO YOU WISH TC 

1 ■ OVERWRITE THE FILE, RESPECIFY THE ENTIRE RESULTS CARD, »,/ 

2 * Oft TERMINATE THE JOB (OVERWRITE, RESPECIFY, TERMINATE! *, 

1 * (MTAPEI'J 

235 COL »■ C 
ERRCOR - 3 

CALL CTLWRO (CARD, CCL, RESCCD, 3, CQCE* KEYBD, f RRC0R1 
IF (ERPCOR .CO 61 CALL RTMAIN 
CO I(J 1300,260.2661, CODE 

260 kRl IE ( 1YPEWR, 92501 
WRITE IKEYBO, 92501 

9250 FORMAT I* IOICO ENTER RESULTS CARD |MMT APE 1 1 1 

CUL * C 
ERRCOR * 3 

26 2 CALL CILVRC (CARO. COL* RESULT, I, COCE. XEYBO, CRHCORJ 
If (ERPCOA .EC. 6) CALL RTMAIN 

261 IF (CCL .EC 721 GO TO 100 

CALL CILPRM [CARD, COL, CLACOD, 6, CODE, (2651 
c CG ia (250 ,260,2 70*2301 , CODE 

C TERMINATE CPTICN 


COCE, XEYBO, CRRCORJ 


! RETYPE RESULTS CARD AFTER ERROR IN ENTRY 


GU TO 262 

C PROCESS PARAMETER ERRORS 

269 WRITE ITYPEUR, 9268) CARO 
9268 FORMAT ISX.20A6) 

CALL EftPRNI IERRMJM, «GO a l 
GO TO 262 

C OECOCE TAPE SPECIFICATION 
250 VECSZ * l _ 

CALL tVAL (CARO, COL, ROTAPE, VECS2, 62521 
IF IVECS2 .EO. 01 GO TO 252 
CO TO 2*3 
252 ERRNUK » 368 
CO TO 268 

C F TLc SPECIFICATION 
C 

260 VECSZ » 1 

CALL IVAL (CARO, COL. RQF1LE, VECSZ, 62521 
IF (VbCSZ EO 01 GO TO 252 
CO TO 263 
265 ERRNUM - 367 
COL * C 
ERRCOR ■ 3 
CO 10 262 


C DISK SPECIFICATION — WHEN DISK APPEARS ON CARD, IT OVERRIDES CTHER 
C PARAMETERS. CANNOT CHECK FCR WHEN BOTH SPECIFIED ON CARD SINCE 

C ROTAPE AND ROFUE ARE ALREADY USED UN MOUNTED TAPE I 

260 ||g*f£| m -1 

„ WRITE JTYPEWR, 9280! 

9260 FORMAT I* 10 10 1 DISK WILL BE USED FCR RESULTS (MHTAPEI ♦ ) 

GO TC 100 

300 WRITE (TYPEHR.9300) 

53CC FORMAT I » 10065 RESULTS TAPE MOUNTED AND POSITIONED*, 
i * (MMiApei*) 

CALL TCPRF (UNIT) 
ig ° Hoi?8 0£ iJ;C. -1» GO TO 600 
MAP5AV ■ RCTAPE 
FILESV - FILNO 
RETURN 

600 COPSE R » RCTAPE 
CUPF1L ■ FILNO 

return 

300 CALL ERPRNT I 365, 'STOP* 1 
510 REWIND KEYBD 
GO 10 280 

530 CALL ERPRNT I?6l,«ST0P') 

END 


LARS Program Abstract 

MODULE IDENTIFICATION 

Module Name: SECHOl Function Name: SECSUP 

Purpose : Performs Supervised ECHO Processing on area. 

System/Language : 

Author : P. D» Alenduff D ate : 4/11/77 

Latest Revisor: Dates 



MODULE ABSTRACT 

This routine is called once for each area to be classified 
and performs all classification on each line of data and 
produces the output file. 


PURDUE UNIVERSITY 

Laboratory for Applications of Remote Sensing 
1220 Potter Drive 
West Lafayette, Indiana 47906 
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1. Module Usage 


SECHOl 

CALL SECHOl (ICLSCV, VECPRD, SCAPRD, DETCOV, CONST, BDATA, 
RDATA, PIXCOR, PIXVAL, CELCOR, CELSUM, CELIKE, 
CELCLS, PIXCLS , STACK, OPEN, CLOSED, FLDSIZ, 
FLDLIK , PRTBUF, BUFFER) 


Input Arguments : 
ICLSCV 

(VARSZ3, NOPOOL) 


VECPRD 

(NOFET3 , NOPOOL) 


SCAPRD {NOPOOL) 


DETCOV { NOPOOL) 


CONST (NOPOOL) 


R*4 An array dimensioned {channels 

* (channels + l)/2) * number of 
classes, containing the upper half 

of the triangular covariance matrices 
for each of the classes used for 
classification . 

R*4 An array dimensioned channels 

* classes containing the vector pro- 
duct of the classes for each channel 
from the inverse covariance matrix 
and the determinant of the covariance 
matrix. 

R*4 An array dimensioned "classes" for 
the scalar product of the mean and 
covariance matrices for each class. 

R*4 An array dimensioned "classes" 
words for the determinant of the co- 
variance matrix for each class. 

R*4 An array dimensioned "classes" 
words containing the constant term used 
consisting of an expression with 
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BDATA 

( (N0FET3*ID { 6) )/2) 


RDATA 

(NVR, CELWTH) 


PIXCOR 

(VARS 2 3 * CELSIZ) 


PIXVAL 

(NOFET3 * CELSIZ) 


CELCOR (VARS Z 3) 


CELS DM (NOFET3) 


CELIKE 

(NWORD, JPTS) 


determinant and scalar product for 
each class. 

1*2 The array used for reading the 
byte information of data values from 
data tape. 

R*4 The array dimensioned "(ID(6)-6) 

* cell width" words for reading actual 
data values from tape and calibrating 
them. 

R*4 The array dimensioned (channels * 

2 

(channels + l)/2) by (cell width) 
words for storing the cross product 
terms of the data values of a cell 
(correlation term) . 

R*4 An array dimensioned "channels * 

2 

(cell width) " for storing actual 
pixel data values for each channel for 
1 cell. 

R*4 An array dimensioned (channels * 
(channels + l))/2 words for saving 
the cross product of the pixel values. 
R*4 An array dimensioned "channels" 
words used to store the sum of data 
values for each pixel in the cell for 
a given channel number. 

R*4 An array dimensioned "NWORD by 
columns of cells" for storing the cell 
likelihoods or pixel values for each 
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CELCLS ( JPTS) 


PIXCLS 

(NWORD, JPTS) 


STACK (JPTS) 


OPEN (JPTS) 


CLOSED (JPTS) 


FLDSIZ (JPTS) 


FLDCLS (JPTS) 


cell in the line. 

1*2 An array dimensioned "column" of 
cells halfwords used to contain the 
most likely class number of all cells 
in a given row. 

1*2 The array dimensioned NWORD 2 by 
"column" halfwords for storing the 
class number of each pixel in a given 
cell for singular cells (equivalent 
to CELIKE) . 

1*2 A stack of field numbers for the 
cell in the given position for a row 
of cells. 

L*1 An array of "column" bytes used 
as logical flags to store whether or 
not a particular column of cells has 
been opened as a field. 

L*1 An array of logical values for each 
column of cells to indicate whether 
the field for each column has been 
opened for annexation. 

1*4 The array of "column" halfwords to 
store the length of the current field. 
1*2 The array of "column" halfwords to 
store the length of the current fields 
as the number of previously annexed 


rows in this column 
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FLDLIK q R*4 The array of likelihoods dimen- 

( NOPOOL, JPTS) 

sioned columns of cells by classes 
words used to store the likelihood 
of each field for each class. 

PRTBUF (JPTS) 1*2 The output print buffer dimensioned 

"columns of cells" halfwords. 

BUFFER 1*2 The buffer array used for annexa- 

(BUFROZ * JPTS) 

tion of neighboring rows of cells 
dimensioned BUFROZ * JPTS. 

SECH01 performs processing of an area to be classified. 

If only an intermediate tape is to be produced, cell information 
is gathered and written to tape. Otherwise all information is 
used in annexation and classification, and the classification 
results are written to the output device. Any maps are also 
produced . 

2. Internal Description 

Initially, column headers are printed if any map is to be 
produced. All buffer pointers and cell arrays are initialized 
so that the arrays OPEN and CLOSE indicate that no annexation 
and classification has taken place. Then a loop is entered and 
executed once for each row of cells (two or more data lines) . 

The routine FILBUF is called to gather statistics on one row 
of cells, either directly from processing the data tape, or from 
reading the intermediate tape. If no annexation is to be per- 
formed, the information is written to the output device and this 
loop is executed until this entire area is completed and a return 
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is executed. Otherwise, a series of calls to LIKRAT for a com- 
parison of the likelihood of the cell being annexed to the cell 
above, left, and right is performed. The cell is annexed to 
the most likely cell if the comparison with the annexation thres 
hold is successful. This is continued for all cells in the row. 
When finished, any fields formed from adjoining cells that are 
closed and no longer continued are classified and then any new 
fields are opened in that statistics from the initial cell are 
copied into the array storing field statistics. Then LIKRAT is 
called to classify closed areas. A print buffer is then filled 
with symbols for either the cell map or classification amp and 
printed out. Statistics are handled for remaining open fields 
and the row is written to the output tape. This is continued 
until the area has been processed. 

3. Input Description 

Not applicable 

4 . Output Description 

The unit PRNTR receives the map and column header produced 
if a map is to be printed. 

The unit OUTPUT is used as the output unit for either inter 
mediate tape - cell information or classified cell lines. These 
lines follow the format of regular classification records. 

5 . Supplemental Information 


This program uses common blocks GLOCOM and SECCOM. 
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FLOWCHART 
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20.7 


5CCH01 


OAF E * 77129 


15021059 


PACE 0001 


FORTRAN IV C LEVEL 


FILE SECH01 


0002 

0003 


S6CH01 LARS XXXX 

SPECIAL VERSION 10 IMPLEMENT EXTRA SCANNING LOGIC. 

subroutine seckoi iiclscv, vecpud, sc apra.ostcov, const. 

1 BOA TA ,RDA TA , P I XCQK , PI XV*L « CElCPR tCELSUM. „ 

2 CCUX£tCELCt.S.P|XCLS r STACA l OPENTCLaSEOtFLDSl2,FLOCLS. 
i FlDLIMPRIBUF, BUFFER* .. 

IMPLICIT I Mt-GtR * A IA-2I 

COMHQN /GLCCOM/ BLANK., CARO (20 1 < CHXOUT, COPFIL, Cl. AS SR, CLASS* , 

1 ClUSTX, CONPUlt CPVOUT, CHORD ft, CRDSEC. DaTAPE, 

2 0UPL7P, 0UPRU6, ERRMSC, F SPfJT , 

J FILESV, FLO0ND, HDATA, HEADIBS), 10(2001, INACEX. 

4 (MARK, KEYBD, MAPTAP, MAXCHA. MAXCLS, 


I r AKK> KtTUU, RAP [Art RAXCHA, R AALL5 , 
PACSli, PNCH, POINT, PRESUX. PRNTR, READ IN 
REJTRT, RUNFUt RUNT AB I 10,31. 

SOATA. SEPARX, SEPTPX, *»*• 



She! 101, TEMPASI301, 


OATEISlt KE01I16I, 8ED2U6*, TIHEI5J 


(OAT SAV, 10(11), I CUftRUN* 1 DI3I I » I FRl 

midffiglttii.'HWUiRSWfeil!* 

iSI?1K:jpH? i : 

ICCPSER. TPSTA7 (SI I , t TRACUF ,TPST AT I* 

(BLANK, BLANX2,0LANK1) 


0012 
OC 13 
001A 


2 NO I 

J PI 5 . 

J FEIVECI30), rt 
*OFLAG, CLSMAP 


;ls, ndfei, 

IL,riHLKU,NKU^U2 f f*VK, LiU 1 FU1 , PUIM AM l 2, 6U t , PREFIXI2I, 
ROFILE, SVMCNf, VARS73. CSELI30I, CS£l 3 ! 30), 
FETVC3I30I, PCLPTRC2.60I, P0L$IK(60), 

llsmapj csEn(3,30), oeJMAp, phase i , phase;, pclnm1(601. 


, prstat! 

REAL * A CSELEC, CSET, CSET3, THRES, ANN l 

INTEGER * 2 CSEl, CSEL3 , FETVEC, FETVC3, PCLPTR, PDtSTK 

LOGICAL • I COFLAG, CL5MAP, CSEU, OBJMAP, PFASE1, PHASE2, PQLNH1 , 


ECHOCAAC 

ECHooAso 

ICH00^60 


0016 

0017 


0019 

0020 
0021 


c define LOCAL VARIABLES 


real * A RC AT A I H VR, CE LWTH ) ,CELlK£jNMORD*JPTSl,FLDLlKl 6CPC0L, JPI SI, 
I IClSC V I VARS 7 3 rNOPOOL I , V£ C PRO (NQF ET3 , NOPCCL I , SC A*RO I NCPCOL I , 

20E ICCV INOPCCLI .CCNSTINOPOCL) . CELSUMC NCFETJ 1 rCELCORT VARS23) , 

1 PIXVALlll , P (XCOFM 1 ) , AlJXL l X 16 O.A ) , SUP .SR l GPT , SLEPT , SUPLUS, ANNt 

INTEGER * A FLDSIZ1JPTSJ 

INTEGER • 2 CELCLSIJPTS) .PUCLS(NtiOR02,JPrS)«eCATA( 11,S7ACX(JPrSl, 

1 BUFFER!!!, FLDCLSt JPTSl * PRT0UF I PIS P .BLOCK! A) , M2, ClaSSZ, 

2 JPTS2 * SIKGUL, VERTI, HORII, PLUS, L2, BLANK?, AUXCLS ( A 1 , KSYH 1 1 1 

LOGICAL » l OPEN! JPTSl * CLOSED (JPTSl, STOP 

FOUI VALENCE |FSTLIN,|NFOI Al J, I LSIL IN, INFOI SI I , I LIMINT, INF0I61 1 

EQUIVALENCE ICURLIN.INFOIAI), (SYHtU.hSYMI 1)1 

DATA SINGUL. VERTI, H0RI2 , PLUS, L2/ '0» , * !• ,'L»/ 


FORTRAN IV G LEVEL 20.7 
FILE SECKOI 


DATE * 77129 


0026 

0027 

0020 

0029 

0030 

0031 

0032 

0033 
OOJA 

00 35 

0036 

0037 
00 JS 
00 39 
OOAO 
00 Al 

8855 

OOA* 

00A5 

OOA6 

885? 

00A9 

0050 

0051 

0052 


005A 

0055 

Soil 

005B 

0059 

0060 
0061 
0062 
0063 
006A 

0065 

0066 
0067 


0060 
00 69 

0070 

0071 


RQWSI2 * JPTS * CELS1 Z 

TURN OFF COLUMN HEADER IF NO PHASE2 PROCESSING 
IF I.NCT. PHASE 2 l 00 TO 90 

transform annexation parameter 

ANNl * -ALGGf 10. I*?HRES 


CBJMAP) CO TO 90 


C PRINT COLUMN NUMBERS ON OUTPUTV 

N - |NF0(9I«CFLMTH 

IF ( NOT CLSMAP AND. .NOT. 

L - INFO f 7 I 

INTVL * (JPTS * 110) f 111 

JPTSl « 

P - L A 
11 * 1 
IF (M .GT 
1 - L 


11 - 2*JPTS1*1 

nn 40 I«l r ii * 

I+JPti 

yr | pou i 

90 au format i* 'jsx.inm 

JMK .l£. 9991 GO TO SO 

D0"a 5 J* 1 , JPTSl 
PRTBUF(J) » KO01 1,101 
AS I * 1 ♦ N 

MR I TE ( PRNT R , 90A 0 ) 1PRTBUH 1 1, J*l. JPTS II 
50 VRITE (PRNTR, 91201 
9120 FORMAT COM 


£ 


90 JPTS2 • JPTS 

£ SKIP ZEROING OLT ARRAYS IF ONLY PHASE 1 C THEY'RE NOT ALLOC ATEO — 
100 IF ( NCT ^PHASEZI GO TO 105 


I IF ( NCT PHA 

STACK? II - I* 

CLOSEOUI - 

110 GPFN(l) » .FALSE. 

10S httli tlNINCCELMTH 
NOLlNt « 0 
PUNE * FSTLIN 
NEXT * t 
PR TRW - 1 
OLORV * O 
NEWRW * 1 
CURL IN • FSTLIN 


; ENTER MAIN LGOP 


ISO J - NOLINE*CELWTK 

IF IMOO (NOLINE ,501 .EO. 0) WRIT6ITYPEWR.9 ISO J 
9150 FORMAT 116, * LINES CLASSIFIEO.M 
NOLINE • NOUNE ♦ 1 
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FORTRAN IV C LEVEL 
FILE SECMOl 


DUE * 77129 


00/* 

outs 

00)6 

0077 

ou/l 

oS®? 

008) 

00»4 

0085 

0086 
008 7 
Q068 


00 9 3 
0014 
0095 

8S3? 


otoo 

OfCl 

0102 

0103 


oio* 

0105 

0106 


8SH 

0112 

0U3 


i ill INI Rrik nr HAt A 


CO TC 156 


155 IF ( NO I CBJnJP) GO TO 157 

&Q J 1jZ l»l , JPTS.IMVL 
PRT8UF1IIJJI i finNX? 

IF ICELCLSCII ,NE.Ol I 

156 F 35 cS« i {^ - *' 1 

157 

NRl lEICUfPLTl PREFIX, CURLIN, 

ff 


. » I A , LUK 

iRUN ■ CLRLIN ♦ GLCR 
- 0 GE- 


II, I-1,JPTSI> 
CURLlNfCELCLSrCEUKE 
RETURN 


C INSHI IN BOEFEft CL^SS NUMBERS FOR EACH PIXEL OF SINGULAR CEltS 

C HC J * NRCV5S 

DO 200 C6LNUK»1,JPTS 
BUFFER (JJ * 0 

IF 1CELCL»ICELNU1*7 .Wfc. 0) GO TO 200 


JP7S2 


200 J p j ♦ C ELS 1 1 


228 N * I 
H * l 

13 * QRDhBS 
JJ - NROwBS 


00 2«0 CELNUH rn 1 T JPTS 

12 « 15 

13 * 13 ♦ CELS I i 
CELSIZ 

r SUFFER! 12) 

IFIFLOUP .LE. JprSl CLOSEOIFLOUP) - .FALSE. I 

IF (SUFFER! J2> .EQ 01 GO TO 229 
N ■ «l 
GO TO 240 

229 SUP « - IE* 70 
SRI GH T m * 1E+70 
SLEPT - --1E*70 

I TRY TO COMPARE CELL WITH FIELD ABOVE. 

.CR FIDUP.GT JPTS) GO TO 230 

C N .NE. M IFF CELL 'CELNUN* HAS ALREACY BEEN COMPARED TO FIELD *FLOUP* 


Echo: 

echo: 

ISE! 

ECHO, 


ISS 


2H 


FCH02100 

ECH021IO 

ECHQ2120 

FCH02130 


FORTRAN IV G 
FILE S6CH01 


0119 

0120 


LEVEL 20 T 


DATE - 77129 


0121 
0122 
□ 123 
0124 


0125 
□ 126 

0127 

0128 

0129 

0130 

0131 


0136 

01J7 

0138 

0139 

0140 

0141 
D 142 
0143 


IN WHICH CASE SETTING K«N IS ALL THAT IS NECESSARY. 

IFfN .EQ Ml CALL LlXRATI SUP.CEllKEl 1 .CElMJH l.f LOL IK( I, PLDUPJ . 
C AUXLlXU.M>,CELCLS(CELhUN|,FlDCLSIFL0UPI.AUXCL5(M)i 

230 M . u 

TRY TO COMPARE CELL WITH FIELD ON LEFT. 

LFT - 0 


IF?£eLNUM .EQ. 1) CO TO 232 
FLDLFT * bli-ferijd 

lPiPLDLFT.Le.JPTS AND. FLCLFT.GT.O . A6C FLOLFT NE.Fl 

flgsiziflolf 

. CALL UKRAT1SLEFT,CELIk£U,CELN0P),FLCLIX11.ELCL?T1, 

; CELCLSiCELNUMI v FlOCLSIFLOLFT),AUXCLS 121 1 


tCH0229C 

ECH02300 

ECH02310 

ECH0232C 

ECH02330 ' 

ECH02J4C 

ECH023S0 

ECHD2360 

ECH02370 

ECH02390 

ECHD231C 


TRY TO COMPARE cell WITH FIELD ON RIGHT. 

232 IFICELNUM .EQ JPfSI GO TO 234 
FLORT * BUFFER! 13 ) 

IFIFL0PT.G7.JPTS .OR FLOSIMFlCRTJ GE 


.OR. FLORT CO FLOLFT) 


..fldup And. 

ILFTI.LT.KAXSI2)ECH0241C 
■ AUXLIK(1,2I* ECH0242C 
ECH02430 

li 8888 

ECH024feC 
ECH0247Q 
ECH0248C 

'^FLORrtEo’FlDUPl CD TO 2346CH025CC 
T) GO TO 234 LCH02510 


IF ( H .EQ. 4) N*l . ECHi 

CALL LlXRAI(SUPLUS.CCUKeUiCELNU^«l I .FLOLlKI 1* FLOAT I , AUXl IK l 1 , Hi .ECU 
C CELCL SICELNUMM I ,PLDCLS( FlDRTI , ALXCL5INI ) ECHI 

IFISUPIUS GE ANN1I CALL LIKRAl(SRlGFT*CEL|RE( l.CELNLK), |CH| 

t AUXLIM1*NI.AUXHRI1 ,31 ,CELCLSICCLNUP1, AUXCLS(N) t 6UXCLSnn ECh 


: FIND GREATEST LIKELIHOOD RATIO. 

234 TF (SUP LT SLEPT .OR SUP. LT.SRIGHT1 GO TO 236 
IFCSUP .LT. ANN 1 1 GO TC 240 
OPEN! FLDUP 1 * .TRUE. 

FIELD - FLDUP 
X m M 
GO TO 239 


236 


s ;®W 2 ia 8 

mrFKtii • m,t - 

K » 3 


tasstss 


0144 

0145 

0146 

81« 


KSRMbti" 

GO TO 239 


23 * ^^FlHfT^" C ° T ° ”° 

K » 2 

. ADD CURRENT CELL TO FlElO 

239 CALL E ACuilFLOLIMl .FIELD) , AUXLIM 1 ,K) .FLOCLSIF JELO) * AUXCLS IK), 
6 FLOSWiFIELOM 

240 Jt ■ J 2 


; close fields That are in gldrm but not in KEWRW 


0153 CO 260 FIEL0-1, JPTS 

0154 IFICLOSEOlFtELDI .CR. 

0155 hEXI > NEXT - 1 

01*6 STACK^EKT) - FtELO 

C CLASSIFY CLOSED FIELDS 
0157 c O.ASS2 - FIO:ISIFIELO) 


OPENIFI ELOJ I CO TO 250 
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FORTRAN IV G LEVEL 20 7 
FILE SECH01 


SEC HO 1C 


date » 77129 


150210S9 


0 1 SB 
0159 
0140 


0161 

0164 

0165 

Silt 

0168 

0169 


iu utiunn 

GO TO iOOO 

IIS S®1SiI{Sl°S Ui'dt- 

B AC F SCAN ACRUSS ROW. BEGIN NEW fIFLDS. 

’ 265 * • NRCWBS ♦ ROWS I Z"C f l S W 
I • JPTS 
DO 2dO N*l ,JPIS 
L * K - CE LSI l 


I =• I-t 
H - BUFJ 


FERIXJ 


Sill 


SiS? 

01 W2 
0163 
0184 
0165 


01 £9 

0190 

0191 

0192 

0193 

0194 

0195 

0196 

0197 

01 SB 
0159 

02 CO 
0201 
0202 
0203 


C IF H GT. JPTS THIS CELL IS SINGULAR 
IF f H .GT. JPTS) GO TO 260 

C PLACE CLASS Of FIELD H TWO EACH PIXEL OF FIELO H 
C It 1 M . NE OJ GO fo 270 

« * STAC*! NEXT) 

K£k* » NEAI ♦ I 
BUFFER (RJ » H 

C CALL BEGFLO TO INITIALIZE STATISTICS FDR NEW FIELD 

CALL BEGFLC(CELCLSIN) .FLDCLSf H) tCELIKEI 1.N).FL0L1R| 1,H) .ELDSIZ (HI ) 
270 IF1I-EC.0 DR- BUFFER CL) NE.O> GO TO 280 

C hOOl F 1 E 0 

CALL LlXRAI(SRtGHTrCELlKEfI«t)tFLCLTRll*H) ,AUXL IR»CELCLSM I « 
t FLDCLSCH tAUXCLS) 

IFISRIGHT LI AHNI) GO TO 280 

Call adozifldlimi.hj .auxux.floclsihj.auxcls.fldsizi fh 


275 BUFFER CD * H 
280 K * L 

1FC0LDRW .60. 0) CO TO 340 
1290 FORMAT \ 19< 1X.L11A1 ) 

340 CLRLIN » CLRLIN » CLCR 

IFlCURLIN GT. LSTLlNl GO TO 400 


C POINT TO NEXT ROW IN BUFFER 
C 

QLORW * NEWRW 

360 NEWRW * NOOlNEWRW.BDFROn » 1 
IFINEWRH .HE. PRTRWJ GO TO ISO 

C CLASSIFY AM) PRINT ONE ROW OF CELLS 

K » 1 ♦ iprtrw-imrowsu 

L » K ♦ R 04 S 1 7 -CELS 1 Z 
00 365 I »K [L.CELSIZ 
H - BUFFER*! > 

IFIM GT. JPTS) GO TO 365 


(FIBUFFERCJ) .ne. 

P ■ J ♦ CELS 2 T 
CO 361 N • J.N 
361 BUFFER (Nl 


Hi GO TO 364 


y- ' £lasS2 

364 CCNTINL6 

365 CONTINLE 

ASSIGN ISO TO RETURN 
GO TO 4000 

C PRINT LAST PRINTER LINE. 


PAGE 0005 


FORTRAN IV G LEVEL 20 7 


DATE • 77129 


FILE sechoi 


0206 
0207 
02 C« 

0209 

0210 
0211 
0212 

0213 

0214 


0215 

0216 
0217 
0210 

0219 

0220 
0221 


C CLASSIFY all OPEN FIELDS 

L - NRCMBS * ROMS12-CELS12 
00 420 I ■NRDWBS.L.CELSI 2 
FIELD - BUFFER ( I ) 


ASSIGN 420 TO RETURN 
GO TO 3000 
420 CONTINLE 

% FINISH PRINTING BUFFER. 

ASSIGN 430 TO RETURN 

430 IMPRTRW .ED* NEVRWI ASSIGN 500 TC RETURN 
GO TO 400Q 
SCO PL INI « O 

PrEFIX(I) » 7 

WRlTE(CUTPLftER**600> PREFIX, PL INE.PRTBUF 
RETURN 


C END OF TAPE OR JO ERROR 


0 FQRNAT 1 *^| C^47 END DF TAPE OR AN 120 ERROR QCCUftED. ». 
t 'CLASSIFICATION TERKINATEO AT LINE* » 16. « f ECH0A5) * ) 

kRlTE UYPEWR.9610) PUNE 

«irBw*i«b»m«. ,B 500 

i « i 

CALL TCPSSlOUTPUT.I. ERROR) 

IF 1 I .LF. Cl RETURN 
BACKSPACE CUTPUT 
GO TO 500 


: PROPAGATE CLASSIFICATION THROUGHOUT BUFFER. 


“II* 

# 

0242 

0243 


0246 
0249 
02 50 
Ok 51 

0252 

0253 

0254 

0255 

0256 


.TRUE 

C PROPAGATE CLASSIFICATION THROUGHOUT RCW. 


LD) €0 TO 3002 


CD . 

IFIL 

STOP « - 

U > u ♦ cELsn 
DC 3001 15*13,14 

3001 BUFFER M5) m CLASS? 

3002 row.EO.PRTRW) GO TO RETURN, 1250,420) 
ROM » ROW - 1 

GO TO 3000 

C PRINT ONE ROW OF CELLS 

4000 II * 1 ♦ IPRTRW-1 I4R0WSIZ 
12 * II *• CELS ! I - CELWTH 
00 4002 13*1 l , iZtCELWTH 
15 * 0 

14 « 13 ♦ RUMSI2-CELSIZ 
DO 4001 16*13, 14, CELS12 
17 ■ 16 ♦ CELWTH - 1 
DO 4001 I6*|6, I 7 

“ .15* i 
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FORTRAN IV G LEVEL 20 1 


DATE • 7fL29 


IF KNOT. CLSKAP) CO TO A 002 

K I Intvlfcelwth 

IHW0D{L1N»kT -he. 01 GO TO 4002 

PRTBUMIT " hSyK^RIBUF I Jl*2-1> 

IF (t.GT.JPTSI) PRT8UFII1 » BLANK? 
I - I ♦ 1 


A0D3 v§JT|iF*N T R*9AOO*»PHNE »(PRTBUF(I J»I»l*JPT$l) 

’jbi man wmi ih» #l 

F* TRW ■ MOOlPRTRR.BlJFRO/l 
CO TO RETURN* 1140, 430, 40i 


END 


0) 




MODULE IDENTIFICATION 


T 


LARS Program Abstract 

ECHO FUNCTIONAL 

Modulo Name : GATHER » Function Name: SUPPORT 

Purposes Gather cell statistics information for processing. 

System/Language : CMS/fortran 

Author i P. D. Alenduff pate : 4/1 3/77 

Latest Revisor: Date : - 


MODULE ABSTRACT 

GATHER computes information about cell statistics for actual 
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1. Module Usage 


GATHER 

CALL GATHER ( RDATA , PIXVAL, PIXCOR, CELSUM, CELCOR, CELWTH , 
VECSIZ, MTXSIZ , NSR, NVR) 


Input Arguments ; 
RDATA (NSR, NVR) 

CELWTH 

VECSIZ 

MTXSIZ 


NSR 


R*4 Calibrated data values from data 
tape with NSR samples and NVR channels. 
1*4 Number of data points in width of 
cell. 

1*4 Number of channels to be used 
in classification. 

1*4 Number of entries in the tri- 
angular covariance matrix. Equal to 
(channels * (channels + l))/2. 

1*4 Number of data samples per line 
in RDATA. 


NVR 


Output Arguments ; 

PIXVAL (N0FET3 , 
CELWTH* CELWTH) 


PIXCOR (MTXSIZ, 
CELWTH* CELWTH) 


1*4 Number of channels of data in 
RDATA array. 

R*4 An array of actual data values in 
the cell dimensioned channels by (cell 
width) 2 . 

R*4 Cross product terms for all data 
points in a cell, dimensioned (cell 
width) 2 by MTXSIZ. 

R*4 Sum of data values in the cell for 


CELSUM (VECSIZ) 


each channel of data (dimensioned 
channels) . 
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CELCOR (MTXSIZ) R*4 Correlation matrix for data 

values in cell consisting of a sum 
of products of data values in cell* 
GATHER is called to return information concerning the 
data characteristics of one cell of data. 

2. Internal Description 

GATHER loops for each channel of data and each pixel with- 
in the cell to produce the data array returned. 

3. Input Description 

Not applicable 

4 . Output Description 

Not applicable 

5. Supplemental Information 


Not applicable 
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LARS XXXX 
GATHER STATISTICS 


ON A GIVEN CELL 


0001 
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0012 

0013 
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0018 
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0023 
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0031 
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00 33 

0034 

0035 

0036 


GATHER lmi\j aaaa GAT00010 

GAT0C020 

ST«iSTi8r«i D J» / H^r...5 S t-KSm6™ gatoooso 


C=GATHER 
C COMPUTE 


yA£r,5J.?L ISTICS F0R THE CELL WHOSE FIRST DATA VALUE IS RDATA(l). GAT00060 
F S MT5§y7 I NlQ G SIoi R(RDATA,PIXVAL,PIXCOR,CELSUM * CELCOR » CELWTH » V6 CSIZ,GAT00070 

IMPLIClf N INtEGER r.^rn^r 


c 

c 


c 

c 

c 


c 

c 


c 

c 

c 


DO A 10 4 L = l*vicSli PIX ^ A ^^' PIXCOR< l) » CEl - SUM(VECSIZ >*CELCGR( MTXSI Z> 

10 CELSUMIL>= 0. 

DO 11 L=1 t MTXSI Z 

11 CELCOR(L) = 0. 

N = C£LWTH*MTXSIZ 
Jl = 0 
J2 = 0 
L2 = 0 

DO 4 J=1 tCELWTH 
K1 = JI 
K2 = J 2 
= (J-1)*NVR 
= (J-L)*N 
DO 3 K=1 »CELWTH 
K1 = K1 + 1 
M2 = 0 
LI = K 1 
K3 = L2 
= Jl + K 

- J2 + ( K-l ) ♦MTXSI Z 

■oA K i i tam i ,J ' U * vEcs 1 ,z * cel "th 

Ml = K3 
L2 = L2 + 1 

= K ♦ IL-1)*NSR + IJ-1)*NVR 

%IXVAial!*^§il*(L{f _l '* V6CS,Z * CELMT H 
CELSUMIU = CELSUMU) 4* PIXVAL(L2) 

DO 1 M -1 • L 
Ml = Ml + l 
M2 * M2 + 1 
M3 * K2 + M2 

M2 = M ♦ L*7l-1)/2 SIZ + < J-1>*VECS1Z*CELWTH 

M3 = M 4- L*(L-l)/2 + (K— 11*MTXSIZ ♦ <J-l)*N 
PIXC0R(M3) = PIXVAUM1)*PIXVALCL2) 

CELCOR (M2 1 = CELC0RIM2) 

LI = LI 4 NSR 
K2 * K2 + MTXSI Z 
Jl = Jl 4- NVR 
J2 * J2 + N 
RETURN 
END 


Jl 

J2 


K1 

K2 

K3 


LI 

L2 


1 

2 

3 


4- PIXCORI M3) 


GAT0008C 

GAT00090 

GAT00100 

GAT00110 

GAT00120 

GAT00130 

GAT00140 

GAT00150 

GAT 00160 

GAT00170 

GAT 00180 

GAT00190 

GAT00200 

GAT00210 

GAT00220 

GAT 00230 

GAT00240 

GAT00250 

GAT00260 

GAT 00270 

GAT00280 

GAT00290 

GAT00300 

GAT0031 0 

GAT00320 

GAT00330 

GAT00340 

GAT00350 

GAT00360 

GAT00370 

GAT0038C 

GAT00390 

GAT00400 

GAT00410 

GAT00420 

GAT 00430 

GAT00440 

GAT00450 

GAT00460 

GAT00470 

GAT00480 

GAT0C490 

GAT00500 

GAT00510 

GAT00520 

GAT00530 
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1 . 


Module Usage 


FILBUF 

CALL FILBUF (ICLSCV, VECPRD, SCAPRD, SETCOV, CONST, BDAT A , 
RDATA, PIXCOR, PIXVAL, CELCOR, CELSUM, CELIKE, 
CELCLS, PIXCLS , CURLIN, *) 


Input Arguments ; 
ICLSCV 

(VARSZ3, NOPOOL) 


VECPRD 

(N0FET3, NOPOOL) 


SCAPRD (NOPOOL) 


DETCOV (NOPOOL) 


CONST (NOPOOL) 


BDATA 

(NVR, ID (6)) 


RDATA 

(NVR, CELWTH) 


R*4 Class covariance matrices from 
statistics dimensioned classes by 
(channels * (channels + l))/2. 

R*4 Vector product matrix that is 
inverse of class mean matrix trans- 
posed times the covariance matrix for 
the class, dimensioned channels by 
classes . 

R*4 Scalar product of the classes for 
all statistics classes. 

R*4 Determinant of the covariance 
matrix for each class. 

R*4 Constant array computed from In 
(2*PI*class covariance matrix) for 
each class . 

L*1 Buffer array used by GADLIN for 
reading data tape dimensioned channels 
on tape by samples per line on tape. 
R*4 Array used to store calibrated 
data values from tape for one row 
of cells (multiple data lines) . 
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PIXCOR ? 

(VARSZ3, CELWTH ) 


PIXVAL 

(NOFET3, CELWTH; 


CELCOR (VARSZ3) 


R*4 Array used to store cross product 
data values in one cell for a correla 
tion type indication. 

R*4 Array used to store pixel data 
values for a cell. 

R*4 Sum of the PIXCOR values for a 
cell being processed. 


Output Arguments : 
CELSUM (N0FET3) 

CELIKE 

(NWORD, JPTS) 
CELCLS (JPTS) 


PIXCLS 

(NWORDS , JPTS) 


R*4 The sum of all data points in 
a cell for each channel. 

R*4 Likelihood values for each class 
for each cell in a row of cells . 

1*2 Class numbers for each cell in 
a row of cells. 

1*2 Class numbers of pixels for all 


cells in a row of cells. 

CURLIN 1*4 Current line number requested. 

FILBUF is called to fill up CELCLS and CELIKE for annexa- 
tion or writing to the intermediate tape. By either reading 
from the intermediate tape for annexation processing, ' or reading 
from the data tape and carrying out processing for the area, 
cell information can be produced. 


2. Internal Description . 

FILBUF produces arrays CELCLS and CELIKE from either reading 
the intermediate tape that is in a classification format or 
from reading the input data tape via GADLIN . The raw informa- 
tion is then processed by producing cell statistics with GATHER, 
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establishing cell classes with SAMCLS, testing for cell homo- 
geneity, and outputing a row of cell information. This routine 
uses common blocks GLOCOM and SECCOM. 

3 . Input Description 

The only input is possible input from the intermediate 
tape which stores the same format used by classification tapes 
and stores cell statistics for later annexation processing. 

4 . Output Description 

Not applicable. 

5. Supplemental Information 

Not applicable. 


6 . Flowchart 
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FORTRAN IV G LEVEL 20.7 FILBUF DATE * 77133 21008045 

FILE FIL0UF 


0001 

0002 

0003 


0004 

0005 

0006 
0007 

oooa 

0009 

0010 


0011 


0012 

0013 

0014 


0015 

0016 


0017 

0010 

0019 


c 

c 


FIL3UF LARS 


XXXX 


F IL 0001 0 
FIL0CQ20 
^ IL 0CC3 0 
L 0004 0 
L00050 
LOCO60 


SUBROUTINE FIL8UF (ICISCV, VECPRD, SCAPRD, DETCOV, CONST. 
t 8DATA. RDATA • PIXCOR, PIXVAL, CELCGR, CELSUM, CELIKEr CELCLS, 
L PIXELS. CURtIN) 

IMPLICIT INTEGER ♦ 4 (A-Z) 

COMMON /GLGCCM/ BLANK 


C** + * ♦♦•***♦****♦*♦**** t *********** ******** ***** ************ ************f I 

C FI 

C RETURNS CELCLS AND CELIKE VALUES FROM EITHER DATA TAPE OR INTER FI 

C MEDIATE TAPE . FI_ ___ 

C FIL00C70 

€<♦•*♦♦♦*♦♦♦****♦♦*♦♦*♦♦♦*********♦♦*** ********* 4***4 *4*********44******F1LCCQ80 

F IL 0 009 0 
FILC0100 
FILOQUO 
F I L 001 20 
FIL00130 
FI L 00 1 40 
FILC0150 
F ILQ0L60 
FIL00170 
F 2 L 001 60 
FIL0019C 
F 1L00200 
F I L 002 1 0 
FIL00220 
FIL0023C 
F IL0C240 
F I L 00250 
F1L00260 
F IL 002 7 0 
F I L 00280 
F IL 0029 0 
FILCOJCO 
FIL00310 
FIL00320 


1 

2 

3 

4 

5 

6 
7 

a 

9 


CLUSTX. CONPUT, C 
DUPLTP 


, CARD (20 ) . CHKOUT. CQPFIL, CLASSR* CLASSX, 
PYOUT , CRDRDR, CRDSEG, DATAPE. 


bUPRUNi ERRMSG, F6PNT, 

FI LESV . FIDGNO, HDAT A « HEAD(88>, ID(200I V IMAGEX. 
IMARK. KEYBD, MAPTAP. MAXCHA, MAXCLS, 

PA GS I L « PNCH . POINT, PRESUX. PRNTR, READ IN t 
RESTRT . RLNFIL, RUNT AB < 10 . 3 J . 

SDATA, SEPARX, SEPTPX, SPAREUO), TEMPAS430), 
TPSTAT (6) . T TPLDX .TYPEWR, 

TOP, ARRAY t 12500) 


REAL * 8 ARRAY 
REAL * 4 F ROC AL (5,30) 

INTEGER * 4 C0MENTU6I, DATEI5), HED1 ( 16) . HE02U61, TIHEC5) 

INTEGER * 2 BLANK2 
LOGICAL ♦ 4 CttKOUT 

LOGICAL ♦ 1 BLANKl 

EQUIVALENCE (OAT SAV . 10U )> , l CURRUN, I DO ) > , IFRQCAU U ♦ ID* 51 )) • F1L00330 

1 (HED l f 1 ) i HE AD ( 8 ) > • (DATE! 1) .HE AD 126) ) t 1HED21 1 ), HE AD( 39 )), F I L 0C34 0 

2 (TIME! 1 ) .HEAD (58) f , ( COME NT (1 ),HEAD(72) >, FILC035C 

3 (MAPSAV, TPSTATClM , FIL0C360 

4 1 SEP SCR ,TPSTAT(2) I , ( OUP I N.TPST AT ( 3 ) ) . ( OASTAT . TPSTAT ( 4 I) , F I L CC3 70 

5 (CCPSER,TPSTAT(5)), C TRAOUT , T PST AT (6 J 5 , FIL0C38C 

6 (BLANK, BLANK2 ,8LANK1) FIL00390 

F IL 00400 

COMMON /SECCQK/ BUFROZ , CSELEC ,CSET l 3 , 30 \ ,CSET3( 3, 30 ) , THRES ,CELS t l , F IL004 10 

1 CELWTH.INFOd 7) , INPUT, INTTAP, INTFIL, JPTS, LINES, NOCLS, NOFET , F I L0042 0 

2 NOFET3.NOPaOL,N*CRD t NWORD2, NVR, OUTPUT, POLNAM 1 2, 60 ) • PREF IX( 2 ) , F IL 004 3 0 

3 PTS* ROTAPE, RQFILE, SYHCNT, VAftSZ3. CSEL(30). CSEL3I 301* 

3 FETVECC30), FETVC3(30>, POLPTRI 2 ,60 1 , POLSTM&oT, 

4 CDFLAG, C L SNAP • CSETl(3 f 30), OBJMAP, PHASE1, PHASE2, P0LNM1 (60) , 

5 PR S TAT, SYM(60), SYHHTX i 60 l 

REAL ♦ 4 CSELEC , CSET, CSET3, THRES 


INTEGER ♦ 2 CSEL, CSEL3* FETVEC, FETVC3, POLPTR, POLSTK 

LOGICAL ♦ 1 CDFLAG, CLSMAP, CSET1, OBJMAP, PHASE1, PHASE2, PGLNNl 
I PRSTAT 


LOCAL VARIABLES 

INTEGER*2 BLOCK 14 ) ,CE LCLS < JPTS) , PIXCLS (NMCRD2 , JPTS ) , 

1 BDATAU) 

REAL • 4 I CLSCVt VARSZ3*N0P00L J , VECPRD (NOFET 3, NOPOOL It 
I SCAPROINOPCOLltDETCOVlNOPOaLltCOASKhOPOaDtCELSUMlNfiFeTa), 
\> ,P|XC0R( 1 ) , CEL1KE INWORD, JPTS ) » 


* ^ W n p i nwr ViU ^ i wu r I iwr v 

2 CELCOR l VARSZ3 > »P l XVAL1 1 > ,F 

3 RDATAtNVR.CELWTH) .CHISQR, J 


INITIALIZE 

ERRFG * 0 
DO 100 1*2,4 

LOO BLOCK (U - INFQU451 

GLCR * NUMBER OF OATA LINES TO SKIP TO GET NEXT CELL DATA LINE 


F I L 00440 
F I L 00450 
F L 00460 
F L 004 7 0 
F IL 0048 0 
FIL004S0 
FILC0500 
F IL 0051 0 
F I L 00520 
F IL 00530 
F I L 0054 0 
F I L 00550 
F I L 00560 
FILQ0570 
FIL00580 
F I L 00590 
F I L 00600 
F I L 0061 0 
FILOC6ZO 
F IL 00630 
F IL 0064 0 
F I L 00650 
FIL00660 
Fit 00670 
F IL 00680 
F IL 00690 
F I L 00700 
F IL 007 1 0 
FIL00720 
FIL00730 
FIL00740 
FIL00750 
F IL00760 


PAGE OOOl 


-J 
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FILE FIL0UF 


20*7 


FXL8UFC 


DATE - 77133 


21008045 


PAGE 0002 


0020 


0021 

0022 

0023 

0024 

0025 

0026 

0027 

0028 


0020 

0030 


0031 

0032 

0033 


0034 

0035 
00 36 
0037 


0036 

0039 

0040 

0041 

0042 

0043 


0044 

0045 

0046 


0047 

0048 

0049 

0050 

0051 


NSR * NVR/NOFET3 


F IL00770 

- - F I L 0078 0 

C F IL 00790 

*♦♦******« + ♦*#***«♦♦*♦*****+*****# ************* + *******************Fm 00 800 

C FIL00810 

C ENTER MAIN LOOP FIL00820 

C FIL00830 

€♦♦♦*♦*♦****♦**♦♦**#♦+♦*♦♦*•♦♦***♦*******♦# ♦♦♦♦*♦♦********♦** ***********f 1 1 00840 


C GET ONE ROW OF DATA FOR CELLS (MAY BE MORE THAN ONE LINE) 

IF J PHASE2 .AND* .NOT. PHASED GO TO 400 
I « CURLIN 
DO 155 3*1 tCELWTH 
BLOCK t 1 ) » I 

CALL GAOL I N( BLOCK, CSE L3, CSET3 , I Dt DAT APE, NQFET3 ,NSR , 00 ATA. 

C RDATA { 1 , J ) . ROLL » ERROR) 

IF(ERRQR.Gl.O) CALL L I NERR( ERROR, I, C8 10, C 155, £ 6 10 ) 

155 I * I ♦ INF016) 

H - 1 
C 

C GET CELL DATA ARRAYS UPDATED FOR EACH CELL IN LINE 
DO 200 CELNUM* 1 . JPTS 

c C ^e^t:^g^5tSt^^I PlXVA1 ' ,PI,tC0R ’ ceL5UM ' ceLC0R ’ ceLWTH ' 

C CLASSIFY EACH CELL WITH MAXIMUM LIKELIHOOD AND RETURN CHI-SCUARE 
C VALLE FOR THRESHOLDING HOMOGENEOUS DEFINITION 

C CALL SAMCLS(CLASSrCHISQR,CEUKEUtCELNUM) r ICLSCV f VECPRD, 

t CELCLS(CELNUM? E «CLiss LSUM,CeLSIZ,NaFET3,VARSZ3tN0P0 ° i " , * TRUE * 
IF (CHISQR • LE. ^CSELEC ) GO TO 200 


C NON HOMOGENEOUS CELL 

L * 1 
M * 1 

CELCL S (CELNUH) * 0 
CO 196 PlXNUH*l,C£LSIZ 

C TRY AGAIN 


CLASSIFY VECTORS SEPARATELY. 


CALL SAHCLS (CLASS rCH I SQR.Al, I CLSCV,VECPRD*SCAPRD, CONST. 

L PIXCOR(M) , PI XVAU L) , 1 » N0FET3 , V ARSZ3 , NOPOCLt .FALSE.) 

P1XCLS (PI XNUM. CELNUH) ■ CLASS 
t - L ♦ N0FET3 

196 M * M ♦ VARSZ3 

200 H * H ♦ CELWTH 

^ CO TO 210 

C •*♦*♦♦♦♦*♦* mo 1280 

C FIL0I290 

C READ FROM INPUT INTERMEDIATE TAPE FIL01300 

C FIL0I3I0 

****«#***#*« ****** ******* ********** **** ****** ***************p I L 0 l 32 G 

™ F I L 0 1 330 


F1L0O85O 
F I L 00860 
FIL008/0 
FIL00830 
F I L 00890 
F IL 00900 
FIL00910 
F IL 00920 
FIL00930 
F IL 00940 
FIL00950 
F I L 0096 0 
FIL00970 
F I L 0090 C 
F IL 00990 
F IL 0 1 000 
FIL01010 
FIL0I020 
F IL 0 1 030 
F IL 0 1 040 
F I L 0 1 05 0 
FIL01Q60 
F I L 0 1 07 0 
FIL0I080 
FIL01090 
FIL01100 
FILOIUO 
FIL01I20 
FIL0U30 
FIL0I140 
FIL0H50 
FIL0II60 
FIL0U70 
FIL0II80 
F1L01I90 
F1L0I2O0 
FIL0I2L0 
FIL01220 
FIL01230 
F I L 0 1 240 
FIL01250 
FIL0126C 
FIL0I270 


400 REAO UNPUT, £RR*01O) PREFIX. I. CELCL S.CELIKE 
. IF U .LI. CURLIN) GO TO 400 
210 RETURN 


^ . FIL01360 

^ ^ — — - — — — L0I370 


ERROR PROCESSING 


L 0 1 340 
L01350 


810 WRI TE( PRNTR f 9810) CURLIN 
9810 FORMAT ( * ORE Jd ERROR FROM GAOL IN* 
C ■ PROCESSING TERMINATED* * ) 

ERRFG * 1 
GO TO 210 
END 


CURRENT LINE NO. IS S 14// 


FIL0138C 
FIL01390 
F I LO 1400 
♦FIL01410 
F I LO 1420 
FR01430 
F IL 0144 0 
FIL01450 
F I LO 1460 
FIL01470 
FILO140O 


CO 
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Author ; P» D. Alenduff Date : 4/20/77 
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MODULE ABSTRACT 

STATS 2 contains entry points LIKRAT, ADD2 , and BEGFLD. LXKRAT 
produces a likelihood test by subtracting log - likelihoods 
and returns a resultant likelihood value. ADD 2 adds a cell 
to the statistics of a field. BEGFLD initializes field sta- 
tistics from a fields statistics. 
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1. ‘ Module Osage 


LIKRAT 

CALL LIKRAT (STAT, CELIKE , FLDLIK , AUXLIK, CELCLS, FLDCLS , 
AUXCLSi) 

Input Arguments : 


CELIKE (NOPOOL) 


FLDLIK (NOPOOL) 


CELCLS 


FLDCLS 


R*4 Log likelihood values for this cell 
for each class in the statistics. 

R*4 Log likelihood values of all cells 
annexed to the field for each class 
in the statistics. 

1*2 Class number of cell being con- 
sidered. 

1*2 Class number of field being tested. 


Output Arguments : 

AUXLIK (NOPOOL) R*4 Resultant log likelihood of an- 

nexation of cell to given field for 
each class. 

AUXCLS 1*2 Resultant class of annexation 

test for this cell. 

STAT R*4 Resulting statistic concerning 

annexation performance for this cell 
to the field being considered. 

LIKRAT produces a resulting likelihood, class and sta- 
tistic if the given cell is annexed to the supplied field. 


ADD 2 

CALL ADD 2 (FLDLIK, AUXLIK, FLDCLS, AUXCLS, FLDSIZ) 
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Input Arguments : 
AUXLIK (NOPOOL) 


AUXCLS 


FLDSIZ 


R*4 Likelihood values formed with 
temporary annexation of cell to field. 
1*2 Class number of cell when annexa- 
tion of cell to field was successful. 
1*4 Number of pixels included in the 
given field. 


Output Arguments : 

FLDLIK (NOPOOL) R*4 Likelihood values of new field 

formed. 

FLDCLS 1*2 Class of field formed. 

FLDSIZ 1*4 (see above) . 

ADD2 adds the cell supplied to the field information sup- 
plied. Temporary likelihood values are assigned as new like- 
lihoods for field. 


BEGFLD 

CALL BEGFLD (CELCLS, FLDCLS, CELIKE, FLDLIK, FLDSIZ) 
Input Arguments : 

CELCLS 1*2 See above. 

CELIKE (NOPOOL) R*4 See above. 


Output Arguments : 
FLDCLS 

FLDLIK (NOPOOL) 
FLDSIZ 


1*2 See above. 
R*4 See above. 
1*4 See above. 


BEGFLD initialize FLDCLS, FLDLIK, and FLDSIZ to be equal 


to CELCLS, CELIKE, and CELSIZ 
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2. Internal Description 

See above. 

3. Input Description 

Not applicable 

4. Output Description 

Not applicable 

5. Supplemental Information 

Not applicable 
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Flowchart 

LIKRAT 


ADD 2 



RETURN 
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STATS2 


DATE « 77X33 


21008002 


PAGE C001 


FILE STATS2 


0001 

0002 

0003 


0004 

0005 

0006 

0007 

0008 
OOC9 

0010 

0011 

0C12 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 
002 3 

0024 

0025 


0026 

0027 

0028 
0029 
00 30 

0031 

0032 


C SIATS2 LARS XXXX ECH00010 
C ECH0002 0 

c ECH00040 
C THIS SUBROUTINE WAS WRITTEN TO HELP IMPLEMENT A MORE COMPLEX SCANNING ECHCC05 0 
C LOOP THAN THE CNE IN THE ORIGINAL ECHC* ESSENTIALLY IT DOES W FAT ECH0CC60 
C ANNEX OCESi BUT IT OPERATES ONLY IN SUPERVISEE MODE OF ANNEXATION. ECHCCCfC 
C ECHCC08C 
C***« ****♦♦*♦♦**♦♦*♦♦*♦♦#♦♦♦*#*♦* ♦**♦♦• **♦*♦♦♦+#*♦♦♦**♦*♦*♦♦*♦*♦ ♦♦♦♦♦♦♦♦ECHCCC90 


SUBROUTINE STATS2 
IMPLICIT INTEGER * 4 <A-Z> 


ECHCOIOO 
ECHOOl I 0 
ECH0C120 
ECHOOl 30 

COMMON /SECCOM/ BUFROZ,CSELEC*CSET< 3. 30), CSET3I3* 30 I , THRES*CELS!Z ,ECH0C140 

1 CELWTH*INF0U7) , INPUT, INTTAP, IMFILt JPTS, LINES* NOCLS# N0F6T* ECHCOl 50 

2 NOFE T3» NO POOL , NWCRD»NW0RD2 , NVR, OUTPUT, POLNAM ( 2 , 60 ) * PREF IX < 2 > , ECHOO 160 

i Hh e ri»RviwM5.CTri».Jrtiai, c sfei«Uo« k »™'' pus; 

4 CDFLAG, CLSMAP, CSETM3.30J, CBJMAP* PHASEl* PHASE2* PQLNK1I60I* ECH0C19C 

5 PRSTAI* SYMI60)* SYMMTX 160) ECH002C0 

ECH002 l 0 

REAL * 4 CSELEC, CSET, CSET3, THRES ECH00220 

INTEGER * 2 CSEL* CSEL3, FETVEC, FETVC3, POLPTR, POLSTK ECH00240 

ECH00250 

LOGICAL » l CDFLAG* CLSMAP, CSET1, OBJMAP* PHASE!* PHA5E2, POLNMl , ECH0C260 
1 PRSTAI 

REAL ♦ 4 ANN1 


" REAL * 4 $TAT,CELIKE( NO POOL ) * FLDL IKI NGPOOL 1 v AUXLI K I NO POOL) 

INTEGfcR*2 CELCLS, FLOCLS, AUXCLS 

ENTRY UKRAT(STAT,CELIKEfFLOUXf AUXL l K* CELCLS* FLOCLS* AUXCL 5) 
CO 20 1*1, NCPOOL 

AUXLIKII) - CELIKEH) ♦ FLOLIKUI 
IF I l . CQ. 1) GO TO 10 
IFMUXLIKm • L E. STAT) GO TO 20 
10 STAT = AUXLlKtn 
AUXCLS * I 
20 CONTINIE 

STAT * STAT - CEL IKE (CELCLS \ - FLCLIK (FLOCLS > 

RETURN 

P* 

ENTRY ADD2 < FLOL I K, AUXLXK, FLOCLS* AUXCLS t FLOS IZ I 
50 DO 60 1=1, NCPCOL 
60 FLDLIKU! * AUXUK(I) 

FLOCLS ■ AUXCLS 

FLDSIZ - FLDSIZ ♦ CELSIZ 

RETURN 

Sentry for field initialization. 

ENTRY 6EGFLD (CELCLS* FLOCLS* CEL IKE* FLOL IK* FLDSIZ > 

FLOCLS - CELCLS 
FLDSIZ » CELSIZ 
DO 40 IpI * NOPODL 
40 FLOLIKUI - CEL1KEU) 

RETURN 

END 


ECH00270 
EC H 00280 
ECH00290 
ECH00300 
EC H 00310 
ECHQ0320 
ECH00330 
ECH0C340 
ECH0035C 
ECH0C360 
ECH00370 
ECH0C38 0 
ECH00390 
ECHOOAOO 
ECH00410 
ECHQC42 0 
ECH0Q43Q 
6CH0044 0 
ECH00450 
ECH0C460 
ECH0C47G 
ECH00480 
ECHQ0490 
ECH00500 
ECH0051 0 
ECH0052 0 
ECH00530 
ECH00540 
ECHO 0550 
ECH00560 
ECH0057 0 
ECH00580 
ECH00590 
ECH00600 
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LARS Program Abstract 

MODULE IDENTIFICATION 

Module Name t SAMCLS Function Name: SECS UP 

Purpose: Classify cell using maximum likelihood 

System/Language: CMS /fortran 

Author: p » D - Alenduff Date: 4/14/77 

Latest Revisor: Date: 


MODULE ABSTRACT 

SAMCLS performs a maximum likelihood classification by calling 
LOGLIK and finding the maximum likelihood class for a cell. 


PURDUE UNIVERSITY 

Laboratory for Applications of Remote Sensing 
1220 Potter Drive 
West Lafayette, Indiana 47906 
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1. Module Usage 


SAMCLS 

CALL SAMCLS (CLASS, CHISQR, LIKELY, COVMTX, VECPRD, SCAPRD, 
CONST, SAMCOR, SAMSUM, SAMPTS, NOFET3, MTXSIZ, 
NOPOOL, RETLIK) 

Input Arguments : 

COVMTX R*4 An array dimensioned (channels 

(VARSZ3 , NOPOOL) 

* (channels + l)/2) * number of 
classes, containing the upper half 
of the triangular covariance matrices 
for each of the classes used for 


VECPRD 

(N0FET3, NOPOOL) 


SCAPRD (NOPOOL) 


CONST (NOPOOL) 


SAMSUM (NOFET3) 

SAMPTS 

NOFET3 

MTXSIZ 


classification. 

R*4 An array dimensioned channels 
* classes containing the vector 
product of the classes for each 
channel from the inverse covariance 
matrix and the determinant of the 
covariance matrix. 

R*4 An array dimensioned classes for 
the scalar product of the mean and 
^covariance matrices for each class. 
R*4 An array dimensioned "classes" 
words containing the constant term 
used consisting of an expression. 

R*4 Sum of data values for all points 
in cell for given channel. 

1*4 Number of points per cell. 

1*4 Number of channels. 

1*4 Number of words needed to store 
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NOPOOL 

RETLIK 

Output Arguments : 

CLASS 

CHISQR 

LIKELY (NOPOOL) 


half of triangular covariance matrix. 
1*4 Number of classification pools. 

L*4 Indication of whether or not like- 
lihoods are to be returned in LIKELY. 

1*4 Class number of cell. 

R*4 Resultant Chi-square value from 
LOGLIK and classification. 

R*4 Log likelihood of cell belonging 
to pools. 


2. Internal Description 

SAMCLS calls LOGLIK once for each class and is returned 
a log likelihood for the cell for each class. If the likeli- 
hood is greater than the previously stored, the new likeli- 
hood is saved with the class number and, when done, this class 
and value are returned. 


3. Input Description 

Not applicable. 

4 . Output Description 

Not applicable. 

5 . Supplemental Information 


Not applicable 
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MODULE IDENTIFICATION 


LARS Program Abstract 


Module Name : LOGLIK Function Name: SECSUP 

Purpose : Produce Log likelihood of a cell 

System/Language : cms/fortran 

Author : P« D» Alenduff Date : 4/13/77 

Latest Revisor: Date: 


MODULE ABSTRACT 

LOGLIK produces a log likelihood for a class in the input 
statistics . 


PURDUE UNIVERSITY 

Laboratory for Applications of Remote Sensing 
1220 Potter Drive 
West Lafayette, Indiana 47906 
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1. Module Usage 
LOGLIK 

CALL LOGLIK (LIKELY 
DIMEN, 

Input Arguments : 
ICLSCV (MTXSIZ) 

VECPRD (DIMEN) 

SCAPRD 

CONST 

DIMEN 

SAMCOR (MTXSIZ) 

SAMSUM (DIMEN) 
SAMPTS 

Output Arguments : 

LIKELY 

CHISQR 


CHISQR, ICLSCV, VECPRD, SCAPRD, CONST, 
SAMCOR, SAMSUM, SAMPTS) 

R*4 Class covariance matrix from 
statistics, dimensioned channels * 
(channels + l)/2. 

R*4 Vector product of class covariance 
matrix determinant and means dimensioned 
"channels" words. 

R*4 The scalar product term for the 
class being considered. 

R*4 The constant term derived in SECPRT. 
1*4 Number of channels of data in 
cell information. 

R*4 The cross product terms from the 
cells data values dimensioned as 
ICLSCV. 

R*4 The sum of the data values for 

each channel for the cell. 

1*4 Number of data points in a cell, 

2 

actually (cell width) . 

R*4 Log likelihood value for the cell. 
R*4 Chi-square value returned for the 


cell 



LOGLIK-3 


LOGLIK is called to return a log likelihood and a Chi- 
square value for the cell being considered. 

2. Internal Description 

The Chi-square value is calculated by the equation: 
CHISQR = trfCj” 1 I Y ± Y i ’ ) - 2M_.' C^ 1 ? + mM.. ' C^ 1 Mj 

and the log likelihood is calculated by: 

LIKELY = -.5 (m • ln|2TT C^ | + CHISQR) 

Where : 

Cj is the class covariance matrix for class j . 

M j is the mean vector for class j . 

Y^ is the channel response vector for pixel i of the cell, 
m is the number of pixels in the cell. 

3 . Input Description 

Not applicable. 

4 . Output Description 

Not applicable. 

5. Supplemental Information 


Not applicable 



6* Flowchart 
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FDR TRAN IV G LEVEL 20.7 


LGGLIK 


DATE = 77129 


15025029 


PAGE 0001 


FILE LOGLIK 


0001 


0002 

0003 

Q0G4 

0005 


0006 


0007 

0008 
00C9 
0010 
0011 
0012 
0013 

0014 

0015 


0016 

0017 

0018 
0019 


C LOGLIK LARS XXXX LOGOOOiO 
C L0G00020 
C***** ********************************* 4************************ ********LQGC0030 


SUBROUTINE LOGLI K ( LI KELY* CHISQR, ICLSCV, VECPRD, SCAPRD, CONST, DIMEN , 
£ SAMCOR , SAMSUM, SAMPTS) 

LOCAL VARIABLES 

DIMEN = NUMBER OF CHANNELS USED 
SAMPTS » NUMBER OF PTS PER CELL 

CONST = LN(2*PI*CIN)) WHERE C(N) IS COVARIANCE MATRIX 
FOR CLASS N { N IS DETERMINED BY SAMCLS) 

SAMSUM = SIM OF Ytl) 

SAMCOR = SUM OF YfI)*Y<I> TRANSPOSED 
ICLSCV = C!J> ** -1 _ t 

SCAPRD = M!J)*C(J> TRANSPOSED *M(J) WHERE C(J) = CLASS 
COVARIANCE MATRIX AND M(J) = CLASS MEAN MATRIX 
VECPRD = M ( J ) TRANSPOSE * C(J) ♦* -1 


LDG00040 
LOGOOOSO 
L0G0C06C 
LOG0007 0 
L0G00080 
L0G00090 
LOGOOiOO 
LOGOOllO 
LQG00120 
LOG0013C 
L0G00140 
LOGQ0I5O 
L0G00160 
LOGOOITO 
LGG00180 
L0G0019G 
L0G0020C 
L0000210 


C *************** ************** ************************* *****************^ogo 0220 


IMPLICIT 
REAL * 4 


K = 1 
LIKELY 
C 

C SET CHISQR 
C 

CHISQR 

lCORll) 


INTEGER * 4 <A-Z) 

L IKELYtCHISQR, ICLSCV (1 ) t VECPRD! DIMEN ), SCAPRD , CONST, 
SAMCOR ( 1 ) , SAMSUM! DI MEN) 


LOG00230 
L0G00240 
L0G0025C 
L0G00260 
LDG00270 
L0G00280 
L0G00290 
L0G00300 
LOG00310 

» .5 * SAMPTS ♦ SCAPRD-VECPRD! 1)*SAMSUM< 1 ) + . 5* ICLSCV t 1 )*SAMLQG00320 

L0G0033C 


- 0 

= .5 * NO. OF POINTS/CELL * LN ( 2*P I*OET ( CL S COVAR MIX)) 


C 

C AT 

C 

C 


DO 2 I 
CHISQR 
L = I - I 
DO 1 J*1,L 
K = K+l 

CHISQR = CHISQR 
K = K ♦ 1 
LIKELY = LIKELY 
CHISQR = CHISQR 


2 ,DI MEN 

CHISQR - VECPRD! I)*SAMSUM( I) 


+ I CLSCV ( K ) *S AMCOR! K ) 


ICLSCV! K ) * 
•5*LI KELY 


SAMCOR ( K ) 


THIS POINT , CHISQR = ,5 * SAMPTS * SCAPRO + *5 * TRACE (ICLSCV* 

SAMCOR) - VECPRD* SAMSUM 

LIKELY * -SAMPTS*CONST - CHISQR 
CHISQR = CHISQR + CHISQR 
RETURN 
END 


L0G00340 
LQG00350 
LOG0036C 
LOG00370 
L0GQG380 
LQG00390 
LOG0040C 
LQG0041 0 
L0G00420 
L0G00430 
LDG00440 
L0G00450 
L0G00460 
L0G00470 
L0G0048 0 
L0G00490 
L0G00500 
LQG00510 
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MODULE IDENTIFICATION 


LARS Program Abstract 


Module Name: SECRD Function Name: SECSUP 

Purpose s Read data arrays from intermediate tape, 

System/Language : CMS/Fortran 

Author : p « D * Alenduff Date . 4/12/77 

Latest Revisor: Date: 


MODULE ABSTRACT 


SECRD reads data arrays from the intermediate tape for 
annexation processing by SECHO. 


PURDUE UNIVERSITY 

Laboratory for Applications of Remote Sensing 
1220 Potter Drive 
West Lafayette, Indiana 47906 
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1. Module Usage 


SECRD 


CALL SECRD (COVMTX, AVEMTX) 
Input Arguments : 


COVMTX 

(VARSZ3, NOPOOL) 


AVEMTX 

(NOPET 3, NOPOOL) 


R*4 An array dimensioned (channels 
* (channels + l))/2 fullwords for 
storing the covariance matrices of 
classes used in classification. 

R*4 An array dimensioned channels 
by classes words for storing the 
class means of classes used in the 
classification. 

Output Arguments : 

Same as above. 

SECRD is called to read one reocrd from the intermediate 
tape. COVMTX and AVEMTX are read in an unformatted fashion 


from the unit INPUT. This program uses common blocks GLOCOM 
and SECCOM. 

2, Internal Description 

See above. 

3. Input Description 

SECRD reads from unit INPUT a record type 4 from the 
classification format file. 


Output Description 


Not applicable 
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5. Supplemental Information 

See LARSYS Systems Manual fox format of classification 
results file. 
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6 . Flowchart 



FORTRAN IV G-LEVEL 20.7 


SECRD 


DATE = 77129 


15024017 


PAGE 0001 


FILE SECRO 


C SECRD LARS XXXX SEC00010 

r SEC00020 

r****************#*****+******************4**4***+*’M*******************SEC0003C 
C SEC00040 

C SECPRD ROUTINE TO READ DATA FROM INPUT INTERMEDIATE TAPE SEC0QO50 

C SEC00060 

r *********** *************************** **** **** ****** *********** ********SECCJ0C7C 
C . SEC00080 

SUBROUTINE SECRD (COVMTX, AVEMTX) SEC00090 

IMPLICIT INTEGER * 4 ( A— Z > SEC00100 

C SEC 001 10 

COMMON /GLOCOM/ BLANK. CARD120), CHKOUT, COPFIL, CLASSR. CLASSX. SEC00120 

1 CLUSTX, CONPUT, CPYOUT, CRDRDR, CRDSEQ* OATAPE, SEC00130 

2 DUPLTP . DUPRUN. ERRMSG. FBPNT . SEC00140 

3 FILESV » FLDBND, HDATA, HEAD188), IDJ200), IMAGEX, SEC00150 

4 IMARK, KEYBD. MAPTAP, MAXCHA, MAXCLS. SEC00160 

5 PAGSIZ, PNCH, POINT, PRESUX. PRNTR, READIN, SEC00170 

6 RESTRT, RUNFIL, RUNTABI 10 ,3) , , SEC00180 

7 SDATA, SEPARX. SEPTPX, SPARE ( 10 ) » TEMPAS(30>, SEC00190 

8 TP STAT (6) » TTFLDX .TYPEHR, SEC00200 

9 TOP, ARRAYU2500) SEC00210 

C SEC00220 

REAL * 8 ARRAY SEC00230 

REAL * 4 FRQCAL(5»30) SEC00240 

INTEGER * 4 C0MENTU6), DATEI5) , HE0K16), HED2116), TIMEC5) SEC0025G 

INTEGER * 2 BLANK2 SEC00260 

LOGICAL * 4 CHKOUT SEC00270 

LOGICAL * 1 BLANK1 SEC00280 

EQUIVALENCE (DATSAV.IDCl)) * ( CURRUN, I D<3 )) , (FRQCALt l )• 10(51) ), SEC00290 


COMMON /GLOCOM/ BLANK, CARD120), CHKOUT, COPFIL, CLASSR, CLASSX, 

1 CLUSTX, CONPUT, CPYOUT, CRDRDR, CRDSEQ, OATAPE, 

2 DUPLTP, DUPRUN, ERRMSG, FBPNT, 

3 FILESV, FLDBND, HDATA, HEADI88), 10(200), IMAGEX, 

4 IMARK, KEYBD, MAPTAP, MAXCHA, MAXCLS, 

5 PAGSIZ, PNCH, POINT, PRESUX. PRNTR, READIN, 

6 RESTRT, RUNFIL, RUNTABI 10 ,3) , 

7 SDATA, SEPARX, SEPTPX, SPARE ( 10 ) » TEMPASC30), 

TP STAT (6) , TTFLDX .TYPEHR, 

TOP, ARRAYU2500) 


IMAGEXi 


REAL * 8 ARRAY 
REAL * 4 FRQCAL(5,30) 

INTEGER * 4 C0MENTU6), DATEI5), HE0K16), HED2116J, TIME(5> 
INTEGER * 2 BLANK2 


* 4 CHKOUT 

* 1 BLANK1 

NCE (DATSAV, ID ( 1 ) ) • 
(HED1 ( 1 ) ,HEAD(8)> . 
(TIME(1 ) .HEAD (58 J I , 
(MAPSAV.TPSTAT < 1 ) ) , 
t SEPSCR ,TPSTAH2 1 ) , 


(CURRUN, IDO) ), (FRQCAL(I), I0(51>), SEC00290 

(DATE( 1 ) * HEAD( 26) ) , (HED2( 1 )»HEAQ(39) ) , SEC 0030 C 
( CQMERTi 1 ) ,HEAD( 72 ) ) , SEC0031C 

SEC00320 

(DUPIN,TPSTAT(3) ), ( DASTAT.TPSTAT ( 4 ) ) , SEC00330 


ISEKbGK.f P5TAMZI ) , ( DUP I N • TPSTAT ( 3 ) ) , 

(CGPSER,TPSTAT(5) ) , < TRAOUT.T PSTAT ( 6 > ) , 
( BLANK, BLANK2, BLANK!) 


SEC00340 

SEC00350 

SEC00360 

SEC00370 


COMMON /SECCOM/ BUFROZ.CSELEC.CSET ( 3, 30) , CSET3 ( 3,30), THRES ,CELSIZ, SEC 0001 0 
1 CELHTH.INFQin), INPUT, INTTAP, INTF1L , JPTS, LINES, NOCLS, NOFET.SEC00020 


N0FET3»NQP00L,NH0RD» NHQRD2 » NVK, UUi kui,_ kulnahi pu », rncruu » 
PTS, RQTAPE , RQFILE , SYMCNT, VARSZ3 , CSELi 30 ) , CSEL3 ( 30 ) , 
FETVECOO), FETVC3 (301, POLPTR< 2 ,60 ) , POLSTK(60)» 

CDFLAG, CLSMAP, CSET1(3,30), OBJHAP, PHASEI, PHASE2, P0LNM1 ( 60 ) , 
PRSTAT, SYM(60), SYMMTX (60) 


TF1L , JPTS, LINES, NOCLS, NOFET.SEC00020 
OUTPUT, POLNAM( 2,60 ), PREFIX( 2) .SEC00030 
Z3. CSEL( 30 ) , CSEL3 ( 30 ) , SEC00040 


REAL * 4 

INTEGER 

LOGICAL 
I PRSTAT 


CSELEC, CSET, CSET 
2 CSEL, CSEL3 , FE 
1 CDFLAG, CLSMAP, 


CSET3, THRES 


FETVEC, FETVC3, POLPTR, POLSTK 
P, CSET1, OBJMAP, PHASEI, PHASE2, 


C OVMTX ( VARSZ3 , NO POOL ) , AVEMTX ( NO FET 3, NOPOOL ) 


READ (INPUT) 

RETURN 

END 


PREFIX, COVMTX, AVEMTX 


SEC0005C 
11(60), SEC00060 
SEC00070 

secoooao 

SEC00090 
SECOOIOO 
SEC00I10 
SEC00120 
P0LNM1 , SEC00130 
SEC00140 
SEC00150 
SEC00150 
SEC00530 
SEC 00540 
SEC00550 
SEC 00560 
SEC0057C 
SEC00580 
SEC00590 
SEC00600 



LARS Program Abstract 

MODULE identification 

Module Name: SECPRT Function Name: SECSUP 

Purpose: Statistics printing and output routine for SECHO 

System/Language : CMS / F0RTRAN . 

Author : p » D » Alenduff Date : 4/12/77 

Latest Revisor: pate: 

• " 1 

- . - .. ■ ■ 11 M 1 1 f 11 m 

MODULE ABSTRACT 

SECPRT produces tabular statistics for user information as 
well as carrying out minor data array processing. 


PURDUE UNIVERSITY 

Laboratory for Applications of Remote Sensing 
1220 Potter Drive 
West Lafayette, Indiana 47906 


' SECPRT-3 


each channel. Then the information concerning classes used 
is printed and WRTMTX is called to print the correlation ma- 
trix, The determinant of the covariance matrices is formed 
by a call to SAMINV, and a check is made to determine if any 
determinants were zero or negative. Then SMMULT is called 
to produce the array SCAPRD, the scalar product, and a calcu- 
lation is performed to produce the constant array. Execution 
is terminated if any of the determinants were zero or negative. 
This routine uses common blocks GLOCOM and SECCOM. 

3 ♦ Input Description 

Not applicable. 

4 . Output Description 

A record type 4 of the classification format output file 
is written to unit OUTPUT. 

Also, statistics information is produced on unit PRNTR. 

5. Supplemental Information 


Not applicable. 
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1 . 


Module Usage 


SECPRT 

CALL SECPRT (COVMTX, AVEMTX , SCAPRD, DETCOV, CONST, WORK) 
Input Arguments : 

COVMTX R*4 Covariance matrix for each 

(VARS2 3, NOPOOL') 

class . 

AVEMTX R*4 Mean matrix for each class. 

(N0FET3, NOPOOL) 


Output Arguments : 


SCAPRD (NOPOOL) 
DETCOV (NOPOOL) 


CONST (NOPOOL) 


WORK (NOFET3 * 
(NOFET3 + 2) ) 


R*4 Scalar product for each class. 
R*4 Determinant of the covariance 
matrix for each class . 

R*4 Constant formed for computa- 
tion in classification for each 
class . 

R*4 Buffer array for inverting 


covariance matrix at least chan- 
nels by channels + 2 words. 

SECPRT is called to print the statistics information to be 
used for classification, if requested, and produce the record 
type 4 of a classification results format file for SECHO. Also 
the determinant, scalar product, and constant arrays, DETCOV, 
SCAPRD, and CONST, are produced for later usage. 


2. Internal Description 

Header information is printed on the Unit PRNTR followed 
by the writing of record type 4 on the unit OUTPUT. This 
classification format record contains COVMTX and AVEMTX, the 
covariance and mean matrices used in the classification for 



Flowchart 
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FORTRAN IV 0 LEVEL 20. ? 
FILE SECPRT 


SECPRI DATE - 77129 15024034 


0001 

0002 

0003 


0006 

0007 

0000 

00C9 

0010 


0011 


0012 

0013 

0016 


ooia 

0019 


SECPRT LARS OGOO 

ROUTINE TO PRINT INFO FOR SUPERVISED ECHO £nD WRITE INITIAL RECORDS 


SUBROUTINE SECPRTlCOVMTX, AVEHTX.SCAPftC.OETCOV.CONST, 

Implicit integer * a u-zi 


" , ?f^!i“?£);p5vf' , ?fYgs?? , ?s!.8oS™?a! E r Eivj..ii* s sr> ci * ssx - 

FllESvJ REBfeWJBKil. 10(2001, IWAGEX, 

J NARK. KE YBD* KAPIAP. NAXCHA, MAxCLS, 


REAL * 8 ARRAY 

INTEGER 4 * V CCK£nI( 16), OATECSl* HEDl ( 161 , HED2U6), TlNETSl 
INTfcCtR * 2 BLANK? 

LOGICAL * 6 CHKOUT 
LOGICAL * 1 8LANK1 

LUUIVALENCfc IO&TSAV. ID! 1 1 » , ( CLIR RON » I Ct 3 > ) , (FRGCALU ) . 10(51 1) , 

I IHE01 I LI ( HEAO 18 1) , (DATE! I1«HEAD(2<1 1* IH€02( lit HE AO (391 ) , 

i ««M.€NTIll,HMO»ill. 


IClPScR.TPA IAM3I l , If 
(BLANK ,SLANK2 tBLANXl 1 


(DASTAT.TPSTAT16)), 


COMMON /SECCCH/ 8UFft02,CSELCC»CSen 3,101 ,C$ET3(3. JO),ThRE$,CELSIZ, 
1 CELmTHiINFOIIH .INPUT, 1NITAP. INTFlL, JPTS, LINES* NOCLS, NOFET * 

! *?S“fa7I8“*8i«8: K s “?KSf , “;S6 s S^T p gIt L I f * x ‘ 2 1 ' 

3 FElvEC(301, FEIVC3I30I. PCLPIRI 2.60 1 * POLSTKUOl, 

*, CDFLAG* CLSHAP, CSETI ( 3 » JO I , PBJFAP, PHASE! r PHA$E2, PCLNH1 [601 « 

5 PRSTAU $YM(60), SVMHTX1601 

REAL ♦ 6 CSELEC. CSET, CSE13. THRES 

INTEGER * 2 CSEL, CSEL3* FETVEC, FEIVC3, PCLPTR, POLSTK 

CL5MAP, CSET 1 * OejXAP, PHASE1, PHASE 2* POLNMl, 


LOCAL VARIABLES DEFINITION 

REAL * 6 COVKTXCVARSZ3.NOPaOLl*AVEKTXINOFET3.NOPGGLJ,VORKm. 
C SCAPROINOPOOL) « DETCOV(NDPCDL) .CONST (NCPGCL) * ALOG.SQRT 
INTEGER * 6 BCDTkiO.CNT.LOCtSTCP 
DATA 0COTNC /*2»/ 

LOCAL VARIABLES DESCRIPTION 

AVEMIX > REDUCED MEAN MATRIX 

CNT * PAGE COUNT 

COVHTX - REDUCED COVARIANCE MATRIX 

LOG ■ INDEX POINTER 

STOP « index pointer,. _ 

WORK - UTILITY VECTOR OF NOF£T3*<N0f ETJ»2I WCROS 
PRINT CLASSES AND CHANNELS CONSIDERED 


PAGE 0001 
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FILE SECPftT 


OATE * 77129 


0022 

0023 

0026 

0025 

0026 
0027 


0030 

8811 


885! 

885! 

0039 
0060 
0Q6 1 
0062 

Soil 

0065 

0066 
D06T 
0068 
0069 

8818 

0052 

881 ! 

0056 

0057 

0058 

882S 

0062 

0063 

006* 

0065 

0066 

0067 

0068 
0069 


0072 

0073 
0076 
0075 


0076 
00 77 
0078 


NRITE IPRH(R, 91001 

9100 format (// rr, ‘T raining classes* *T66, 

1 IT.IM'-M, t66 ,24 I *-• i 

1 . _ //TI2, 'CLASS* *T 45, 

2D* , T96 ,«CAL. CODE* I 
J - NAXOINCPOCL * NCFEf 3 > 

DO 150 I * I * J 
WRITE (PRNTR.91501 

9150 FORMAT I* •» 

IF II LE Nl 

iso if ii.le r 

l FRQCALI2 r 

9151 FORMAT I *+*»Tll,2A4i 

9152 FORMAT r*'* 


CHANNELS FROM STATISTICS*/ 
CHANNEL NO* * »T70» 


SUP0076C 
SUP0C770 
SUPCJ076C 
SPECTRAL BANSUP0C79C 
SUP0080C 
SUPOOfllC 
SUP 00 620 
5UP00830 

NDPOOLI WRITE (PRNTR.9I511 POLNAMT 1, 1 l.POLNAHI 2, M luPOOBSO 


*T51, 12*169. F6. 2, 3X,F6*2,T95. 151 
C WRI IE RECORD TYPE 4 


AVERT X 


PREFIXIIl » * 

WRI It ioutputi prefix, covmtx 

IF I NOT-PRSTATI GO TO 200 
C PRINT MEAN AND COVARIANCE 

sf : %!a«Ar 4 °* NQFCT3,#nNa,:6Tj4 “ ,/i21 

INC « CNT _ 

00 1B0 ICLAS « 1 , NOPOOL 
JK * O 

DO 179 J * l.NOFCU 

WORK! J 1^*^ SORT ICO VMT X ( JK4J T JClA$| 1 

JX ■ JK*I ' 

W0RKINDFET3+JK) » O. 

IF (WQRK(KJ*WORK(J) .LT. I -0E-2S 1 GO TO 179 
179 C0N^NLe E:t3 * JK, " CCVMiXCJK.ICLASl/lWCRKI JMWORKIKM 
IF I INC.Lf.CNTl GO TO 160 
CALL GETIMEITJMEI 


5UP009JO 

SUP0D94C 

SUP0O95C 


5{tN c !KS«iH6», 

INC * 0 


160 WRITE (PRNTR.91601 f POLNAMI I , 
9L60 FORMAT I *0*/IodLASS.... • ,2A*J 
CO 170 tOC « 1.NQFET3.I2 
STOP * M|NOINOfET3,LOC*ll) 


(POLNAMI I. JCLASM-1,2) 


_ „.3,L0L 

VRlth IPRN1R, 91611 (FETVC3I I ) * I * LDC.STOPJ 
9161 FORMAT I'O'/' CHARNEL*,!?, lllil 

IP8NTR.9I62I IFROCALtl, FETvCJl 1 11 , 


1.R 


9162 FOI 
5163 FUI 

Wl„ 

9)64 FORMAT 
WRITE 

9165 FORMAT , 

170 WRITE IPRNIR, 91661 

9166 FORMAT I 


91631 IFRQCAL(2,FEFVC3( 111. 
,CTRAL',12(F6.2,* - Ol 
BAND* ,6X»l2(Fb.* 


RAT llXt'BAND* ,6X »I2(F6.£,3X)) 

• icc.jropi 

iir 'imi'&lL !;&!!{ f^,: 2 v cc ' s,QP ‘ 


WRITE (PRNTR, 91701 
9170 FORMAT (*0*/« CORRELATION MA 
CALL WRTMT*(WCRMNCFET3Fll,h 


180 INC - INC «■! 


C COMPUTE INVERSE OF COVARIANCE, DETERMINANT * VECTOR AND SCALAR PRODUCTsf UPOI 370 
200 DO 210 I • UNOPOOL 

81V8 6 ?c 4 SW?! J ?!8! x 8i i '!4 N ?| E J ! j 0E ,MV 1 " • - cax ' 

2 .0 |Sfelr?': 7 u t T iiSTEJS!U!lASrIsi 5 iAfiT 4 Fl?iElf«ii??? ,c ’ raFe,i * 

c COVARIANCE matrix is SINGULAR - PROCESSING TERMINATED 

C TERMINATE OUTPUTS TAPE 
220 CALL ERPRHT (160**00*) 

WAITE I PANTR, 92221 I 

9222 FORMAT (• SINGULAR FOR CLASS *,12) 


00 c 
.010 
1020 


"8 


i?8 

ZBC 

290 

300 

310 

.32C 

1330 


360 
390 
400 
.410 
42 C 
430 
440 
450 
46C 
47 C 
ABO 
470 
500 


PAGE 0002 


SECPRT-6 


FORTRAN IV C LEVEL 
FILE SECPftT 


CATE » 77129 


i!u*TPuV* SSR ’ 


V NO - F TLESVt I 

lift kin isskhp 

CALL TDPBF (OUTPUT 
CALL RJKAIN 
RETURN 


ROTAP6, NO, 13, K. SERIAL, (10*1 » l*7» 



PACE OOOJ 





